| 查看: 321 | 回复: 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') |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
博士读完未来一定会好吗
已经有27人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
