《数学实验》实验报告——用MATLAB绘制sin曲线

课程名称数学实验成绩评定实验项目名称曲线绘制【实验目的】1．了解曲线的几种表示方式。

2．学习、掌握MA TLAB软件有关的命令。

【实验内容】绘制下列四种曲线：1．以直角坐标方程y=sin x,y=cos x表示的正、余弦曲线。

2．以参数方程x=cos t,y=sin t,t∈[0,2π]表示的平面曲线（单位圆）。

3．以参数方程x=e−0.2t cosπ2t,y=π2e−0.2t sin t,z=t,t∈[0,20]表示的空间曲线。

4．作出摆线的图形。

5.做出以参数方程x=e−0.25t cosπ2t,y=e−0.25t sinπ2t,z=t,t∈[0,30]表示的空间曲线。

6．以极坐标方程r=a(1+cosϕ),a=1,ϕ∈[0,2π]表示的心脏线。

7．绘制极坐标系下曲线 ρ=acos (b+nθ)的图形，讨论参数a、b和n对其图形的影响。

8．（曲线族绘制）三次抛物线的方程为y=ax3+cx，讨论参数a和c对其图形的影响。

【实验方法与步骤】练习1做出函数y=sin x,y=cos x的图形，并观察它们的周期性。

MATLAB代码及结果如下：>> x=0:0.01*pi:4*pi;y1=sin(x);y2=cos(x);plot(x,y1,'b',x,y2,'r');legend('y=sin(x)','y=cos(x)','location','best');axis([0 4*pi -1 1])绘制结果如下图：y=sin x,y=cos x的图形如上图，两个函数的周期皆为2π练习2设y=√32e−4t sin(4√3t+π3)，要求以0.01秒为间隔，求出y的151个点，绘出y及其导数的图形。

MATLAB代码及结果如下：dt=0.01;t=0:0.01:1.5;w=4*sqrt(3); %设定频率y=sqrt(3)/2*exp(-4*t).*sin(w*t+pi/3);Dy=diff(y)/dt; %求导for i =1:length(t)-1t1(i)=t(i);endsubplot(2,1,1);plot(t,y);xlabel('时间t');ylabel('y(t)');gridsubplot(2,1,2);plot(t1,Dy);xlabel('时间t');ylabel('Dy(t)'' ');grid绘制结果如下图：练习3做出以参数方程x=cos t,y=sin t,t∈[0,2π]表示的平面曲线（单位圆）。

数学实验报告院系：西安交通大学软件学院软件工程系；班级：软件11；项目：MATLAB软件与基础数学实验；指导教师：张芳；日期：2012年6月11日星期一；学生姓名：贺翔；学号：2111601006；题目【一】在同一坐标系下画出函数y=sin x, y=cos x, y=0.2e0.1x sin (0.5x)和y=0.2e0.1x cos(0.5x)在区间[0，2pi]的曲线图，并对该图进行修饰。

(1)解题思路：首先按步长赋值法生成x向量，则生成相应函数值向量；然后运用plot命令，再添加网格或者其他修饰等。

(2)算法设计：x=0:0.07*pi:2*pi;y1=sin(x);y2=cos(x);y3=0.2.*exp(0.1.*x).*sin(0.5.*x);y4=0.2.*exp(0.1.*x).*cos(0.5.*x);plot(x,y1,'r--',x,y2,'k:',x,y3,'g.',x,y4,'b+','linewidth',3,'markersize',5); grid;xlabel('variable\it{x}')ylabel('variable\it{y}')title('four cruves')text(2.6,0.7,'sin(x)')text(3.5,0.3,'0.2.*exp(0.1.*x).*sin(0.5.*x)')text(5.8,0.8,'cos(x)')text(4.1,-0.4,'0.2.*exp(0.1.*x).*cos(0.5.*x)')(3)结果截图：题目【二】某农夫有一个半径10m的圆形牛栏，长满了草。

他要将一头牛拴在牛栏边界的栏桩上，但只让牛吃到一半草，问栓牛鼻的绳子应为多长？(1)解题思路：设R 为牛栏的半径，而栓牛绳长为r; 则根据数学公式：S=12R 2·4arcsin(r 2R )+ 12r 2·2arccos(r 2R )-2×12r √R 2−r 24；以及令S=12πR 2，即可解出方程的解。

利用Matlab绘制正弦信号的频谱图并做相关分析一、作业要求：1、信号可变（信号的赋值、相位、频率可变）；2、采样频率fs可变；3、加各种不同的窗函数并分析其影响；4、频谱校正；5、频谱细化。

