- 积分
- 0
- 注册时间
- 2010-4-9
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2011-12-3 13:23:11
|
显示全部楼层
来自 四川成都
guojunhang 发表于 2011-12-3 13:00 ![](static/image/common/back.gif)
不一样,具体原因可以参考《数值方法》。
令附:
从x=[1:3],一共有3个点,最多可拟合出2次多项式,哪来的5 ...
多谢指导~~
x=[0:0.5:10];
y=[0.008368214
0.007191997
0.006001203
0.004942747
0.004060792
0.003399394
0.00279685
0.002293626
0.001823075
0.001484878
0.001220212
0.000955543
0.000705578
0.000485028
0.000279182
0.000146868
-7.36809E-05
-0.000279526
-0.000147135
-7.35691E-05
0]';
plot(y,x,'y');
x1=x(:,1:5);
y1=y(:,1:5);
x2=x(:,5:10);
y2=y(:,5:10);
x3=x(:,10:15);
y3=y(:,10:15);
x4=x(:,15:20);
y4=y(:,15:20);
p1=polyfit(x1,y1,5);
p2=polyfit(x2,y2,5);
p3=polyfit(x3,y3,5);
p4=polyfit(x4,y4,5);
hold on
y1=polyval(p1,x1);
y2=polyval(p2,x2);
y3=polyval(p3,x3);
y4=polyval(p4,x4);
x=[x1,x2,x3,x4];
y=[y1,y2,y3,y4];
set(gca,'YDir','reverse');
plot(y,x,'r')
plot(y,x,'r',y1,x1,'s',y2,x2,'d',y3,x3,'r*',y4,x4,'-+')
帮忙看看这个分段拟合有没有问题?谢谢 |
|