实验程序40.docxVIP

MATLAB数学实验程序ex1. 单位圆的正多边形逼近alpha=linspace(0,2*pi,13);x=sin(alpha); y=cos(alpha); %计算正十二边形顶点坐标I2=1:2:13;I1=1:4:13;x1=x(I1);y1=y(I1); %提取正三边形坐标x2=x(I2);y2=y(I2); %提取正六边形坐标plot(x1,y1,x2,y2,x,y)axis offex2 抛射曲线的最大射程计算v0=100; %设定初速度g=9.8;alpha=pi/4;T=2*v0*sin(alpha)/g; %计算飞行时间t= linspace(0,T,17); %设定飞行时刻x=v0*t*cos(alpha); %计算坐标点y=v0*t*sin(alpha)-g*t.^2/2;plot(x,y,x,y,r*)xmax=x(end) %显示最远射程ex3. 巴拿马草帽.图形[x,y]=meshgrid(-8:0.5:8); %创建网格点矩阵r=sqrt(x.^2+y.^2)+eps; %计算网格点到原点距离z=sin(r)./r; %计算二元函数值mesh(x,y,z) %二元函数绘图colormap([1,0,0]) %设置红色axis off %去掉坐标架ex4.三叶玫瑰线实验theta=linspace(0,2*pi,10000); %创建区间[0，2]上离散点n=3; %设置三叶玫瑰线参数r=cos(n*theta); %计算玫瑰线坐标数据figure(1) %创建第一个图形窗polar(theta,r,k) %极坐标绘图x=r.*cos(theta); %直角坐标转换y=r.*sin(theta);figure(2) %创建第二个图形窗口comet(x,y) %用慧星绘图做动态模拟ex5.抛射曲线簇绘图程序alpha=linspace(0.1,pi/2,20); %创建20个发射角数据的行向量v0=100;g=9.8;Taim=2*v0*sin(alpha)/g; %计算对应的飞行时间T=(0:16)*Taim/16; %计算飞行时刻矩阵TX=v0*T*diag(cos(alpha)); %对角矩阵右乘飞行时刻矩阵Y=v0*T*diag(sin(alpha))-g*T.^2/2;plot(X,Y,r)ex6. 黎曼曲面绘图程序[r,theta]=meshgrid(s,t); %创建网格点z=r.*exp(i*theta); %设置复数自变量f=sqrt(r).*exp(i*theta/2); %计算复变函数值x=real(z);y=imag(z); %分离实部和虚部figure(1),mesh(x,y,real(f)),axis offfigure(2),mesh(x,y,imag(f)),axis offexp7. 宝石切割问题实验程序P=perms(1:3); %创建全排列操作数a=[19,14,12];b=[5,4,2]; %设置石料数据for ki=1:6pk=P(ki,:);L0=a; %提取操作数 for kj=1:3 L=L0;II=pk(kj);L(II)=[]; %模拟切割 A(ki,kj)=prod(L); %计算单次切割面积 L0(II)=b(II); endendQ=A*[1;1;1]; %计算每种操作切割面积di