二、采用matlab编写如下程序：clear;clf;fs=100;N=1024; %采样频率和数据点数A=20;B=30;C=0.38;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,1),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');grid on;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,2),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图2：fs=100,N=1024，两种信号叠加');grid on;%加噪声之后的图像x=A*sin(2*pi*B*t+C)+28*randn(size(t));y=fft(x,N);yy=abs(y);yy=yy*2/N; %幅值处理subplot(3,3,3),plot(f(1:N/2.56),yy(1:N/2.56));xlabel('频率/\itHz');ylabel('振幅');title('图3：fs=100,N=1024混入噪声');grid on;%改变采样点数N=128N=128;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,4),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图4：fs=100，N=128');grid on;%改变采样频率为200Hz时的频谱fs=400;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,5),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图5：fs=400，N=1024');grid on;%加三角窗函数fs=100;N=1024; %采样频率和数据点数n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号window=triang(N);%生成三角窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,6),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图6：fs=100,N=1024,加三角窗函数');grid on;%加海明窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hamming(N);%生成海明窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,7),plot(f(1:N/2.56),1.852*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图7：fs=100，N=1024,加海明窗函数');grid on;%加汉宁窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hanning(N);%生成汉宁窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,8),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图8：fs=100，N=1024,加汉宁窗函数');grid on;三、运行结果如下：四、分析与结论：1）从所做图像可以看出，信号的幅值均小于真实值，说明在截断信号时存在泄露。

利用Matlab绘制正弦信号的频谱图并做相关分析一、作业要求：1、信号可变（信号的赋值、相位、频率可变）；2、采样频率fs可变；3、加各种不同的窗函数并分析其影响；4、频谱校正；5、频谱细化。

二、采用matlab编写如下程序：clear;clf;fs=100;N=1024; %采样频率和数据点数A=20;B=30;C=0.38;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,1),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');grid on;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,2),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图2：fs=100,N=1024，两种信号叠加');grid on;%加噪声之后的图像x=A*sin(2*pi*B*t+C)+28*randn(size(t));y=fft(x,N);yy=abs(y);yy=yy*2/N; %幅值处理subplot(3,3,3),plot(f(1:N/2.56),yy(1:N/2.56));xlabel('频率/\itHz');ylabel('振幅');title('图3：fs=100,N=1024混入噪声');grid on;%改变采样点数N=128N=128;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,4),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图4：fs=100，N=128');grid on;%改变采样频率为200Hz时的频谱fs=400;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,5),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图5：fs=400，N=1024');grid on;%加三角窗函数fs=100;N=1024; %采样频率和数据点数n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号window=triang(N);%生成三角窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,6),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图6：fs=100,N=1024,加三角窗函数');grid on;%加海明窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hamming(N);%生成海明窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,7),plot(f(1:N/2.56),1.852*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图7：fs=100，N=1024,加海明窗函数');grid on;%加汉宁窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hanning(N);%生成汉宁窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,8),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图8：fs=100，N=1024,加汉宁窗函数');grid on;三、运行结果如下：四、分析与结论：1）从所做图像可以看出，信号的幅值均小于真实值，说明在截断信号时存在泄露。

matlab正弦函数Matlab正弦函数是工程师和学者在数学、物理、工程领域中使用Matlab计算机软件编程应用最重要的基础函数之一，它主要用于解决圆锥变换，平面定向，坐标转换，代数概念解析，特征分析，轨迹计算，信号处理，模型拟合等方面的问题。

Matlab正弦函数是定义在正弦曲线上的一个函数，其表达式为y = sin(x)，其中x表示弧度，即用来表示单位圆上某点到圆心的线段与x轴正半轴的夹角大小。

y值表示弧度值x相对于圆心的线段，其长度为1。

由于sin(x)函数涉及三维曲线常见的正弦曲线，因此它通常用于描述曲线的各个参数，这些参数主要包括周期、振幅、相移、平均值和中间值。

正弦函数的应用非常广泛，比如，在介绍三角函数时，我们可以用正弦函数来描述三角函数的参数，这样便可以很方便地计算出三角函数的值。

另外，正弦函数也可以用来描述各种微分方程的解，如伯努利微分方程、波动方程等。

正弦函数也常用于各种信号处理的软件程序，它可以提供精确的参数分析，以便更准确地定位出特定信号的部分，这些部分就可以作为分析或处理信号的依据。

正弦函数在数值分析、有限元与有限差分分析等物理仿真中也有广泛的应用，它可以处理复杂的物理问题，从而使我们能够更容易地理解计算机程序的工作原理。

此外，Matlab正弦函数还可以用于统计分析与模型拟合，它可以建立合理的模型，分析数据，并能够以良好的准确度应用于实际中的问题。

例如，它可以用来分析人口统计数据，估测投资回报，建立温度预测模型等。

总之，Matlab正弦函数是一个非常有用的函数，它能够处理复杂的数学模型，从而方便我们进行计算机编程应用程序。

它被广泛应用于数学、物理和工程等领域，用来处理各种复杂的数学问题，使得计算机程序的编程工作更加顺利、高效。

