许莹莹MATLAB实验2

合集下载

Matlab实验指导书(含答案)详解

Matlab实验指导书(含答案)详解

实验一:Matlab操作环境熟悉一、实验目的1.初步了解Matlab操作环境。

2.学习使用图形函数计算器命令funtool及其环境。

二、实验内容熟悉Matlab操作环境,认识命令窗口、内存工作区窗口、历史命令窗口;学会使用format命令调整命令窗口的数据显示格式;学会使用变量和矩阵的输入,并进行简单的计算;学会使用who和whos命令查看内存变量信息;学会使用图形函数计算器funtool,并进行下列计算:1.单函数运算操作。

➢求下列函数的符号导数(1) y=sin(x);(2)y=(1+x)^3*(2-x);➢求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3)y=1/sqrt(1—x^2);(4) y=(x-1)/(x+1)/(x+2);➢求反函数(1) y=(x—1)/(2*x+3);(2) y=exp(x);(3)y=log(x+sqrt(1+x^2));➢代数式的化简(1) (x+1)*(x—1)*(x—2)/(x—3)/(x-4);(2) sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x—3*cos(x)+4*x*sin(x);2.函数与参数的运算操作。

➢从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2)y2=(x+2)^2(3) y3=2*x^2(4)y4=x^2+2(5)y5=x^4(6)y6=x^2/23.两个函数之间的操作➢求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5➢乘积(1) exp(—x)*sin(x)(2)sin(x)*x➢商(1)sin(x)/cos(x);(2)x/(1+x^2);(3)1/(x-1)/(x-2);➢求复合函数(1) y=exp(u) u=sin(x)(2)y=sqrt(u) u=1+exp(x^2)(3)y=sin(u)u=asin(x)(4) y=sinh(u)u=—x三、设计提示1.初次接触Matlab应该注意函数表达式的文本式描述。

(完整word版)MATLAB实训实验讲解

(完整word版)MATLAB实训实验讲解

2015/2016学年下学期《信号与系统》实验报告班级:学号:学生姓名:指导教师:2016年3月8 日实验一 基本函数仿真实验项目: 基本函数仿真实验时间: 2016年 3 月 8 日 星期 二 第 34 节课 实验地点: 1501实验室 实验目的:1、 学习使用MATLAB 软件2、 学习MATLAB 中各种函数,并应用函数分析3、 对MATALB 的进一步的学习了解,熟练掌握MATALB 的各种操纵,学会使用MATALB 解决复杂的运算并学会用MATALB 解决平时学习4、 了解MATALB 的数值运算5、 了解MATALB 的基本函数和命令6、 学习掌握MATALB 有关命令 实验内容: 1、(1) 题目:应用MA TLAB 方法实现单位阶跃信号和矩形脉冲。

(2) 程序清单(源程序)解:对于阶跃函数,MATLAB 中有专门的stairs 绘图命令。

