贵大matlab实验四
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。
矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。
MATLAB实验报告
MATLAB实验报告实验报告课程名称 MATLAB基础及应用专业班级电子xxxx姓名学号电气与信息学院实验二 MATLAB 数值计算（一）一实验目的：1．掌握数组的创建与运算方法；2. 掌握矩阵的创建与运算方法；3. 掌握数组的运算方法和矩阵运算方法的区别；4．掌握线性方程的求解方法二实验装置：计算机三实验内容：1．数组的创建和运算创建两个含5个元素的一维数组，并求这两个数组的四则运算。
2．矩阵的创建和运算（1）创建两个3×3的矩阵，并求这两个矩阵的四则运算。
（2）创建一个4×4的矩阵，并求这个矩阵的行列式值。
（3）线性方程的求解解方程组=6613753467294x 。
四实验要求:写出实验程序与仿真结果。
1. >> a=linspace(2,10,5) a =2 4 6 8 10>> b=linspace(1,9,5) b =1 3 5 7 9>> a+bans =3 7 11 15 19>> a-bans =1 1 1 1 1>> a.*bans =2 12 30 56 90>> a./bans =2.0000 1.3333 1.2000 1.1429 1.1111 >> a.\b ans =0.5000 0.7500 0.8333 0.8750 0.9000 2.(1) >> a=[1 1 1;2 2 2;3 3 3];>> b=[4 4 4;5 5 5;6 6 6];>> a+bans =5 5 57 7 79 9 9>> a-bans =-3 -3 -3-3 -3 -3-3 -3 -3>> a*bans =15 15 1530 30 3045 45 45>> a/bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNNaN NaN NaN>> a\bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNInf Inf Inf(2).>> a=[1 1 1 1;2 2 2 2;3 3 3 3;4 4 4 4];>> det(a)ans =(3).>> a=[4 9 2;7 6 4;3 5 7];>> b=[13;6;6];>> x=inv(a)*bx =-0.53081.7109-0.1374实验三 MATLAB数值计算（二）一实验目的：1．掌握多项式的创建与运算方法；2. 掌握基本的数据分析方法；二实验装置：计算机三实验内容：1．多项式创建输入系数矢量，创建多项式x^3-2*x^2+5*x+3。
MATLAB实验报告(四份!!!)
0380160815880001449451490168779983287916534749293454515731] >> double(x) ans = -33.5325 60.6078 -54.7098 27.3192 实验三 MATLAB 绘图
一 实验环境 计算机 MATLAB 软件 二 实验目的 1. 掌握 MATLAB 的基本绘命令。 2. 掌握运用 MATLAB 绘制一维、二维、三维的图形方法。 3. 给图形加以修饰。 三 预备知识 特殊的二维图形函数有：直方图、柄图、阶梯图、饼图、频数累计柱状 图、极坐标图 四 实验内容和步骤
3.创建一个 3x3 矩阵，并求其转置，逆矩阵。 >> A=magic(3) A= 8 3 4 >> A' ans = 8 1 6 >> inv(A) 3 5 7 4 9 2 1 5 9 6 7 2
ans = 0.1472 -0.0611 -0.0194 -0.1444 0.0222 0.1889 0.0639 0.1056 -0.1028
3. 画出 y=的曲线（x(-5,5)） 。在这曲线上加入相同区间里的 y=的曲线，并且采用绿色折 线标识。 >> x=linspace(-5,5,100); >> y=x.^2; >> plot(x,y) >> hold on >> z=x.^(1/3); >> plot(x,z,'g--')
4. 在同一窗口不同坐标系里分别绘出 y1=sinx,y2=cosx,y3=cinh(x),y4=cosh(x)4 个图形。 >> x=linspace(0,2*pi,30); >> subplot(2,2,1);plot(x,sin(x)); >> subplot(2,2,2);plot(x,cos(x)); >> subplot(2,2,3);plot(x,sinh(x)); >> subplot(2,2,4);plot(x,cosh(x));
2013MATLAB原理及应用实验报告第四章
《MATLAB原理及应用》实验报告第四章MATLAB程序设计一．实验目的1、掌握脚本文件的建立2、掌握条件语句和程序语句的使用3、掌握MATLAB的程序设计方法二．实验内容1.关于M脚本文件和M函数文件M文件是在M文件编辑器窗口中编写的。
在MATLAB的桌面上单击新建按钮，就可以打开M文件编辑器窗口，也可以通过依次单击【File】/【New】/【M-File】打开文件编辑器。
【实验4-1】脚本文件在M文件编辑器窗口输入一下内容：N=3;for m=1:Nfor n=1:Nif m==nA(m,n)=1;elseA(m,n)=0;endendend单击M文件编辑器窗口中的保存按钮，以“ex1.m”为文件名保存在当前工作目录下。
在命令窗口中输入：>> ex1运行后可以在命令窗口中看到变量A的图标，继续在命令窗口输入：>> AA =1 0 00 1 00 0 1这一脚本文件创建了一个3阶单位阵【实验4-2】 M 函数在M 文件编辑器窗口输入一下内容：function [mean,stdev]=stat(x)x=input('请输入x 的值：')%[mean,stdev]=stat(x)计算输入向量的均值和平均差%输入参数x 是向量%第一个输入参数mean 是向量各元素的平均值%第二个输入参数stdev 是向量的均方差%例如，取向量x=[1,2,3,4,5];%调用[mean,stdev]=stat(x)，计算可得%均值mean=3%均方差stdev=1.4142n=length(x) %计算向量长度mean=sum(x)/n %计算向量平均值stdev=sqrt(sum((x-mean).^2/n)) %计算均方差输入完毕后，单击保存按钮，把文件保存在当前工作目录下，文件名为“stat.m ”2.MATLAB 程序流程控制【实验4-3】一个简单的for 循环事例。
MATLAB实验报告
MATLAB实验报告专业班级学号姓名实验一Matlab环境及命令窗口的使用实验目的与要求：熟练掌握MA TLAB的启动和退出；熟悉MA TLAB的命令窗口；熟悉常用选单和工具栏；熟悉MA TLAB桌面的其它窗口。
实验内容：启动MA TLAB，熟悉各环境窗口及作用，熟悉file、edit、View等菜单下各选项的功能；熟悉工具栏中各项的功能；打开演示程序（demo），开始其中的演示程序；退出MA TLAB（exit、quit）。
在命令窗口中，输入简单的命令操作；打开函数编辑窗，编写简单的程序并保存；程序：>> x1=sqrt(100),x2=5.48,y=9/x2x1 =10x2 =5.4800y =1.6423>>demo演示框左侧是库目录，选定demo下的图形类（Graphics）按钮，双击演示框右边的栏里的function of complex variable选项，就会出现如图中的演示框。
此例为复数Z的三次方：Z^3。
实验二Matlab数值计算实验目的与要求：熟练掌握matlab变量的使用；熟练掌握矩阵的创建；熟悉掌握matlab 的矩阵和数组的运算；熟悉matlab多项式的运算。
实验内容：1、给2×3阶矩阵变量a赋值（注意逗号与分号的使用），取其第2行第2个元素，取其第2行元素，取其第1列和第3列的元素；计算a/3、5*a、a+2等的值；程序：>> a=[7,9,5;5,9,4]a =7 9 55 9 4>> a=[7,9,5;5 9 4]a =7 9 55 9 4>> b=a([2],[2])b =9>> c=a([2],:)c =5 9 4>> d=a(:,[1,3])d =7 55 4>> e=a/3e =2.33333.0000 1.66671.6667 3.0000 1.3333>> 5*aans =35 45 2525 45 20>> a+2ans =9 11 77 11 62、生成2×3阶的全1矩阵x和3×3魔方矩阵y，计算a+x,x*y，x/y，y^2,2^a;a.*x, a.\x，y.^2,2.^a的值。
Matlab实验报告
实验结果及分析实验1：程序如下x=1:10y=2*x;plot(x,y)仿真结果：实验结果分析：仿真结果是条很规则的直线，X轴和Y轴一一对应，清楚明了，而序又特别简单。
所以用Maltab 软件很方便地画出规则的直线，方便研究。
实验结果及分析1、A=2、A=1A=实验结果及分析实验三 Matlab在信号与系统中的应用实验名称实验1、掌握信号与系统课程中基本知识的Matlab编程、仿真方法目的实验原理实验1程序：b=[1];a=[1 1];p=;t=0:p:5;x=exp(-3*t);subplot(1,2,1);impulse(b,a,0:p:5);title('冲激响应');subplot(1,2,2);step(b,a,0:p:5);title('阶跃响应');实验内容<设计性实验>1、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)= exp(-3t)ε(t)的冲激响应、阶跃响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
<设计性实验>（选做）2、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)=(1+exp(-3t))ε(t)的冲激响应、阶跃响应，要求用conv编程实现系统响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
实验结果及分析实验1仿真结果：simulink仿真环境下冲激响应阶跃响应实验名称实验四 Matlab在数字信号处理中的应用实验结果及分析实验1仿真结果:6khz12kHZ。
MATLAB实验报告四
实验四文件操作一、实验目的掌握MATLAB文件操作的各个函数。
二、实验内容1. 验证本章所举的例子例 5.1>> Fid=fopen('std.dat','w');>> Fid=fopen('std.dat','r')Fid =3>> A=fread(Fid);>> Sat=fclose(Fid)Sat =例 5.2>> magic=[1 2 3 4;5 6 7 8]magic =123 45678>> Fid=fopen('magic5.bin','w') Fid =3>> fwrite(Fid,magic,'int32') ans =8>> Sat=fclose(Fid)Sat =>> A=[-0.6515 -0.2727 -0.4354 -0.3190 -0.9047;-0.7534 -0.4567 -0.3212 -0.4132 -0.3583;-0.9264 -0.8173 -0.7823 -0.3265 -0.0631;-0.1735 -0.7373 -0.0972 -0.3267 -0.6298;-0.4768 -0.6773 -0.6574 -0.1923 -0.4389]A =-0.6515-0.2727-0.4354-0.3190-0.9047-0.7534-0.4567-0.3212-0.4132-0.3583-0.9264-0.8173-0.7823-0.3265-0.0631-0.1735-0.7373-0.0972-0.3267-0.6298-0.4768-0.6773-0.6574-0.1923-0.4389>> Fid=fopen('test.dat','w')Fid =3>> cnt=fwrite(Fid,A,'float')25>> fclose(Fid)ans =程序段将矩阵A的数据以二进制浮点数格式写入文件test.dat 中。
matlab实验四
实验四 MATLAB语言绘图一、实验目的：1、学习MATLAB图形绘制的基本方法；2、熟悉和了解MATLAB图形绘制程序编辑的基本指令；3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形的各种标注；4、掌握plot、subplot的指令格式和语法。
二、实验基本知识：1.plot(x,y)------绘制由x,y所确定的曲线；2.图形窗口的分割；3.图形编辑窗口的使用。
三、实验内容1.用图形表示连续调制波形 y=sint*sin(9t)及其包络线。
>> t=(0:pi/100:pi)';>> y1=sin(t)*[1,-1];>> y2=sin(t).*sin(9*t);>> t3=pi*(0:9)/9;>> y3=sin(t3).*sin(9*t3);>> plot(t,y1,'r:',t,y2,'-bo')>> hold on>> plot(t3,y3,'s','MarkerSize',10,' MarkerEdgeColor ',[0 10],'MarkerFaceColor',[1,0.8,0])>> axis([0,pi,-1,1])>> hold off00.51 1.52 2.532.画出t e y t sin 5.0-=的曲线，在图上用“小红圈”标出第一次使y=0.6的那点位置。
>> fun=inline('exp(-0.5*t)*sin(t)-0.6'fun =Inline function:fun(t) = exp(-0.5*t)*sin(t)-0.6>> a=fzero(fun,0)a =-3.2595>> fplot(fun,[-2*pi,2*pi])>> hold on>> plot(a,0.6,'ro')-6-4-20246-4-2246810123.在同一个图形窗里画出实验三中第四题的三个输出，给出图题、横坐标、纵坐标名称。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验四 MATLAB 绘图
一、实验目的
1．掌握二维图形的绘制。
2．掌握图形的标注
3．了解三维曲线和曲面图形的绘制。
二、实验的设备及条件
计算机一台（带有MATLAB7.0以上的软件环境）。
设计提示
1．Matlab 允许在一个图形中画多条曲线：plot(x1，y1，x2，y2，……)指令绘制y 1 = f 1(x 1), y 2 = f 2 (x 2 )等多条曲线。
Matlab 自动给这些曲线以不同颜色。
标注可用text 函数。
2．绘图时可以考虑极坐标和直角坐标的转换。
3．三维曲线绘图函数为plot3，注意参考帮助中的示例。
三、实验内容
1．生成1×10 维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、杆图、阶梯图和条形图，并分别标出标题“连线图”、“杆图”、“阶梯图”、“条形图”。
2、绘制函数曲线，要求写出程序代码。
(1) 在区间[0:2π]均匀的取20个点，构成向量t
(2) 在同一窗口用黑色绘制曲线y1=sin(2*t-0.3); 蓝色绘制曲线小y2=3cos(t+0.5)；
(3)用红色圆圈标注两曲线的交叉点(使用ginput 函数获取交叉点的近似值)
（4）为你所绘制的图形添加相应的图形标注。
3.用mesh 或surf 函数，绘制下面方程所表示的三维空间曲面，x 和y 的取值范围设为[-3，3]。
10102
2y x z +-=
4.在[-2π,2π]的区间上，分别用plot 和fplot 函数绘制函数y=sin(1/x)的曲线，分析两曲线的差别。
思考题：
1.编写一个mcircle(r)函数，调用该函数时，根据给定的半径r，以原点
为圆心,画一个如图所示的红色空心圆。
（图例半径r=5）；左图参考
polar函数的用法，右图绘制圆形的参数方程为x=sin（t），y=cos
（t）。
其中，t的区间为0~2*pi，步长为0.1。
2．（1）绘一个圆柱螺旋线（形似弹簧）图。
圆柱截面直径为10，高度为5，每圈上升高度为1。
如左图所示。
（2）利用（1）的结果，对程序做少许修改，得到如右图所示图形。
四、实验报告要求（包含预习报告要求和最终报告要求）
1.实验名称
2.实验目的
3.实验设备及条件
4.实验内容及要求
5.实验程序设计
指程序代码。
6.实验结果及结果分析
实验结果要求必须客观,有数据的可以记录数据，没有数据的简单描述实验现象。
结果分析是对实验结果的理论评判。
7.实验中出现的问题及解决方法
8. 思考题的回答
五．实验代码
clc
clear all
clf
%1
y=rand(1,10);
x=linspace(0,10,10);
subplot(2,2,1)
plot(x,y,'.-')
subplot(2,2,2)
stem(x,y)
subplot(2,2,3)
stairs(x,y)
subplot(2,2,4)
bar(x,y)
%2
figure
t=linspace(0,2*pi,20);
y1=sin(2*t-0.3);
y2=3*cos(t+0.5);
plot(t,y1,'k',t,y2,'b')
xlabel('t')
ylabel('y')
[x,y] = ginput(2);
hold on
plot(x(1),y(1),'ro',x(2),y(2),'ro')
gtext('y1=sin(2t-0.3)')
gtext('y2=3cos(t+0.5)')
%3
figure
x=-3:0.1:3;
y=-3:0.1:3;
[X,Y]=meshgrid(x,y);
Z=-X.^2/10+Y.^2/10;
mesh(X,Y,Z)
%4
figure
a=linspace(-2*pi,2*pi,100);
b=sin(1./a);
subplot(1,2,1)
plot(a,b)
subplot(1,2,2)
fplot('sin(1/x)',[-2*pi,2*pi])
用fplot函数能够自动在斜率大的点处多取参考点，使图像更准确
思考题：
%1
mcircle函数:
function [y,t]=mcircle(r)
figure
subplot(1,2,1);
t=0:2*pi/64:2*pi;
y=r*ones(size(t));
polar(t,y,'*r')
[X,Y]=pol2cart(t,y);
t = 0:0.1:2*pi;
x=r*sin(t);
y=r*cos(t);
subplot(1,2,2);
axis square;
axis equal;
plot(x,y,'*r')
调用时：
r=5;
mcircle(r);
% 2.(1)
t=0:pi/100:2*pi*5;
r1=10/2;
x1=r1.*cos(t);
y1=r1.*sin(t);
z=t/(2*pi);
figure
subplot(1,2,1);
plot3(x1,y1,z);
%2.(2)
r2=linspace(5,0,length(t));
x2=r2.*cos(t);
y2=r2.*sin(t);
subplot(1,2,2)
plot3(x2,y2,z)
六.实验结果及结果分析。