MATLAB实验报告（绘图）MATLAB实验报告——MATLAB绘图学号：学院：班级：姓名：1、绘制y = e x 3sin 3x (x ∈ 0,4π ) 的图像，要求用蓝色的星号画图；并且画出其包络线y = ±e x 3 的图像，用红色的点划线画图. MATLAB Code: x = 0:pi/50:4*pi;y = exp(x/3).*sin(3*x);z(1,:) = exp(x/3);z(2,:) = -exp(x/3);plot(x,y,'b*',x,z(1,:),'r-.',x,z(2,:),'r-.');2、用fplot 和ezplot 命令绘出函数y = e2t 3sin ?(1+2t ) 在区间[1,10]上的图像.MATLAB Code:subplot(1,2,1),fplot(@(t)exp(-2*t/3).*sin(1+2*t),[1,10]); title('fplot');xlabel('t');subplot(1,2,2),ezplot('exp(-2*t/3).*sin(1+2*t)',[1,10]); title('ezplot');3、在同一图形窗口画三个字图，要求使用指令gtext,axis,legend,title,xlabel,ylabel:(1)y=xcosx,x∈(?π,π)(2) y=xtan1xsin x3,x∈(π,4π)(3) y=e 1sin x,x∈[1,8]MATLAB Code:subplot(2,2,1),x1 = -pi:pi/50:pi;y1 = x1.*cos(x1);plot(x1,y1,'r'),axisequal,legend('x*cosx'),title('f1'),xlabel('x'),ylabel('y'); subplot(2,2,2),x2 = pi:pi/50:4*pi;y2 = x2.*tan(1./x2).*sin(x2.^3);plot(x2,y2,'b'),axisequal,legend('x*tan(1/x)*sin(x^3)'),title('f2'),xlabel('x') ,ylabel ('y');subplot(2,2,3:4),x3 = 1:0.01:8;y3 = exp(1./x3).*sin(x3);plot(x3,y3,'g'),legend('exp(1/x)*sinx'),title('f3'),xlabel( '时间'),ylabel('位移');gtext('x*cosx');4、使用合适的单轴对数坐标函数绘制函数y=e x2的图像（其中1≤x≤10）MATLAB Code:x = 1:0.01:10;y = exp(x.^2);semilogy(x,y,'b-.');xlabel('x'),ylabel('y');5、绘制圆锥螺线的图像并添加各种标注，圆锥螺线的参数方程为：x =tcos πt y =tsin π6t z =2t(0≤t ≤20π) MATLAB Code:t = 0:pi/50:20*pi;x = t.*cos(pi/6*t);y = t.*sin(pi/6*t);z = 2*t;plot3(x,y,z,'b'),grid on;title('圆锥螺线');xlabel('x = tcost');ylabel('y = tsint');zlabel('z = 2t');6、在同一个图形窗口画半径为1的球面、柱面x2+y2=1以及极sin4t,t∈[0,2π].坐标图形ρ=12MATLAB Code:subplot(1,2,1),sphere(100);hold on;cylinder;hold on;t = 0:pi/50:2*pi;r = 0.5*sin(4*t);subplot(1,2,2);polar(t,r);title('r = 0.5*sin4t');7、用mesh与surf命令绘制三维曲面z=x2+3y2的图像，并使用不同的着色效果及光照效果.MATLAB Code:t = -3:0.1:3;[x,y] = meshgrid(t);z = x.^2+3*y.^2;subplot(1,2,1),mesh(x,y,z),title('网格z = x^2+3y^2'),shading flat;light('position',[10,10,2]);subplot(1,2,2),surf(x,y,z),title('表面z = x^2+3y^2');shading interp; light('position',[5,-4,-2]);8、绘制由函数x 29+y216+z24=1形成的立体图，并通过改变观测点获得该图形在各个坐标平面的投影.MATLAB Code:t = 0:pi/50:2*pi;[x,y] = meshgrid(t,t);X = 3*sin(y).*cos(x);Y = 4*sin(y).*sin(x);Z = 2*cos(y); subplot(2,2,1);mesh(X,Y,Z);title('x^2/9+y^2/16+z^2/4=1');subplot(2,2,2);mesh(X,Y,Z);view(90,0);title('x^2/9+y^2/16+z^2/4=1在yoz面投影'); subplot(2,2,3);mesh(X,Y,Z);view(0,0);title('x^2/9+y^2/16+z^2/4=1在xoz面投影'); subplot(2,2,4);mesh(X,Y,Z);view(0,90);title('x^2/9+y^2/16+z^2/4=1在xoy面投影');9、画三维曲面z=5?x2?y2?2≤x,y≤2与平面z=3的交线. MATLAB Code:t = -2:0.01:2;[x,y] = meshgrid(t);z1=(5-x.^2-y.^2);z2 = 3*ones(length(t));r0 = (abs(z1-z2)<=0.03);zz = r0.*z2;yy = r0.*y;xx = r0.*x;plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'b.');xlabel('x'),ylabel('y'),zlabel('z');title('z = 5-x^2-y^2 与 z = 3的交线');10(附加)、利用迭代复函数z n+1=z n2+c,c∈C构造Mandelbrot集. Mandelbrot集定义为：M: =c∈C| |c|≤2,|c2+c|≤2,|(c2+c)2+c|≤2…….MATLAB Code:%最大迭代数网格细度%本次作图取maxloopNum = 100,Eps = 0.05function drawMandelbrot(maxloopNum,Eps)t = -2:Eps:2;[x,y] = meshgrid(t);z = x+i*y;isMandel = zeros(length(t));for j = 1:length(z)for k = 1:length(z)iter = 0;current = z(j,k);temp = current;while abs(temp)<=2 && iter<maxloopnum< p="">temp = temp^2+current;iter = iter + 1;endif iter == maxloopNumisMandel(j,k) = 1;endendendplot(z.*isMandel,'.');axis equal;title('MandelBrot set');legend('loop = 100,eps = 0.05','Location','northwest'); xlabel('x'),ylabel('iy');end</maxloopnum<>。

