| 查看: 357 | 回复: 0 | ||
[求助]
求助！！！MATLAB进行多元分线性分析，小弟初次接触matlab，还请各位多多指教，谢谢！
|
|
求助！！！MATLAB进行多元分线性分析，自编的程序问题在哪，，小弟初次接触matlab，还请各位多多指教，谢谢！！！ function yy=myfun(beta,x) a=beta(1); %待定系数 b=beta(2); c=beta(3); x1=x(:,1); %每一列 x2=x(:,2); x3=x(:,3); yy=beta(1)*((2*x1.*(1-sqrt(1-x2))-9*x1.^2)-(x1.^2-(x1-0.003*(beta(2)+beta(3)*x3./x1))^2)); %自定义函数 x=[8 0.1163 15;8 0.0786 15;8 0.0762 15;16 0.0366 15;16 0.0285 15;16 0.0571 15;20 0.0295 15;20 0.0218 15;20 0.0178 15;20 0.0205 20; 20 0.0182 20;20 0.0286 20;8 0.113 20;8 0.113 20;8 0.1316 20;16 0.0299 20;16 0.0273 20;16 0.0335 20;16 0.0458 25; 16 0.0461 25;16 0.0263 25;20 0.0295 25;20 0.0224 25;20 0.0388 25;8 0.1029 25;8 0.1453 25;8 0.126 15]; y=[0.17 0.21 0.23 0.22 0.18 0.15 0.2 0.19 0.24 0.15 0.17 0.21 0.15 0.19 0.23 0.22 0.24 0.18 0.2 0.15 0.23 0.22 0.17 0.15 0.15 0.19 0.18]'; beta0=[0.06 7.6 9]; %对待定系数赋初值 [beta,r,J]=nlinfit(X,y,'myfun',beta0); plot(x,y,'+'); hold on; xx=2:0.05:25; yy=model(beta,xx);%用户自编函数 plot(xx,yy,'r') |
回复此楼» 猜你喜欢
2026国自然放榜时间
已经有6人回复
-
据悉今年马上要出结果了
已经有3人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有7人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有15人回复
-
面上项目filecode邪修
已经有7人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
100