计算三角多项式.docVIP

计算三角多项式的MATLAB主程序 (串行) 输入的量： 输出的量：A1,B1,Y11,Rm1 function [A,B,Y1,Rm]=sanjiao(X,Y,X1,m) n= length(X)-1;max1=fix((n-1)/2); if m max1 m=max1; end A=zeros(1,m+1);B=zeros(1,m+1); Ym=(Y(1)+Y(n+1))/2; Y(1)=Ym; Y(n+1)=Ym; A(1)=2*sum(Y)/n; for i=1:m B(i+1)=sin(i*X)*Y; A(i+1)=cos(i*X)*Y; end A=2*A/n; B=2*B/n; A(1)=A(1)/2;Y1=A(1); for k=1:m Y1=Y1+A(k+1)*cos(k*X1)+ B(k+1)*sin(k*X1); Tm=A(1)+A(k+1).*cos(k*X)+ B(k+1).*sin(k*X); k=k+1; end Y;Tm; Rm=(sum(Y-Tm).^2)/n; 计算三角多项式的MATLAB主程序 (并行) 输入的量： 输出的量：A1,B1,Y11,Rm1 并行代码function [A,B,Y1,Rm]=sanjiao_spmd(X,Y,X1,m,p) n= length(X)-1;max1=fix((n-1)/2); if m max1 m=max1; end A=zeros(1,m+1);B=zeros(1,m+1); Ym=(Y(1)+Y(n+1))/2; Y(1)=Ym; Y(n+1)=Ym; A(1)=2*sum(Y)/n; fix_m=fix(m/p); spmd(p) max_size=fix_m*labindex; min_size=(labindex-1)*fix_m+1; if labindex==p for i=min_size:m B(i+1)=sin(i*X)*Y; A(i+1)=cos(i*X)*Y; end else for i=min_size:max_size B(i+1)=sin(i*X)*Y; A(i+1)=cos(i*X)*Y; end end end if p==4 A=A{1}+A{2}+A{3}+A{4}; B=B{1}+B{2}+B{3}+B{4}; elseif p==3 A=A{1}+A{2}+A{3}; B=B{1}+B{2}+B{3}; elseif p==2 A=A{1}+A{2}; B=B{1}+B{2}; else A=A{1}; B=B{1}; end A(1)=2*sum(Y)/n; A=2*A/n; B=2*B/n; A(1)=A(1)/2;Y1=A(1); for k=1:m Y1=Y1+A(k+1)*cos(k*X1)+ B(k+1)*sin(k*X1); k=k+1; end spmd(p) max_size=fix_m*labindex; min_size=(labindex-1)*fix_m+1; if labindex==p for k=min_size:m TM=A(k+1).*cos(k*X)+ B(k+1).*sin(k*X); k=k+1; end else for k=min_size:max_size TM=A(k+1).*cos(k*X)+ B(k+1).*sin(k*X); k=k+1; end end end for k=1:p Tm=A(1)+TM{k}; end Y;Tm; Rm=(sum(Y-Tm).^2)/n; 注意：下例子仅用于测试程序正确性与并行时的加速比，故没有使用例子中全部测试代码。 例7.7.1 根据上的个等距横坐标点 和函数. （1）求的6阶三角多项式逼近，计算均方误差; （2）将这三个三角多项式分别与的傅里叶级数 的前6项进行比较; （3）利用三角多项式分别计算Xi= -2, 2.5的值; （4）在同一坐标系中，画出函数，的三角多项式和数据点的图形. 解 （1）输入程序 X1=-pi:2*pi/13:pi;Y1=2*sin(X1/3);X1i=[-2,2.5]; [A1,B1,Y11,R