MATLAB程序设计与应用实验报告5 下载本文

内容发布更新时间 : 2024/6/9 0:54:16星期一 下面是文章的全部内容请认真阅读。

-0.2725 0.2820 -0.2857 -0.1332 0.3075 0.7015 0.3502 -0.0200 -0.0348 0.7914 1.0984 0.5080 0.5201 1.0187 1.6555 1.3790 1.1275 0.9642 1.3514 2.5260 ans =

1.3790 -0.0348 -1.1201 -1.5771 -1.7502 2.5260 0.5080 -0.2857 -0.7982 -1.0582 1.6555 1.0187 0.2820 -0.4686 -0.8314 0.3075 0.0335 -0.1332 -0.2725 -0.9792 1.0984 -0.7145 -1.1564 -1.2571 -1.3337 1.3514 1.1275 -0.2779 -0.5336 -0.8655 0.7015 0.3502 -0.1765 -0.2248 -2.0026 0.9642 0.7914 -0.2991 -0.5890 -2.0518 0.5201 0.0229 -0.2938 -0.3538 -1.3320 -0.0200 -0.2620 -0.8236 -0.8479 -2.3299

2. 按下表用3次样条插值方法计算0~90°范围内整点数的正弦值和0~75°范围内整数点的正切值,然后用5次多项式拟合方法计算相同的函数值,并将两种计算结果进行比较。

表1 特殊角的正弦与正切值表

α(度) sinα tanα 0 0 0 15 0.2588 0.2679 30 0.5000 0.5774 45 0.7071 1.0000 60 0.8660 1.7320 75 0.9659 3.7320 90 1.0000 a=[0 15 30 45 60 75 90]; A=0:90;

sing=[0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000]; b=[0 15 30 45 60 75];

B=0:75;

tang=[0 0.2679 0.5774 1.0000 1.7320 3.7320]; SC=interp1(a,sing,A,'spline'); TC=interp1(b,tang,B,'spline'); P1=polyfit(a,sing,5); SN=polyval(P1,A); P2=polyfit(b,tang,5); TN=polyval(P2,B);

plot(A,SC,'.',A,SN,'r-',B,TC,'.',B,TN,'b-')

223. 已知多项式P2(x)?5x?x?2,P3(x)?x?0.5,求: 1(x)?3x?2,P(1)P(x)?P1(x)P2(x)P3(x) (2)P(x)?0的全部根。 (3)计算xi?0.2i(i?0,1,2,P1=[3 2]; P2=[5 -1 2];

,10)各点上的P(xi)。

P3=[1 0 -0.5]; P1=[0 P1]; P12=conv(P1,P2); P3=[0 P3]; P=conv(P12,P3) X=roots(P) xi=0.2*(0:10); PZ=polyval(P,xi) P =

1 至 5 列

0 0 15.0000 7.0000 -3.5000

6 至 8 列

0.5000 -2.0000 -2.0000 X =

0.7071 + 0.0000i 0.1000 + 0.6245i 0.1000 - 0.6245i -0.7071 + 0.0000i -0.6667 + 0.0000i PZ =

1 至 5 列

-2.0000 -2.3920 -2.6112 -1.7024 2.7104

6 至 10 列

15.0000 42.1120 94.1408 184.9056 332.5264 11 列

560.0000

n?1?f1?1?f?2n?2?24. 已知?,求f1~f100中:

f?0n?3?3?n?3?fn?2fn?1?fn?2?fn?3(1) 最小数与总和。 (2) 中值与均值。 f(1)=1; f(2)=2; f(3)=0; for n=4:100;

f(n)=2*f(n-1)-f(n-2)+f(n-3); end min=min(f) sum=sum(f) median=median(f) mean=mean(f) min =