matrix dimensions do not agree?

2 次查看（过去 30 天）
john birt
john birt 2012-4-26
how can I add more observed variables into a function 'f'
x = [1.4 1.5 1.6]; z=[4.1 5.3 6.7]; y = [0.2 0.31 0.43];
c = 0.1; %initial guess
t=0.1:0.1:1;
f = @(c,x) arrayfun(@(x)trapz(t,x.*c.*t.^2),x);
cfit = nlinfit(x,y,f,c)
works well, but trying to adapt the function to use more observed information in the form of the vector 'z' by coding
f = @(c,x) arrayfun(@(x)trapz(t,z+x.*c.*t.^2),x);
gives an error
??? Error using ==> nlinfit at 120
Error evaluating model function
'@(c,x)arrayfun(@(x)trapz(t,z+x.*c.*t.^2),x)'.
Error in ==> test88 at 6
cfit = nlinfit(x,y,f,c)
Caused by:
Error using ==> plus
Matrix dimensions must agree.

请先登录，再回答此问题。

采纳的回答

Andrei Bobrov
Andrei Bobrov 2012-4-26
x = [1.4 1.5 1.6]; z=[4.1 5.3 6.7]; y = [0.2 0.31 0.43];
c = 0.1; %initial guess
t=0.1:0.1:1;
f = @(c,x) arrayfun(@(x,z)trapz(t,x.*c.*t.^2),x,z);
cfit = nlinfit(x,y,f,c)
  2 个评论
john birt
john birt 2012-4-26
there is no 'z' in the trapz(t,x.*c.*t.^2) I was looking for trapz(t,z+x.*c.*t.^2)
john birt
john birt 2012-4-26
it works!!!
clear;
x = [1.4 1.5 1.6]; z=[4.1 5.3 6.7]; y = [0.2 0.31 0.43]; n=size(x);
c = 0.1; %initial guess
t=0.1:0.1:1;
f = @(c,x) arrayfun(@(x,z)trapz(t,z+x.*c.*t.^2),x,z);
[cfit, r] = nlinfit(x,y,f,c)
error = sum(r.^2)/n(2)

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by