华中科技大学数值分析 数值分析实验程序.doc

《华中科技大学数值分析 数值分析实验程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.....................最新资料整理推荐.....................Functiont_charpt1_1%数值试验1.1病态问题%输入：[020]之间的扰动项及小的扰动常数%输出：加扰动后得到的全部根Result=inputdlg({‘请输入扰动项：在[020]之间的整数:’},’charpt1_1’,1’{‘19’});Numb=str2num(char(result));if((Numb>20)

2、(Numb<0))errordlg(‘请输入正确的扰动项：[020]之间的整数!’);return;endre

3、sult=inputdlg({‘请输入(01)之间的扰动常数：’},’charpt1_1’,1,{‘0.00001’});ess=str2num(char(result));ve=zeros(1,21);ve(21-Numb)=ess;root=roots(poly(1:20)+ve);disp([‘对扰动项’,num2str(Numb),’加扰动’,num2str(ess),’得到的全部根为:’]);disp(num2str(root));8.....................最新资料整理推荐.....................

4、functioncharpt3%数值实验三：含“实验3.1”和“实验3.2”%子函数调用:dlsa%输入：实验选择%输出：原函数及求得的相应插值多项式的函数的图像以及参数alph和误差rresult=inputdlg({‘请选择实验，若选3.1，请输入1，否则输入2:’},’charpt_3’,1,{‘1’});Nb=str2num(char(result));if(Nb~=1&(Nb~=2)errordlg(‘实验选择错误！’))；return;endx0=-1:0.5:2y0=[-4.447-0.4520.5510.048-0.4470.

5、5494.552];n=3;%n为拟合阶次if(Nb==1)alph=polyfit(x0,y0,n);y=polyval(alph,x0);r=(y0-y)*(y0-y)’;%平方误差x=-1:0.01:2;y=polyval(alph,x);plot(x,y,’k-’);8.....................最新资料整理推荐.....................xlabel(‘x’);ylabel(‘y0*andployfit.y-’);holdon;plot(x,y,’k-’);title(‘离散数据的多项式拟合’);grid

6、on;elseresult=inpurdlg({‘请输入权向量w:’},’charpt_3,1,{‘[1111111]’});w=str2num(char(result));[a,b,c,alph,r]=dlsa(x0,y0,w,n);enddisp([‘平方误差:’,sprint(‘%g,r’)]);disp([‘参数alph:’,sprint(‘%t’,alph)])%-------------------------------------------------------------------------------------

7、-----------------------------------------function[a,b,c,alph,r]=dlsa(x,y,w,n)%功能：用正交化方法对离散数据作多项式最小二乘拟合。%输入：m+1个离散点（x,y,w）,x,y,w分别用行向量给出。%拟合多项式的次数n,0

8、...........if(n<1

9、n>=m)errordlg(‘错误：n<1或者n>=m!’);return;end%求三项递推公式的参数a,b,拟合多项式s(x)的系数c,其中d(k)=(y,sk);s1=0;s2=ones(1,m+1);v2=sum(w);d(1)=y*w’;c(1)=d(1)/v2;fork=1:nxs=x.*s2.^2*w’;a(k)=xs/v2;if(k==1)b(k)=0;elseb(k)=v2/v1;ends3=(x-a(k)).*s2-b(k)*s1;v3=s3.^2*w’;d(k+1)=y.*s3*w’;

10、c(k+1)=d(k+1)/v3end%求平方误差rr=y.*y*w’-c*d’;%,求拟合多项式s(x)的降幂系数alph8...................