I am receiving this error "Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is 5-by-5. Error in fminsearch (line 200) "

I found this code on MATHWORKS site, which I am using find the order of a chemical reaction. But for some reason the code is not working. I am using MATLAB 2020 while the code was coded in MATLAB 2010. Would be much appreciated if someone would help me out.

[KSSV]

t=[0 1 2] ;
C=[20 40 70] ;
reorder(t,C)

Incorrect number or types of inputs or outputs for function 'reorder'.

[Aynkaran]

thanks for the feedback. But still I cannot figure out the error.

[KSSV]

Running reorder gives different error. Are you talking about the same function or it is a function written by the user?

[Torsten]

@Aynkaran

You must show us your function "rorder" for that we are able to help.

[Aynkaran]

[Torsten]

You see the problem ?

t = [1 2 3 4 5];
C = [200 134 120 110 89];
corr(t,C)

ans = 5×5

NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN

t = t.';
C = C.';
corr(t,C)

ans = -0.9235

[Aynkaran]

t=[1 2 3 4 5];

C=[200 134 120 110 89];

dispf (rorder(t,C))

function data = rorder(t,C)

%DATA = ORDERSEEK(T,C)

%Searches for the least-squares optimum solution for the reaction order

%from inputs time (t) and concentration (C). Outputs a struct (data) which

%contains the (1) reaction order (2) rate constant (3) rsquared value (4)

%half-life.

%t and C are vectors.

n = fminsearch(@(x) os_obj(t,C,x),0);

best_fit = polyfit(t,1./(C.^(n-1)),1);

k = best_fit(1)/(n-1);

rsquared = os_obj(t,C,n);

halflife = (2^(n-1) - 1)/...

((n-1)*k*C(1)^(n-1));

data.order = n;

data.k = k;

data.rsquared = -rsquared; %Revert back to positive value.

data.halflife = halflife;

end

function rsquared = os_obj(t,C,n)

%Order seek objective function. Returns the correlation coefficient for the

%line fit of the transformed data.

C_xform = 1./(C.^(n-1));

rsquared = -corr(t,C_xform).^2; %Negative, because we want the minimizer to maximize it.

%(c) February 8, 2010, Bradley James Ridder, University of South Florida.

%References:

%http://en.wikipedia.org/wiki/Rate_equation#Summary_for_reaction_orders_0.2C_1.2C_2_and_n

end

t=[1 2 3 4 5];
C=[200 134 120 110 89];
dispf (rorder(t,C))
function data = rorder(t,C)
%DATA = ORDERSEEK(T,C)
%Searches for the least-squares optimum solution for the reaction order
%from inputs time (t) and concentration (C). Outputs a struct (data) which
%contains the (1) reaction order (2) rate constant (3) rsquared value (4)
%half-life.
%t and C are vectors.
n = fminsearch(@(x) os_obj(t,C,x),0);
best_fit = polyfit(t,1./(C.^(n-1)),1);
k = best_fit(1)/(n-1);
rsquared = os_obj(t,C,n);
halflife = (2^(n-1) - 1)/...
           ((n-1)*k*C(1)^(n-1));
data.order = n;
data.k = k;
data.rsquared = -rsquared; %Revert back to positive value.
data.halflife = halflife;
end

function rsquared = os_obj(t,C,n)
%Order seek objective function. Returns the correlation coefficient for the
%line fit of the transformed data.
C_xform = 1./(C.^(n-1));
rsquared = -corr(t,C_xform).^2; %Negative, because we want the minimizer to maximize it.
%(c) February 8, 2010, Bradley James Ridder, University of South Florida.
%References:
%http://en.wikipedia.org/wiki/Rate_equation#Summary_for_reaction_orders_0.2C_1.2C_2_and_n
end 

[Aynkaran]

Please help me out.

[Torsten]

You have to transpose t and C !!

[Aynkaran]

Great!!!! Thanks Sir

[Aynkaran]

How can I rate you with stars?? is there an option like that in this platform??

[Torsten]

I turned my comment to an answer.

If you like, you can accept it. Thanks.

Answer 1

[Torsten]

You have to transpose t and C !!