利用Matlab绘制正弦信号的频谱图并做相关分析一、作业要求：1、信号可变（信号的赋值、相位、频率可变）；2、采样频率fs可变；3、加各种不同的窗函数并分析其影响；4、频谱校正；5、频谱细化。

二、采用matlab编写如下程序：clear;clf;fs=100;N=1024; %采样频率和数据点数A=20;B=30;C=0.38;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,1),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');grid on;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,2),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图2：fs=100,N=1024，两种信号叠加');grid on;%加噪声之后的图像x=A*sin(2*pi*B*t+C)+28*randn(size(t));y=fft(x,N);yy=abs(y);yy=yy*2/N; %幅值处理subplot(3,3,3),plot(f(1:N/2.56),yy(1:N/2.56));xlabel('频率/\itHz');ylabel('振幅');title('图3：fs=100,N=1024混入噪声');grid on;%改变采样点数N=128N=128;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,4),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图4：fs=100，N=128');grid on;%改变采样频率为200Hz时的频谱fs=400;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,5),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图5：fs=400，N=1024');grid on;%加三角窗函数fs=100;N=1024; %采样频率和数据点数n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号window=triang(N);%生成三角窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,6),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图6：fs=100,N=1024,加三角窗函数');grid on;%加海明窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hamming(N);%生成海明窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,7),plot(f(1:N/2.56),1.852*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图7：fs=100，N=1024,加海明窗函数');grid on;%加汉宁窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hanning(N);%生成汉宁窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,8),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图8：fs=100，N=1024,加汉宁窗函数');grid on;三、运行结果如下：四、分析与结论：1）从所做图像可以看出，信号的幅值均小于真实值，说明在截断信号时存在泄露。

matlab实验报告实验报告：Matlab实验分析1. 实验目的本实验旨在通过Matlab软件完成一系列数值计算和数据分析的任务，包括绘制曲线、解方程、矩阵运算等，以加深对Matlab软件的理解和掌握。

2. 实验内容2.1 绘制函数曲线首先，我们通过在Matlab中输入函数的表达式来绘制函数曲线。

例如，我们可以输入y = sin(x)来绘制正弦函数的曲线。

另外，我们还可以设置曲线的颜色、线型和坐标轴范围等。

2.2 解方程接下来，我们使用Matlab来解方程。

对于一元方程，我们可以使用solve函数来求出方程的解。

例如，我们输入syms x; solve(x^2 - 2*x - 8)来解方程x^2 - 2x - 8 = 0。

而对于多元方程组，我们可以使用solve函数的向量输入形式来求解。

例如，我们输入syms x y; solve(x^2 + y^2 - 1, x - y - 1)来求解方程组x^2 + y^2 - 1 = 0和x - y - 1 = 0的解。

2.3 矩阵运算Matlab也可以进行矩阵运算。

我们可以使用矩阵相乘、相加和取逆等运算。

例如，我们可以输入A = [1 2; 3 4]和B = [5 6;7 8]来定义两个矩阵，然后使用A * B来计算它们的乘积。

3. 实验结果与分析在本实验中，我们成功完成了绘制函数曲线、解方程和矩阵运算等任务。

通过Matlab软件，我们可以快速、准确地进行数值计算和数据分析。

使用Matlab的高级函数和工具箱，我们可以更方便地处理复杂的数值计算和数据分析问题。

4. 实验总结通过本次实验，我们进一步加深了对Matlab软件的理解和掌握。

Matlab提供了丰富的函数库和工具箱，适用于各种不同的数值计算和数据分析任务。

在日常科研和工程实践中，Matlab是一个非常强大和方便的工具，可以帮助我们更高效地完成任务。

⽤matlab画漂亮的sin曲线Technorati 标记: ,相信⼤部分⽤过matlab的⼈都画过sin曲线，直接plot就可以了，不过呢，plot出来的曲线⾃然不那么好看，本着绳命在于折腾的原则，⼩弟学习了下sin曲线的⾼级画法，mark⼀下，还是先上图看看是如何华丽丽的完整转变的吧。

实现代码如下：x = -pi:0.1:pi;y = sin(x);axescenterfid1 = plot(x,y);set(findobj('Type','Line'),'Color','r','LineWidth',2)set(gca,'ytick',[-1,-0.5,0,0.5,1],'LineWidth',0.1)xlabel('-\pi \leq \Theta \leq \pi')ylabel('f(\Theta)','fontname','')set(gca,'ylim',[-2,2],'xtick',[-pi,-pi/2,0,pi/2,pi],'xticklabel',{'-p','-p/2','0','p/2','p'},'fontname','symbol')text(-pi/2,sin(-pi/2),'\uparrow f(\Theta) = sin(\Theta)','VerticalAlignment','top','Fontsize',14)hold on;pp = plot(x,cos(x));set(pp,'Color','b','LineWidth',2)text(-0.85,cos(-0.85),'f(\Theta) = cos(\Theta)\rightarrow','VerticalAlignment','bottom','Fontsize',14,'horizontalAlignment','right')其中为了把坐标轴移动到原点（发现这个才是⼤问题），使⽤了⽹上⼀份代码（对⽐了⽹上⼏个实现⽅案，个⼈认为这个代码最佳）其他参考matlab帮助： axes, set, gca, xdir, findobj。

