Question

matlab非线性拟合问题

h=[2632.23 2624.30 2620.67 2610.29 2606.61 2599.59 2587.60 2582.05 2579.57 2575.44 2569.46 2564.12 2559.83 2548.47 2539.63 2528.01 2521.63 2510.23 2508.17 2500.07 ];
y=[60448.88 60311.43 60248.03 60065.11 59999.69 59874.06 59657.02 59555.51 59509.94 59433.77 59322.85 59223.17 59142.66 58927.69 58758.61 58534.01 58409.58 58185.31 58144.52 57983.36 ];
拟合方程y=(1/3)*((150^2*3.1416-150^2*(1-(h+220*c(1)-150)*(1-...
c(2))/150)+(h+200*c(1)-150)*(1-c(2))*sqrt(150^2-((h+200*c(1)-150)^2)*((1-...
c(2))^2)))*(h/c(1)+200)-(150^2*3.1416-150^2*(1-(h-600*c(1)-150)*(1-...
c(2))/150)+(h-600*c(1)-150)*(1-c(2))*sqrt(150^2-((h-600*c(1)-150)^2)*((1-...
c(2))^2)))*(h/c(1)-600))...
+3.1416/2*(150^2*(h+200*c(1)-150)-((h+200*c(1)-150)^3)/3+2/3*150^3)...
+3.1416/2*(150^2*(h-600*c(1)-150)-((h-600*c(1)-150)^3)/3+2/3*150^3);
拟合参数 c(1) c(2)
怎样通过这些数据 求出 c(1) c(2)的值呢？
需要matlab 编码 望高手支招~~~
一楼的代码用后报错??? Error using ==> mpower
Matrix must be square.

Error in ==> myfun at 2
y=(1/3)*((150^2*3.1416-150^2*(1-(h+220*c(1)-150)*(1-...

Error in ==> lsqcurvefit at 209
initVals.F =
feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});

Caused by:
Failure in initial
user-supplied objective
function evaluation.
LSQCURVEFIT cannot
continue.

是怎么回事呢..?

Answer 1

h=[2632.23 2624.30 2620.67 2610.29 2606.61 2599.59 2587.60 2582.05 2579.57 2575.44 2569.46 2564.12 2559.83 2548.47 2539.63 2528.01 2521.63 2510.23 2508.17 2500.07 ];
y=[60448.88 60311.43 60248.03 60065.11 59999.69 59874.06 59657.02 59555.51 59509.94 59433.77 59322.85 59223.17 59142.66 58927.69 58758.61 58534.01 58409.58 58185.31 58144.52 57983.36 ];
F=@(c,h)(1/3)*((150^2*3.1416-150^2*(1-(h+220*c(1)-150)*(1-c(2))/150)+(h+200*c(1)-150)*(1-c(2))*sqrt(150^2-((h+200*c(1)-150)^2)*((1-c(2)）^2)))*(h/c(1)+200)-(150^2*3.1416-150^2*(1-(h-600*c(1)-150)*(1-
c(2))/150)+(h-600*c(1)-150)*(1-c(2))*sqrt(150^2-((h-600*c(1)-150)^2)*((1-c(2))^2)))*(h/c(1)-600))+3.1416/2*(150^2*(h+200*c(1)-150)-((h+200*c(1)-150)^3)/3+2/3*150^3)+3.1416/2*(150^2*(h-600*c(1)-150)-((h-600*c(1)-150)^3)/3+2/3*150^3);
c=lsqcurvefit(F,[1 1],h,y)
应该是向量乘除的问题，应该在乘除后面加点号，例如：h.^2,h./3,h.*(h+1)等等，楼主自己改一下F函数。