例如,实现)(t 和矩形脉冲的程序如下:t=-1:2; % 定义时间范围向量t x=(t>=0);subplot(1,2,1),stairs(t,x);axis([-1,2,-0.1,1.2]); grid on % 绘制单位阶跃信号波形 t=-1:0.001:1; % 定义时间范围向量t g=(t>=(-1/2))-(t>=(1/2));subplot(1,2,2),stairs(t,g);axis([-1,1,-0.1,1.2]); grid on % 绘制矩形脉冲波形(3) 运行结果(截图)00.20.40.60.8100.20.40.60.81图1 例1图(4)函数解析Subplot:使用方法:subplot (m,n,p )或者subplot (m n p )。

是将多个图画到一个平面上的工具。

其中,m 表示是图排成m 行,n 表示图排成n 列,也就是整个figure 中有n 个图是排成一行的,一共m 行,如果m=2就是表示2行图。

Matlab第二次实验报告

Matlab第二次实验报告

Matlab第二次实验报告实验目的:1、了解plot函数和subplot函数的基本用法和matlab绘图的基本原理。

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

比如画图的颜色,画图采用的线性标识符等。

二:实验基本知识1:1. 单窗口单曲线绘图;2. 单窗口多曲线绘图;3. 单窗口多曲线分图绘图;4. 多窗口绘图;5.可任意设置颜色与线型;6.图形加注功能;7.fplot——绘制函数图函数;8.ezplot——符号函数的简易绘图函数9:subplot—将画图区域分块函数。

实验内容:将高数课本后的18个图画出上机练习程序第1——6图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1);x=-50:1:50;y=a.*x.^3;subplot(3,2,1);plot(x,y,'r');title('y=a*x^3');xlabel('X');ylabel('Y');x=0:1:50;y=sqrt(a.*x.^3);subplot(3,2,2);plot(x,y,x,-y,'r');title('y^2=a*x^3');xlabel('X');ylabel('Y');x=-3:0.1:3;y=exp(-x.^2);subplot(3,2,3);plot(x,y,'r');title('y=e^x^2');xlabel('X');ylabel('Y');x=-3:0.1:3;y=8*a^3./(x.^2+4*a^2);subplot(3,2,4);plot(x,y,'r');title('y=8*a^3./(x.^2+4*a^2)'); xlabel('X');ylabel('Y');x=0:0.01:5;y=sqrt(x.^3./(a*2.-x)); subplot(3,2,5);plot(x,y,x,-y,'r');title('y^2*(2a-x)=x^3'); xlabel('X');ylabel('Y');%t=-1:0.01:5;%x=3*a.*t./(1+t.^3);%y=3*a*t.^2./(1+t.^3);subplot(3,2,6);ezplot('x.^3+y.^3-3*3*x.*y'); %plot(x,y);title('x^3+y^3-3axy=0'); xlabel('X');ylabel('Y');输入a=3:显示结果第7——10图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0:0.1:2*pi;x=a*cos(t).^3;y=a*sin(t).^3;subplot(2,2,1);plot(x,y,'b')title('x^1.5+y^1.5=a^1.5'); xlabel('X');ylabel('Y');t=-2*pi:0.1:2*pi;x=a.*(t-sin(t));y=a.*(1-cos(t));subplot(2,2,2);plot(x,y,'b')title('°ÚÏß');xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*(1-cos(t));subplot(2,2,3);plot(p.*cos(t),p.*sin(t),'b'); title('p=a(1-cos£¨t£©)'); xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*t;subplot(2,2,4);plot(p.*cos(t),p.*sin(t),'b'); title('p=at');xlabel('X');ylabel('Y');输入a=1,第11题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º') figure(1);t=-2*pi:pi/100:2*pi;p=exp(a.*t);plot(p.*cos(t),p.*sin(t)); title('p=e^at');xlabel('X');ylabel('Y');输入a=0.1,显示结果:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0.2*pi:pi/100:100*pi;p=a./t;plot(p.*cos(t),p.*sin(t));title('p=at');xlabel('X');ylabel('Y');输入a=1,显示结果:第13——14题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º');figure(1);t=0:pi/100:2*pi;p=(a^2.*sin(2*t)).^0.5;subplot(1,2,1);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*sin(2t)');xlabel('X');p=(a^2.*cos(2*t)).^0.5;subplot(1,2,2);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第15-16题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1)t=0:pi/200:pi;p=a.*cos(3*t);subplot(1,2,1);plot(p.*cos(t),p.*sin(t),'r--');title('p=a*cos(3t)');xlabel('X');p=a.*sin(3*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'r--'); title('p=a*sin(3t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第17——18题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1);t=-pi:pi/200:pi;subplot(1,2,1);p=a.*sin(2*t);plot(p.*cos(t),p.*sin(t),'k-.'); title('p=a*sin(2t)');xlabel('X');ylabel('Y');p=a.*cos(2*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'k-.');title('p=a*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:实验心得:这次实验主要是学习matlab软件的二维绘图功能:matlab软件还有丰富的图形修饰功能,如改变线条的形式和颜色;除此之外其还可以以多种形式如在一个图形中显示多个函数图形,一个figure中显示多个坐标系,还有一个程序中就可以用多个figure一起来表示不同的函数……在画图过程中你可能因为角度的区间取值不一样,画出来的图形可能跟书本上给出的图形有很大的差异,但没关系,只要你思路是对的,区间可以自己慢慢改动知道跟书本上给出的图形想接近。

MATLAB实验二傅里叶分析与应用

MATLAB实验二傅里叶分析与应用

实验二傅里叶分析及应用一、实验目的(一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义2、学会使用Matlab分析周期信号的频谱特性(二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质1、学会运用Matlab求连续时间信号的傅里叶变换2、学会运用Matlab求连续时间信号的频谱图3、学会运用Matlab分析连续时间信号的傅里叶变换的性质(三)掌握使用Matlab完成信号抽样并验证抽样定理1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化3、学会运用MATLAB对抽样后的信号进行重建二、实验条件Win7系统,MATLAB R2015a三、实验内容1、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。

符号运算法数值运算法012345t(20 π ex p(-3 t) heaviside(t) - 8 π ex p(-5 t) heaviside(t))/(2 π)2、试用Matlab 命令求的傅里叶反变换,并绘出其时域信号图。

两个单边指数脉冲的叠加3、已知门函数自身卷积为三角波信号,试用Matlab 命令验证FT 的时域卷积定理。

Code :syms t ; fw =sym('10/(3+i*w)-4/(5+i*w)');ft = ifourier(fw, t); ezplot(ft), grid on ;Code :f = sym('heaviside(t+1) - heaviside(t-1)'); fw = simplify(fourier(f)); F = fw.*fw; subplot(211);ezplot(abs(F), [-9, 9]), grid on title('FW^2') tri =sym('(t+2)*heaviside(t+2)-2*t*heaviside(t)+(t-2)*heaviside(t-2)'); Ftri = fourier(tri); F = simplify(Ftri); subplot(212);ezplot(abs(F), [-9, 9]), grid on ; title('tri FT')4、设有两个不同频率的余弦信号,频率分别为,;现在使用抽样频率对这三个信号进行抽样,使用MATLAB 命令画出各抽样信号的波形和频谱,并分析其频率混叠现象> > > > > > > > > >Code :f1 = 100; % f1 = 100 hz ts = 1/4000;% sample = 4000hz dt = 0.0001;t1 = -0.007:dt:0.007; ft = cos(2*f1*pi*t1);subplot(221); plot(t1, ft), grid on ;axis([-0.006 0.006 -1.5 1.5]) xlabel('Time/s'),ylabel('f(t)') title('Cosine curve');N = 5000; k = -N:N; w = 2*pi*k/((2*N+1)*dt); fw = ft*dt*exp(-1i*t1'*w); subplot(222);plot(w, abs(fw)); grid on ; axis([-20000 20000 0 0.005]); t2 = -0.007:ts:0.007;fst = cos(2*f1*pi*t2);subplot(223);plot(t1, ft, ':'),hold onstem(t2, fst), grid on ;axis([-0.006 0.006 -1.5 1.5])xlabel('Time/s'),ylabel('fs(t)')title('Sample signal'); hold offfsw=ts*fst*exp(-1i*t2'*w);subplot(224); plot(w, abs(fsw)),grid onaxis([-20000 20000 0 0.006])xlabel('\omega'),ylabel('fsw')title(' Sample freq spectrum');-505x 10-3-11Time/s f (t )Cosine curve-2-1012x 104012345-3ωf (w )Cos freq spectrum-505x 10-3-11Time/sf s (t )Sample signal-2-1012x 1040246-3ωf s wSample freq spectrumx 10-3-101Time/s f (t )Cosine curvex 104012345-3ωf (w )Cos freq spectrum-505x 10-3-11Time/sf s (t )Sample signal-2-1012x 1040246-3ωf s wSample freq spectrumf1 = 100Hz将代码中f1设为3800即可↓f2 = 3800Hz-0.500.51Sa(t)0.511.5Sa(t) freq spectrum-0.500.51Sampling signal0.511.5spectrum of Sampling signal5、结合抽样定理,利用MATLAB 编程实现信号经过冲激脉冲抽样后得到的抽样信号及其频谱[建议:冲激脉冲的周期分别取4*pi/3 s 、pi s 、2*pi/3 s 三种情况对比],并利用构建信号。

MATLAB实验报告2

MATLAB实验报告2

实验二用MATLAB实现线性系统的时域分析[实验目的]1.研究线性系统在典型输入信号作用下的暂态响应;2.熟悉线性系统的暂态性能指标;3.研究二阶系统重要参数阻尼比ξ对系统动态性能的影响;4.熟悉在MATLAB下判断系统稳定性的方法;5.熟悉在MATLAB下求取稳态误差的方法。

[实验原理]MATLAB中有两类用于求解系统时域响应的方法。

一、用MATLAB函数(命令)进行暂态响应分析1 求取线性连续系统的单位阶跃响应的函数——step基本格式为:step(sys)、step(num,den)、step(A,B,C,D)、step(sys,t) 、step(sys1,sys2,…,t) 、y=step(sys,t)、[y,t]=step(sys)、[y,t,x]=step(sys)其中模型对象的类型如下:sys = tf(num,den) 多项式模型、sys = zpk(z,p,k) 零点极点模型、sys = ss(a,b,c,d) 状态空间模型参数无t,表示时间向量t的范围自动设定。

参数有t,表示给定时间向量t,应该有初值,时间增量,末值,如t=0:0.01:2。

前5种函数可以绘出阶跃响应曲线;后3种函数不绘阶跃响应曲线,而是返回响应变量y,时间向量t,以及状态变量x。

2 求取线性连续系统的单位脉冲响应的函数——impulse基本格式为:impulse(sys)、impulse(num,den) 、impulse (sys,tf)、impulse (sys,t) 、impulse (sys1,sys2,…,t) 、y=impulse(sys,t) 、[y,t]=impulse(sys) 、[y,t,x]=impulse (sys)3 求取线性连续系统的单位斜坡响应MATLAB没有直接求系统斜坡响应的功能函数。

在求取控制系统的斜坡响应时,通常用阶跃响应函数step()求取传递函数为G (s)/s的系统的阶跃响应,则其结果就是原系统G (s)的斜坡响应。

实验二

实验二

实验报告(一)实验名称:实验二MATLAB 语言矩阵运算___ 班级:计科1班学号:201440700031 姓名:侯雯娟实验时间:2015年11月9日星期:一得分:批改时间:年月日实验教师(签名):一、实验目的1.掌握建立矩阵的方法;2.掌握MATLAB 各种表达式的书写规则以及常用函数的使用;3.能用matlab 进行基本的数组、矩阵及符号运算;4.掌握矩阵分析的方法以及能用矩阵求逆法解线性方程组。

二、实验环境MATLAB7.0三、实验内容,步骤及结果A =1 2 34 5 6B =2 4 -11 3 5C =1-2D =1 4 78 5 23 6 01、下列运算是否合法,为什么?如合法,结果是多少?(1) result1 = a'result =1 42 53 6(2) result2 = a * b不合法(3) result3 = a + bans =3 6 25 8 11(4) result4 = b * dans =31 22 2240 49 13(5) result5 = [b ; c' ] * dans =31 22 2240 49 13-5 -8 7(6) result6 = a . * bans =2 8 -34 15 30(7) result7 = a . / bans =0.5000 0.5000 -3.00004.0000 1.6667 1.2000 (8) result8 = a . * c不合法(9) result9 = a . \ bans =0.5000 0.5000 -3.00004.0000 1.6667 1.2000 (10) result10 = a . ^2ans =1 4 916 25 36(11) result11 = a ^2不合法(12) result12 = 2 . ^ aans =2 4 816 32 642、使用diag 等函数产生下列矩阵。

(完整版)实验报告第2章参考答案yangh1

meshz(x,y,z),rotate3d
xlabel('X'),ylabel('Y'),zlabel('Z')
pause
figure(2)
contour(x,y,z)
pause
figure(3)
contour3(x,y,z)
思考

习题
结合实验结果及相关理论完成思考及习题内容
1.怎样对隐函数的图形进行显示?
(1) ;
>> syms a b x
>>(3*a*x^2+4*b*x^(1/2))/(x-1)
ans=(3*a*x^2+4*b*x^(1/2))/(x-1)
(2) ;
>> syms x
>> (sin(2*x+pi/4)-log(3*x))/sqrt(x^2+1)
ans=(sin(2*x+1/4*pi)-log(3*x))/(x^2+1)^(1/2)
end
end
end
在命令窗口中运行以下部分:
score =[71 80 87 89 69;85 92 74 94 65;74 90 90 81 92;81 84 82 91 71;70 90 76 70 74;86 85 86 94 85;76 75 81 73 74;91 74 78 73 79;91 75 86 91 67;83 81 84 87 95];
>> a=sqrt(exp(exp(1))); b=tan(pi^2/3);>> a=sqrt(exp(exp(1))); b=tan(pi^2/3);
>> 2*a^2+3*a*b^3-5*a^3*b^5>> sec(atan(a))

MATLAB应用实验指导书结果(终审稿)

M A T L A B应用实验指导书结果公司内部档案编码:[OPPTR-OPPT28-OPPTL98-OPPNN08]MATLAB语言实验指导书中国矿业大学信息与电气工程学院2014年3月实验一 MATLAB 工作环境熟悉及基本运算一、实验目的:熟悉MATLAB 的工作环境,学会使用MATLAB 进行一些简单的运算。

掌握基本的矩阵运算及常用的函数。

二、实验内容:MATLAB 的启动和退出,熟悉MATLAB 的桌面(Desktop ),包括菜单(Menu )、工具条 (Toolbar )、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等;完成一些基本的矩阵操作;学习使用在线帮助系统。

三、实验步骤:1、启动MATLAB ,熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算,观察workspace 的变化,记录运算结果。

(1)(365-522-70)3 =(2)area=pi*^2 =(3)已知x=3,y=4,在MATLAB 中求z :()232y x y x z -== 576 (4)将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1( 2 , 3 )=10>>m1( 11 )=6>>m1( : , 3 )= 3 10 6 15>>m1( 2 : 3 , 1 : 3 )=[ 5 11 10;9 7 6]>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)=34(5)执行命令>>help abs查看函数abs 的用法及用途,计算abs( 3 + 4i )=5(6)执行命令>>x=0::6*pi;>>y=5*sin(x);>>plot(x,y)(7)运行MATLAB 的演示程序,>>demo ,以便对MATLAB 有一个总体了解。

南华大学《MATLAB及应用》实验报告2

核科学技术学院实验报告实验项目名称MATLAB符号计算所属课程名称MATLAB及应用实验类型上机实验实验日期12月日指导教师谢芹班级学号姓名成绩一、实验名称MATLAB符号计算二、实验目的(1)掌握定义符号对象的方法(2)掌握符号表达式的运算法则以及符号矩阵运算(3)掌握求符号函数极限及导数的方法(4)掌握求符号函数定积分和不定积分的方法三、实验原理1. 函数极限及导数的方法(1)函数极限:limit(F,x,a) 求符号函数f(x)的极限值。

即计算当变量x趋近于常数a时,f(x)函数的极限值。

(2)limit(f):求符号函数f(x)的极限值。

符号函数f(x)的变量为函数findsym(f)确定的默认变量;没有指定变量的目标值时,系统默认变量趋近于0,即a=0的情况。

(3)limit(f,x,a,'right'):求符号函数f的极限值。

'right'表示变量x从右边趋近于a。

(4)limit(f,x,a,‘left’):求符号函数f的极限值。

‘left’表示变量x从左边趋近于a。

2. 微分:diff(s):没有指定变量和导数阶数,则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。

diff(s,'v'):以v为自变量,对符号表达式s求一阶导数。

diff(s,n):按findsym函数指示的默认变量对符号表达式s求n阶导数,n为正整数。

diff(s,'v',n):以v为自变量,对符号表达式s求n阶导数。

3. 函数定积分和不定积分的方法:int(s):没有指定积分变量和积分阶数时,系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。

int(s,v):以v为自变量,对被积函数或符号表达式s求不定积分。

int(s,v,a,b):求定积分运算。

a,b分别表示定积分的下限和上限。

梯形法:trapz(x,y):x为分割点构成的向量,y为被积函数在分割点上的函数值构成的向量;抛物线法:quad(f,a,b,tol),f 是被积函数,[a,b]是积分区间,tol 是精度。

MATLAB实验二 傅里叶分析及应用复习课程

M A T L A B实验二傅里叶分析及应用实验二傅里叶分析及应用一、实验目的(一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义2、学会使用Matlab分析周期信号的频谱特性(二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质1、学会运用Matlab求连续时间信号的傅里叶变换2、学会运用Matlab求连续时间信号的频谱图3、学会运用Matlab分析连续时间信号的傅里叶变换的性质(三)掌握使用Matlab完成信号抽样并验证抽样定理1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化3、学会运用MATLAB对抽样后的信号进行重建二、实验条件Win7系统,MATLAB R2015a三、实验内容1、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。

Code:ft = sym('(t+2)*(heaviside(t+2)-heaviside(t+1))+(heaviside(t+1)-heaviside(t-1))+(2-t)*(heaviside(t-1)-heaviside(t-2))');fw = simplify(fourier(ft));subplot(2, 1, 1);ezplot(abs(fw)); grid on;title('amp spectrum');phi = atan(imag(fw) /real(fw));subplot(2, 1, 2);ezplot(phi); grid on;符号运算法Code:dt = 0.01;t = -2: dt: 2;ft = (t+2).*(uCT(t+2)-uCT(t+1))+(uCT(t+1)-uCT(t-1))+(2-t).*(uCT(t-1)-uCT(t-2));N = 2000;k = -N: N;w = pi * k / (N*dt);fw = dt*ft*exp(-i*t'*w);fw = abs(fw);plot(w, fw), grid on;axis([-2*pi 2*pi -1 3.5]);t(20 π ex p(-3 t) heaviside(t) - 8 π ex p(-5 t) heaviside(t))/(2 π)数值运算法2、试用Matlab 命令求ωωωj 54-j 310)F(j ++=的傅里叶反变换,并绘出其时域信号图。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

MATLAB 的绘图功能
姓名 学号 班级___________ 说明:
(1)要求在Notebook 环境下完成;
(2)完成后请以姓名(班级学号)实验报告名称命名并存盘;
(3)在下一次上课之前由学习委员收齐打包以附件形式后发到 lvximing@ ;
(5)请自行保存备份,以备后用。

一.图形修饰
1、在同一个窗口于[]
02π, 做出1sin(2);y x = 2cos(3)y x =图形,并给出标题、图例及有关标注
【程序代码】
x=0:0.01:2*pi;
y1=sin(2*x);y2=cos(3*x);
plot(x,y1,'+r',x,y2,'pb')
title('图1正弦和余弦曲线');
legend ('sin(2*x)','cos(3*x)')
xlabel('横坐标x'),ylabel('纵坐标y')
2、画出y=cos(sin(tanx(x))+3)在[0,2*pi]的极坐标图象【程序代码】
x=0:0.1:2*pi;
y=cos(sin(tan(x))+3);
polar (x,y)
二.特殊图形绘图
x
=
、在[-,]上用蓝色星形线绘制的阶梯图. 114-
y e2
【程序代码】
x=-1:pi/20:4;
y=exp(-x.^2);
stairs(x,y)
、绘制向量的饼图,要求把最大和最小的两部分突出显示2x=[1 4 5 2 6 3 ]
【程序代码】
x=[1,4,5,2,6,3];
y=[1,0,0,0,1,0];
pie(x,y)
3、y x x ππtan(sin )sin(tan )=-在[-,]上绘制的火柴杆图
【程序代码】
x=-pi:pi/20:pi;
y=tan(sin(x))-sin(tan(x));
stem(x,y)
三.从网络或参考资料上收集并整理1~2个与所学知识相关的MATLAB程序并运行其结果(添加必要的注释)
1.绘制彗星曲线图
t= -pi:100:pi;
y=tan(sin(t))-sin(tan(t));
comet(t,y)
2.绘制r=sin(t)cos(t)的极坐标图,并标记数据点。

程序如下:
t=0:pi/50:2*pi;
r=sin(t).*cos(t);
polar(t,r,'-*');。

相关文档
最新文档