Matlab实验报告（三）MATLAB绘图实验⽬的1．掌握MATLAB的基本绘图命令。

2．掌握运⽤MATLAB绘制⼀维、⼆维、三维图形的⽅法。

3．给图形加以修饰。

⼀、预备知识1．基本绘图命令plotplot绘图命令⼀共有三种形式：⑴plot(y)是plot命令中最为简单的形式，当y为向量时，以y的元素为纵坐标，元素相应的序列号为横坐标，绘制出连线；若y为实矩阵，则按照列绘出每列元素和其序列号的对应关系，曲线数等于矩阵的列数；当y为复矩阵时，则按列以每列元素的实部为横坐标，以虚部为纵坐标，绘出曲线，曲线数等于列数。

⑵ plot(x,y,[linspec])其中linspec是可选的，⽤它来说明线型。

当x和y为同维向量时，以x为横坐标，y为纵坐标绘制曲线；当x是向量，y是每⾏元素数⽬和x维数相同的矩阵时，将绘出以x为横坐标，以y中每⾏元素为纵坐标的多条曲线，曲线数等于矩阵⾏数；当x为矩阵，y为相应向量时，使⽤该命令也能绘出相应图形。

⑶plot(x1,y1,x2,y2,x3,y3……)能够绘制多条曲线，每条曲线分别以x和y为横纵坐标，各条曲线互不影响。

线型和颜⾊MATLAB可以对线型和颜⾊进⾏设定，线型和颜⾊种类如下：线：—实线：点线—.虚点线——折线点：.圆点 +加号 *星号 x x型 o 空⼼⼩圆颜⾊：y 黄 r 红 g 绿 b 蓝 w ⽩ k ⿊ m 紫 c 青特殊的⼆维图形函数表5 特殊2维绘图函数[1] 直⽅图在实际中，常会遇到离散数据，当需要⽐较数据、分析数据在总量中的⽐例时，直⽅图就是⼀种理想的选择，但要注意该⽅法适⽤于数据较少的情况。

直⽅图的绘图函数有以下两种基本形式。

·bar(x,y) 绘制m*n 矩阵的直⽅图。

其中y 为m*n 矩阵或向量，x 必须单向递增。

·bar(y) 绘制y 向量的直⽅图，x 向量默认为x=1:m close all; %关闭所有的图形视窗。

MATLAB实验报告1MATLAB实验报告1实验题目：使用MATLAB绘制正弦函数图像一、实验目的通过本次实验，我们将学会使用MATLAB软件绘制函数图像，加深对正弦函数性质的理解，并掌握MATLAB基本绘图函数的使用。

二、实验原理正弦函数是数学中常见的一种函数，其表达式为y = A * sin(ω *t + φ)，其中A为振幅，ω为角速度，t为时间，φ为相位角。

三、实验步骤1.打开MATLAB软件，在命令窗口输入以下指令创建一个时间向量t，并以0.01为间隔，从0到2π生成值：t = 0:0.01:2*pi;2.输入以下指令生成正弦函数序列y，并设定振幅A为1，角速度ω为1，相位角φ为0：y = sin(t);3. 使用plot函数绘制正弦函数图像，输入以下指令：plot(t, y);title('Sin(x) Function');xlabel('Time');ylabel('Amplitude');xlim([0, 2*pi]);5.运行程序，即可获得正弦函数的图像。

四、实验结果分析（插入正弦函数图像）从图中可以看出，正弦函数以x轴为中心呈现周期性的波动，振幅为1、当x取0时，正弦函数的值为0；当x取π/2时，正弦函数的值为1；当x取π时，正弦函数的值为0；当x取3π/2时，正弦函数的值为-1；当x取2π时，正弦函数的值为0，周期为2π。

五、实验心得体会通过本次实验，我对MATLAB的基本绘图函数的使用有了初步了解和实践，提高了自己的实际操作能力，并对正弦函数的性质有了深入的认识。

在今后的实验中，我将继续深入学习和应用MATLAB软件，提高自己的数据分析和图像处理能力。

matlab符号运算正弦函数
MATLAB是一款功能强大的数学计算软件，它提供了符号运算功能，可以用于计算各种数学表达式，包括正弦函数。

在MATLAB中，可以使用符号运算来计算正弦函数。

下面是一个简单的示例：matlab
syms x
y = sin(x);
在上面的代码中，syms命令用于声明符号变量x，然后使用sin函数计算正弦函数。

如果你想要绘制正弦函数的图像，可以使用plot函数。

下面是一个示例：matlab
x = linspace(-pi, pi, 100);
y = sin(x);
plot(x, y);
在上面的代码中，linspace函数用于生成一个包含100个元素的向量x，范围从-π到π。

