MATLAB践作品.docVIP

MATLAB实践 材料1209班 0607120901 杨建华 第一题： （1） M文件部分： % diyiti.m fid=fopen(d:/Data.txt); head=fscanf(fid,%s,40); for i=1:19 ch=fscanf(fid,%f,1); d=fscanf(fid,%f,[19,19]); data(:,:,i)=d; [phi1,phi,phi2]=meshgrid(0:5:90,0:5:90,0:5:90); slice(phi1,phi,phi2,data,[45,90],[45,90],[0,45]) 程序部分： diyiti （2） M文件部分： % dierti.m fid=fopen(d:/Data.txt); head=fscanf(fid,%s,40); [phi1,phi]=meshgrid(0:5:90,0:5:90); for i=1:19 ch=fscanf(fid,%f,1); da=fscanf(fid,%f,[19,19]); subplot(5,4,i) contour(phi1,phi,da) %同理注意转置 end 程序部分： dierti （3） M文件部分： disanti.m fid=fopen(d:/Data.txt); head=fscanf(fid,%s,40); for i=1:19 d=fscanf(fid,%f,[19,19]); data(:,:,i)=d; end b=data(10,:,1) %注意此处是第一个数为10而不是第二个 a=0:5:90; plot(a,b,-o) 程序部分： disanti 开始 开始 载入Date 载入Date.txt 读数据ch 读数据 ch为单个φ2数据 d为每个φ2对应的19*19数据 画坐标图利用subplot打开五行四列的图形窗口meshgrid建三维矩阵 画坐标图 利用subplot打开五行四列的图形窗口 meshgrid建三维矩阵 画层片图画切片图 画层片图 画切片图 ·第二题 编程实现以下科学计算算法，并举一例应用之。（参考书籍《ＭＡＬＡＢ与科学计算》，王沫然著，电子工业出版社，２００９年） “Lagrange插值” M文件部分： % lagrange insert function y=lagrange(x0,y0,x); n=length(x0); m=length(x); for i=1:m r=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(r-x0(j))/(x0(k)-x0(j)); end end s=s+p*y0(k); end y(i)=s; end 程序部分： 举例：求sin（pi/13） x0=[0:pi/20:2*pi]; y0=sin(x0); lagrange(x0,y0,pi/13) ans = 0.2393 ２、编程解决以下科学计算和工程实际问题。 M文件部分： % sumF.m clear N=input(输入力的数目N=); for i=1:N f(i,:)=input(输入了的两个分量[fi(x),fi(y)]=); end F=sum(f) 程序部分： sumF 输入力的数目N=2 输入力f(i)的两个分量[fi(x),fi(y)]=[1,3] 输入力f(i)的两个分量[fi(x),fi(y)]=[4,6] F = 5 9 开始 开始 输入力的数目N 输入力的数目N i=1 i=1 输入力的分量[fi(x),fi(y)] 输入力的分量[fi(x),fi(y)] i=i+1 i=i+1 i=N? i=N? 否是 否 是 F=sum F=sum(f) 2) （1） syms x; p1=sym2poly(x^2+x+1); r11=roots(p1) r11 = -0.5000 + 0.8660i -0.5000 - 0.8660i polyval(p1,r11) ans = 1.0e-15 * %误差小