数值分析相关实验代码.docx

下载文档

7
0
约1.53千字
约 5页
2018-05-17 发布于河南
举报
版权申诉
保障服务

数值分析相关实验代码.docx

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

数值分析相关实验代码

连续复合梯形和辛普森公式：function y=hanshu(x)y=(2*exp(-x))/sqrt(pi);………………………………x=0:0.01:1;h=0.01;futix=0;for n=1:100 y1=hanshu(x(n)); y2=hanshu(x(n+1));futix=futix+(y1+y2)*h/2;endfutix=vpa(futix,6)fusim=0;for n=1:100 y1=hanshu(x(n)); y2=hanshu(x(n+1)); y3=hanshu((x(n)+x(n+1))/2);fusim=fusim+(y1+y2+4*y3)*h/6;endfusim=vpa(fusim,6)离散复合梯形和辛普森公式：x=[0,1,4,9,16,25,36,49,64];y=[0,1,2,3,4,5,6,7,8];plot(x,y,ro)x1=0:0.1:64;y1=interp1(x,y,x1,spline);hold onplot(x1,y1,b.)hold offftix=0;h=0.1;for n=1:length(x1)-1ftix=ftix+(y1(n)+y1(n+1))*h/2;endftix=vpa(ftix,10)fsim=0;h=h*2;for n=1:2:length(x1)-2fsim=fsim+(y1(n)+y1(n+2)+4*y1(n+1))*h/6;endfsim=vpa(fsim,10)改进的欧拉方法：functiondy=hanshu(x,y)dy=1/x^2-y/x;………………………………x=1:0.05:2;y(1)=1;h=0.05;for n=1:length(x)-1 y1=y(n)+h*hanshu(x(n),y(n)); y(n+1)=y(n)+h*(hanshu(x(n),y(n)))+hanshu((x(n)+h),y1)/2;endy龙格库塔方法：functiondy=h(t,y)dy=zeros(3,1);dy(1)=-y(1);dy(2)=y(1)-10*y(2);dy(3)=-10*y(2);…………………………tspan=0:0.01:1;[T,Y] = ode45(@h,tspan,[1 1 1]);plot(T,Y(:,1),-,T,Y(:,2),-.,T,Y(:,3),.)LU分解：A=[1,2,3;2,5,2;3,1,5];b=[14;18;20];ifdet(A)==0error();else [m,n]=size(A); L=zeros(m,m); L=L+diag(ones(m,1)); U=zeros(m,m);U(1,:)=A(1,:);L(2:end,1)=A(2:end,1)/U(1,1);for r=2:m U(r,r)=A(r,r)-L(r,1:r-1)*U(1:r-1,r);for i=r+1:m U(r,i)=A(r,i)-L(r,1:r-1)*U(1:r-1,i); L(i,r)=(A(i,r)-L(i,1:r-1)*U(1:r-1,r))/U(r,r);endendendLUy(1)=b(1);for i=2:ny(i)=b(i)-L(i,1:i-1)*y(1:i-1);endx(n)=y(n)/U(n,n);for i=n-1:-1:1 x(i)=(y(i)-U(i,i+1:n)*x(i+1:n))/U(i,i);endx