然后使用sin函数计算每个元素的正弦值，并将结果存储在向量y中。

最后使用plot函数绘制出正弦函数的图像。

如果你想要计算正弦函数的特定值，可以使用数值计算方法。

下面是一个示例：
matlab
x = 0.5;
y = sin(x);
disp(y);
在上面的代码中，x是一个数值变量，其值为0.5。

然后使用sin函数计算其正弦值，并使用disp函数显示结果。

总的来说，MATLAB的符号运算功能可以用于计算各种数学表达式，包括正弦函数。

你可以使用不同的方法来计算正弦函数的值或绘制其图像。

matlab三角函数写法Matlab是一款广泛应用于数学和工程领域的数学软件，它提供了丰富的三角函数库，使得我们可以轻松地编写和求解三角函数问题。

本篇文章将详细介绍Matlab中常用的三角函数，以及如何正确使用它们进行绘图和计算。

一、正弦函数与余弦函数在Matlab中，正弦函数和余弦函数的符号分别为sin和cos。

它们的输入参数为角度（弧度制），输出为该角度的正弦值或余弦值。

例如，以下代码将绘制一个以0到2π为范围的圆，并标出每个角度对应的正弦值和余弦值：```matlabtheta = 0:pi/50:2*pi; % 生成角度向量sin_val = sin(theta); % 计算正弦值cos_val = cos(theta); % 计算余弦值plot(theta, sin_val); % 绘制正弦函数图像hold on; % 保持当前图像，以便绘制其他图形plot(theta, cos_val); % 绘制余弦函数图像legend('sin(x)', 'cos(x)'); % 添加图例说明```二、正切函数与余切函数Matlab中，正切函数和余切函数的符号分别为tan和cot。

它们的输入参数也为角度（弧度制），输出为该角度的正切值或余切值。

需要注意的是，在Matlab中，角度的范围通常默认为度数制，但为了与角度为弧度制的计算相匹配，我们在编写代码时需要注意单位转换。

三、其他三角函数除了正弦、余弦和正切函数外，Matlab还提供了其他一些常用的三角函数，如正弦的平方函数sin^2、反正切函数atan等。

这些函数的符号分别为sin^2、acos等。

在使用这些函数时，需要注意输入参数的单位和范围，以确保计算结果的准确性。

四、绘图与计算在Matlab中，我们可以使用plot函数绘制三角函数的图像，使用相关函数进行数值计算和统计分析。

例如，以下代码将绘制一个以0到π为范围的三角函数的图像，并使用相关函数进行数值计算：```matlabtheta = linspace(0, pi, 100); % 生成角度向量sin_val = sin(theta); % 计算正弦值cos_val = cos(theta); % 计算余弦值atan_val = atan(theta); % 计算反正切值plot(theta, sin_val, 'r'); % 绘制正弦函数图像，颜色为红色hold on; % 保持当前图像，以便绘制其他图形plot(theta, cos_val, 'b'); % 绘制余弦函数图像，颜色为蓝色plot(theta, atan_val, 'g'); % 绘制反正切函数图像，颜色为绿色legend('sin(x)', 'cos(x)', 'atan(x)'); % 添加图例说明```以上就是Matlab中三角函数的常用写法及绘图方法。

利用Matlab绘制正弦信号的频谱图并做相关分析一、作业要求：1、信号可变（信号的赋值、相位、频率可变）；2、采样频率fs可变；3、加各种不同的窗函数并分析其影响；4、频谱校正；5、频谱细化。

二、采用matlab编写如下程序：clear;clf;fs=100;N=1024; %采样频率和数据点数A=20;B=30;C=0.38;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,1),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');grid on;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,2),plot(f,yy); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图2：fs=100,N=1024，两种信号叠加');grid on;%加噪声之后的图像x=A*sin(2*pi*B*t+C)+28*randn(size(t));y=fft(x,N);yy=abs(y);yy=yy*2/N; %幅值处理subplot(3,3,3),plot(f(1:N/2.56),yy(1:N/2.56));xlabel('频率/\itHz');ylabel('振幅');title('图3：fs=100,N=1024混入噪声');grid on;%改变采样点数N=128N=128;n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,4),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图4：fs=100，N=128');grid on;%改变采样频率为200Hz时的频谱fs=400;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,5),plot(f(1:N/2.56),yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图5：fs=400，N=1024');grid on;%加三角窗函数fs=100;N=1024; %采样频率和数据点数n=0:N-1;t=n/fs; %时间序列x=A*sin(2*pi*B*t+C); %信号window=triang(N);%生成三角窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy*2/N; %幅值处理f=n*fs/N; %频率序列subplot(3,3,6),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图6：fs=100,N=1024,加三角窗函数');grid on;%加海明窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hamming(N);%生成海明窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,7),plot(f(1:N/2.56),1.852*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图7：fs=100，N=1024,加海明窗函数');grid on;%加汉宁窗函数后的频谱fs=100;N=1024;n=0:N-1;t=n/fs;x=A*sin(2*pi*B*t+C); %信号window=hanning(N);%生成汉宁窗函数x=x.*window';%加窗函数y=fft(x,N); %对信号进行快速傅里叶变换yy=abs(y); %求取傅里叶变换的振幅yy=yy*2/N; %幅值处理f=n*fs/N;subplot(3,3,8),plot(f(1:N/2.56),2*yy(1:N/2.56)); %绘出随频率变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图8：fs=100，N=1024,加汉宁窗函数');grid on;三、运行结果如下：四、分析与结论：1）从所做图像可以看出，信号的幅值均小于真实值，说明在截断信号时存在泄露。

