24小时热门版块排行榜

返回列表

yjstudy

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 90
帖子: 3
在线: 10小时
虫号: 4007648
注册: 2015-08-05
专业: 结构工程

[求助] 求助！！！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')

回复此楼

» 猜你喜欢

1楼 2015-12-16 19:34:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 yjstudy 的主题更新

返回列表