数值计算方法上机2.docVIP

1．复化求积公式 设，取等距节点，，，作积分和式，则。 上，用作为的近似（即用零次插值函数作为的近似），用作为的近似，那么，可作为的近似，而且，当节点个数增加时，这种近似的效果会越来越好。这种现象，我们可以通过以下动画实验来直观地展示。 1】设，用动画摸拟。 1 编写函数文件f.m function y=f(x) y=exp(x); 步2 编写程序文件 clear all %清除命令空间中所有变量 %输入积分区间的左右端点、区间等分数 disp(用阶梯形面积逼近曲边梯形面积的数值试验与动画试验) a=input(请输入积分区间左端点a=); b=input(请输入积分区间右端点b=); n=input(请输入对区间等分的正整数数组n=[n1,n2,...](各分量是单增的)n=); index=max(size(n)); %绘制阶梯形逼近曲线梯形的动画,计算阶梯形面积 x=a:0.0005:b; y=f(x); for i=1:1:index plot(x,y,*r) %作被积函数图象 hold on sum(i)=0.0; step(i)=(b-a)/n(i); for j=1:1:n(i) u(j)=a+step(i)*(j-1); v(j)=f(u(j)); u(j+1)=a+step(i)*j; sum(i)=sum(i)+v(j)*step(i);%计算阶梯形面积 fill([u(j),u(j),u(j+1),u(j+1)],[0,v(j),v(j),0],b) %用蓝色填充阶梯形 end pause(2) %让图象显示暂停2秒 if iindex fill(x,y,w) %用白色擦去前面的阶梯形 end hold off end title(阶梯形逼近曲边梯形)%为图象加上标题 %输出区间等分数与阶梯形面积 table=[n,sum]; disp( 等分数 阶梯形面积) disp(table) 请输入积分区间左端点a=0 请输入积分区间右端点b=1 请输入对区间等分的正整数数组n=[n1,n2,...](各分量是单增的)n=[2,4,6,8,10] 梯形递推化公式的实现算法。 算法分析与设计： 1 的确定 这里，，， 2 的算法 引入变量，先置 一般地 引入变量，， ， ，，直到 3 计算 4 判断误差，确定输出或继续将步长两分 ，回到第2步。 具体算法： 1 2 3 4 5 6 7 8 作业：1.计算积分值）。 2．蒙特卡罗方法求积分 内容见书中P76 例题算法见P166 作业题：P79 5.4 3．用编程实现求解非线性方程求根 作业： P109 7.1（再加入牛顿迭代法求解） 补充内容：FFT与IFFT Xf=fft(Xt,N,DIM) 计算N点Xt序列的N点离散fourier变换Xf Xt=ifft(Xf,N,DIM) 计算N点Xf序列的N点离散fourier逆变换Xt 例3 利用fft和ifft指令重新计算两序列的卷积。所给已知序列为 ， 产生给定序列a（n），b（n） clear a=ones(1,13); a([1,2,3])=0; b=ones(1,10); b([1,2])=0 直接卷积 c=conv(a,b); 通过变换求卷积 M=32； AF=fft(a,M); BF=fft(b,M); CF=AF.*BF; %必须采用点乘 cc=real(ifft(CF)); %过滤因截断误差所引起的虚部 图示计算结果 nn=0:(M-1); %保证绘图时，时间从0时刻开始 c(M)=0; %使直接卷积所得序列通过尾部补0，与cc长度相同 error=c-cc; subplot(2,1,1),stem(nn,c,’fill’),grid,axis([0,31,0,9]) xlabel(‘nn’),ylabel(‘cc’); subplot(2,1,2),stem(nn,error,’fill’),axis([0,31,-1,1]) ylabel(‘error’)