南昌航空大学数学与信息科学学院实验报告课程名称：数学实验实验名称： MATLAB基本操作实验类型：验证性■综合性□ 设计性□实验室名称：数学实验室班级学号： 10 学生姓名：钟 X任课教师（教师签名）：成绩：实验日期： 2011-１０- 1０南昌航空大学数学与信息科学学院实验报告一、实验目的1、熟悉MATLAB基本命令与操作2、熟悉MATLAB作图的基本原理与步骤3、学会用matlab软件做图二、实验用仪器设备、器材或软件环境计算机MATLAB软件三、实验原理、方案设计、程序框图、预编程序等问题1：在区间【0,2π】画sinx实验程序:>> x=linspace(0,2*pi,30);>> y=sin(x);>> plot(x,y)问题2：在【0,2π】用红线画sinx,用绿圈画cosx,实验程序:第 1 页>> x=linspace(0,2*pi,30);>> y=sin(x);>> z=cos(x);>> plot(x,y,'r',x,z,'co')>>问题3：在【0，π】上画y=sinx的图形。

实验程序:>> ezplot('sin(x)',[0,pi])>>问题4：在【0，π】上画x=cos3t,y=sin3t星形图形。

南昌航空大学数学与信息科学学院实验报告实验程序:>> ezplot('cos(t).^3','sin(t).^3',[0,pi])>>问题5：[-2,0.5]，[0,2]上画隐函数实验程序:>> ezplot('exp(x)+sin(x*y)',[-2,0.5,0,2])>>问题6：在[－２，２]范围内绘制tanh的图形。

实验报告2课程 数学实验与数学软件 实验名称 matlab 的二维绘图 第 1 页专业 数学与应用数学 班级_08级数学3班__ ______ 学号_182__ 姓名 郑贞贞 实验日期：2010 年 03 月18 日 评分一、实验目的1. 了解plot 函数的用法和matlab 绘图的基本原理。

2. 了解图形属性的设置。

二、实验内容1．用 plot 函数在一个坐标系下绘制以下函数的图形：y 1 用黑色间断线点标记为星号 y 2 用红色实线点标记为小方格 y 3 用蓝色虚线点标记为小圆圈2.在不同的图形子窗口中绘制如下曲线：Y=sin(x)的线型为虚线，点为菱形，点标记表面色为红； Y=cos(x)的线型为－.，点为左三角形，点标记表面色为黑； Y=1.2^x 的线型为实线，点为o 形，点标记表面色为黄; 字体大小为15； 线粗细为2；123sin(),cos(),sin(2), 02y x y x y x x π===≤≤3．使用text, title, xlabel和plot绘制下面图形：4、如何使用subplot命令创建下面的图形窗口：5、在同一个窗口中画出三条曲线，并给出图例和栅格线三、实验环境（1）硬件：PC机（2）软件：Windows 操作系统、matlab7.0四、实验步骤1、x=0:pi/10:2*pi;y1=sin(x);y2=cos(x);y3=sin(2*x);plot(x,y1,'k-.*',x,y2,'r-s',x,y3,'b:o')2、x=-2*pi:pi/10:2*pi;subplot(2,2,1)y=sin(x);plot(x,y,'d:','linewidth',2, 'markerfacecolor','r') title('y=sin(x)函数图像','fontsize',15);xlabel('x\in[-2*pi,2*pi]')legend('sin(x)')subplot(2,2,3)x=0:pi/10:2*pi;y=cos(x);plot(x,y,'<-','linewidth',2, 'markerfacecolor','k') title('y=cos(x)函数图像','fontsize',15)xlabel('x\in[0,2*pi]')legend('cos(x)')subplot(2,2,[2 4])x=-10:1:100;y=1.2.^x;plot(x,y,'o-','linewidth',2, 'markerfacecolor','y') title('y=1.2^x函数图像','fontsize',15)xlabel('x\in[-10,100]')ylabel('y=1.2^x','fontsize',15)3、x=0:pi/10:2*pi;y=sin(2*x-1);plot(x,y,'ro-','markerfacecolor','b')xlabel('x\in[0,2*pi]')title('{\omega}=2,{\beta}=-1')text(5,0.4,'\leftarrowy=sin({\omega}x+{\beta})') 4、x=-4:0.5:4;y=x.*x;figuresubplot(2,2,1);plot(x,y)y=2.^x;subplot(2,2,2);plot(x,y)y=-x.*x;subplot(2,2,[3,4]);subplot('position',[0.29,0.1,0.5,0.38]);plot(x,y)5x=0:3:60;y=2*x-1;figure;plot(x,y,'color','b')hold ony=-2*x;plot(x,y,'color','r');hold ony=x;plot(x,y,'color','g');legend('y=2*x-1','y=-2*x','y=x') grid on五、实验结果与讨论实验结果如下：x∈[-2*pi,2*pi] x∈[0,2*pi]7y=1.2x函数图像x∈[-10,100]y=1.2xx ∈[0,2*pi]-4-2024-4-20240102030405060-150-100100150六、总结（说明实验过程中遇到的问题及解决办法；新发现或个人的收获；未解决/需进一步研讨的问题或建议新实验方法等）问题：没有注意到点乘，出现了一些差错，后来通过翻书了解了，偶尔用中文输入法忘了改标点符号收获：学会在同一个窗口显示几个图像，及安排各个子窗口位置不足：需要对各种函数多加了解，对各种线型点型符号需要了解。

