《科学计算与MATLAB》期末大作业.doc

杭州电子科技大学信息工程学院 《科学计算与MATLAB》期末大作业 题 目 专 业 学 期 学 号 姓 名 任课教师 完成日期 2016年12月31日 给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m的MATLAB函数，用以计算下述的值： 绘制t关于函数f(t)的图形，其中t的取值范围为，间距为。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)sin(pi/4))+(sin(-pi/4)).*(sin(t)sin(-pi/4))+(sin(t)).*((sin(t)=sin(pi/4))(sin(t)=sin(-pi/4)));%函数表示 plot(t,y); %画图 end 2.解以下线性方程组 A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵 求： (1)A(2:3,2:3) A(:,1:2) A(2:3,[1,3]) [A,[ones(2,2)；eye(2)]] A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵 A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵 数学函数定义在区域[-8,8]×[-8,8]上。建立一函数文件，调用该函数时传递参数k,当k等于1时画函数z的三维网格图，当k等于2时画函数z的三维曲面图，当k等于3时画函数z的三维等高线图。 function z=T4(k)%定义函数 if k==1 %判断形参 x=-8:0.01:8; %定义变量范围 y=-8:0.01:8; %定义变量范围 z = (sin((x.^2+y.^2).^0.5))./((x.^2+y.^2).^0.5);%函数表示 plot3(x,y,z) %画图 grid on %画网格 xlabel(x) %x轴名字 ylabel(y) %y轴名字 zlabel(z) %z轴名字 elseif k==2 x=-8:0.01:8; %定义变量范围 y=x; %y赋值 [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;%函数表示 Z = sin(R)./R;%函数表示 mesh(X,Y,Z); %画图 xlabel(x) %x轴名字 ylabel(y) %y轴名字 zlabel(z) %z轴名字 elseif k==3 x=-8:0.5:8; %定义变量范围 y=x; %y赋值 [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;%函数表示 Z = sin(R)./R;%函数表示 surfc(X,Y,Z); %画图 colorbar xlabel(x) %x轴名字 ylabel(y) %y轴名字 zlabel(z) %z轴名字 end end 三维网格图 三维曲面图 三维等高线图 5.按以下要求分别用不同颜色、不同线型绘出sinx、cosx、0.003ex三条函数曲线，并对图名、坐标名及图例进行标注。 （1）分别在三个窗口内画出三条函数曲线。 在同一窗口的不同坐标系下画出三条函数曲线； 题（1） x=0:pi/20:4*pi; %x范围 a=sin(x); %函数表示 plot(x,a,ro) %画线（红，o） figure %生成figure框 b=cos(x); %函数表示 plot(x,b,gx) %画线（绿，x） figure %生成figure框 c=0.