数学建模实验——举重模型的matlab实现

模型一：线性模型代码如下：W=[56,62,69,77,85,94,105];Y=[305,327,358,380,394,418,436];p=polyfit(W,Y,1);Y1=polyval(p,W);plot(W,Y,'rx');hold on;plot(W,Y1,'b');legend('原数据','模型一')各个参数：运行结果：模型二：幂函数模型代...

@玉面小蛟龙

8066人浏览 · 2020-05-07 08:51:38

@玉面小蛟龙 · 2020-05-07 08:51:38 发布

模型一： 线性模型
代码如下：

W=[56,62,69,77,85,94,105];
Y=[305,327,358,380,394,418,436];
p=polyfit(W,Y,1);
Y1=polyval(p,W);
plot(W,Y,'rx');
hold on;
plot(W,Y1,'b');
legend('原数据','模型一')

各个参数：
在这里插入图片描述

运行结果：
在这里插入图片描述
模型二： 幂函数模型
代码如下：

W=[56,62,69,77,85,94,105];
Y=[305,327,358,380,394,418,436];
f1=@(X,W)X*W.^(2/3);
[X,cansha]=lsqcurvefit(f1,1,W,Y);
Y2=f1(X,W);
plot(W,Y,'rx');
hold on;
plot(W,Y2,'b');
legend('原数据','模型二');

各个参数：
在这里插入图片描述

运行结果
在这里插入图片描述

模型三： 幂函数改进模型
代码如下：

W=[56,62,69,77,85,94,105];
Y=[305,327,358,380,394,418,436];
a=[1,2];
f2=@(a,W)a(1)*(W-a(2)).^(1/3);
[a,cansha2]=lsqcurvefit(f2,a,W,Y);
Y3=f2(a,W);
plot(W,Y,'rx');
hold on;
plot(W,Y3,'m');
legend('原数据','模型三');