matlab sin函数从古至今，数学家们一直在努力认识数学中的三角形，Sin函数就是精确描绘出数学三角形的有效方法之一。

MATLAB中的Sin函数就是以数学技巧将正弦函数的图像反映在数学三角形上，使我们能够更为精确的对其进行研究和分析。

一、MATLAB中的Sin函数MATLAB中的Sin函数是一种精确描绘数学三角形的功能，它基于正弦函数来实现，用来计算平面上一点的坐标与圆环之间的关系。

其中，正弦函数部分反映数学中三角形的轴线关系，可以通过它来推算三角形的其他特征，如角度、面积等。

因此，MATLAB中的Sin函数就是一个可以精确描绘一个数学三角形的图像的工具。

它可以按照输入的数据，计算出三角形的角度、面积、边长等特性，并以图形的方式进行显示。

二、MATLAB中Sin函数的使用MATLAB中的Sin函数可以用于不同的科学研究，只要输入有关数学三角形的数据，就可以精确的描绘出其图像。

其使用方法如下：（1）首先，定义函数，输入数学三角形的数据，例如顶点坐标，内角大小，边长等；（2）接着，选择函数解析方法，如MATLAB的syms功能，用来求解每个变量的值；（3）最后，根据函数参数，将数学三角形以图形的方式显示出来。

三、MATLAB中Sin函数的应用MATLAB中的Sin函数可以应用于不同领域，用于不同科学研究。

（1） MATLAB中的Sin函数可以应用于测绘学中，用来描绘地理形状，包括山脊等山区、海洋与陆地的边界、不规则的川流曲线等；（2）MATLAB中的Sin函数可以应用于天文学领域，用来表示行星运行的轨道，以及星空图案；（3）MATLAB中的Sin函数可以应用于几何学领域，用来求解平面或空间中图形的轮廓线，如圆、正方形、三角形等。

四、总结MATLAB中的Sin函数既可用于学术研究，也可用于工程应用，只要输入正确的数学参数，就可以精确描绘出数学三角形的图像，使我们能够更为精确的对其进行研究和分析。

MATLAB中的Sin函数的使用方法简单易行，它的应用也非常广泛，可以应用于各个领域，为解决各式各样的科学问题提供便利。

实验7：图形绘制一、实验目的1、 掌握绘制二维图形的常用函数。

2、 掌握绘制三维图形的常用函数。

3、 掌握绘制图形的辅助操作。

二、实验内容1、 已知2*13),2cos(2,12y y y x y x y ===，完成下列操作：（1） 在同一坐标系下用不同的颜色和线型绘制三条曲线。

%homework_7_1_1.mx=0:pi/100:2*pi;y1=x.*x;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'r--',x,y2,'k:',x,y3,'b-.'); （2） 以子图形式绘制三条曲线。

%homework_7_1_2.mx=0:pi/100:2*pi;y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2, 2, 1);plot(x,y1,'r--');subplot(2, 2, 2);plot(x,y2,'k:');subplot(2, 2, 3);plot(x,y3,'b-.');（3） 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。

%homework_7_1_3.m%ÌõÐÎͼ¡¢½×ÌÝͼ¡¢¸ËͼºÍÌî³äͼx=0:pi/100:2*pi;y1=x.*x;y2=cos(2*x);y3=y1.*y2;%µÚÒ»ÐÐsubplot(4, 3, 1);bar(x,y1,'r');subplot(4, 3, 2);bar(x,y2,'k');subplot(4, 3, 3);bar(x,y3,'b'); %µÚ¶þÐÐstairs(x,y1,'r');subplot(4, 3, 5);stairs(x,y2,'k');subplot(4, 3, 6);stairs(x,y3,'b');%µÚÈýÐÐsubplot(4, 3, 7);stem(x,y1,'r');subplot(4, 3, 8);stem(x,y2,'k');subplot(4, 3, 9);stem(x,y3,'b');%µÚËÄÐÐsubplot(4, 3, 10);fill(x,y1,'r');subplot(4, 3, 11);fill(x,y2,'k');subplot(4, 3, 12);fill(x,y3,'b');2、 绘制极坐标曲线)sin(θρn b a +=，并分析参数a ，b ，n 对曲线形状的影响。