matlab综合实验4
河北工业大学MATLAB实验四
2014秋2012级《MATLAB 程序设计》实验报告一、实验目的：1、掌握MATLAB 二维图形绘制命令及其图形控制；2、熟悉MATLAB 三维图形绘制命令及其图形控制；3、熟悉特殊二维图形、三维图形的绘制方法。
二、实验内容：1、在同一坐标系下绘制下面三个函数在[0,40]t ∈的图象，要求三种曲线采用不同颜色以及不同的线型，同时对每条曲线利用函数进行标注，并添加图例。
120.134sin()t y t y y e t π-===2、编写程序，选择合适的步距，绘制下面函数在区间[-6,6]中的图象，并对分段的曲线进行标注，同时添加x 轴和y 轴的说明。
sin ,0(),036,3x x y x x x x x ≤⎧⎪=<≤⎨⎪-+>⎩3、利用subplot 函数在同一绘图窗口中用不同颜色和线型绘制以下两个函数在t ∈[-2π,2π]范围内的图象。
0.50.21222t t y y e -==4、某学校有30位教师,其中教授5人,副教授8人,讲师12人,助教5人.试画出职称比例结构饼图，并强调图形的教授部分同时添加标注。
5、生成1×10维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、脉冲图、阶梯图和条形图，并分别标出标题“连线图”、“脉冲图”、“阶梯图”、“条形图”。
6、使用subplot 函数，把图形窗口分成两个部分，分别绘制sin2t 和3cos3t 曲线,t 范围:0-5.要求前者用红色实线，数据点形状为五角星,有网格线，x 轴加标注“x ”，y 轴加标注“y=sin2t ”,加题目“y=sin2t 的曲线”；后者用蓝色虚线，数据点形状为圆形，无网格线，x 轴加标注“x ”，y 轴加标注“y=3cos3t ”,加题目“y=3cos3t 的曲线”；7、绘制下列三维曲线:(1)/20/20cos sin ,02t t x e ty e t t z t π--⎧=⎪=≤≤⎨⎪=⎩ (2) 23,01x t y t t z t =⎧⎪=≤≤⎨⎪=⎩8、绘制下列曲面图，并调整三维图的视角、背景色、着色以及透视效果。
matlab实验3 4 5
实验四MATLAB 程序设计一、实验目的掌握MATLAB 程序设计的主要方法，熟练编写MATLAB 函数。
二、实验内容(1) M 文件的编辑。
(2) 程序流程控制结构。
(3) 子函数调用和参数传递。
(4) 局部变量和全局变量。
三、实验结果练习：(1) 请把exp2.m 函数文件用while 循环改写。
(2) 用π/4≈1-1/3+1/5-1/7+…公式求π的近似值，直到最后一项的绝对值小于10-6为止，试编写其M 脚本文件。
(1)解：x=input('输入一个自然数:')s=0n=1while n<=xs=s+nn=n+1end(2)解：s=1;k=1;n=1;while(k*(-1)^(n+1)>10^(-6))k=k*(1-2*n)/(1+2*n);s=s+k;n=n+1;endpai=s*4练习：(1) 编写求矩形面积函数rect，当没有输入参数时，显示提示信息；当只输入一个参数时，则以该参数作为正方形的边长计算其面积；当有两个参数时，则以这两个参数为长和宽计算其面积。
(2) 编写一个字符串加密函数nch=my_code( ch , x)，其中ch 是字符串参数，x 为整数；加密方法是：把ch 的每一个字符的ASCII 码值加上x，得到的即为加密后的新的字符串nch。
由于可显示ASCII 码值是有范围的(32，126)，因此当得到的ASCII码值大于126 时，需要减去93 。
同理，再编写一个解码函数nch=my_dcode( ch , x)。
提示：char(32:126) 可获得ASCII 码值为32～126 的字符。
7(1)解：function s=rect(a,b)switch(nargin)case 0error('请输入一个或两个数');case 1s=a*a;case 2s=a*b;end(2)加密：function nch=my_code(ch,x)nch=ch+x;if nch>126nch=nch-93;char(nch)elsechar(nch)end解码：function nch=my_code(ch,x)nch=ch-x;if nch<32nch=nch+93;char(nch)elsechar(nch)end实验五MATLAB 数据可视化一、实验目的掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。
matlab实验三、四、五、六(2012-2013-1)
实验三 图像的代数运算一、实验目的1．了解图像的算术运算在数字图像处理中的初步应用。
2．体会图像算术运算处理的过程和处理前后图像的变化。
二、实验原理图像的代数运算是图像的标准算术操作的实现方法，是两幅输入图像之间进行的点对点的加、减、乘、除运算后得到输出图像的过程。
如果输入图像为A(x,y)和B(x,y)，输出图像为C(x,y)，则图像的代数运算有如下四种形式：C(x,y) = A(x,y) + B(x,y)C(x,y) = A(x,y) - B(x,y)C(x,y) = A(x,y) * B(x,y)C(x,y) = A(x,y) / B(x,y)图像的代数运算在图像处理中有着广泛的应用，它除了可以实现自身所需的算术操作，还能为许多复杂的图像处理提供准备。
例如，图像减法就可以用来检测同一场景或物体生产的两幅或多幅图像的误差。
使用MATLAB的基本算术符(+、-、*、/ 等)可以执行图像的算术操作，但是在此之前必须将图像转换为适合进行基本操作的双精度类型。
为了更方便地对图像进行操作，MATLAB图像处理工具箱包含了一个能够实现所有非稀疏数值数据的算术操作的函数集合。
下表列举了所有图像处理工具箱中的图像代数运算函数。
表2-1 图像处理工具箱中的代数运算函数函数名功能描述Imabsdiff 两幅图像的绝对差值Imadd 两幅图像的加法Imcomplement 补足一幅图像Imdivide 两幅图像的除法Imlincomb 计算两幅图像的线性组合Immultiply 两幅图像的乘法imsubtract 两幅图像的减法使用图像处理工具箱中的图像代数运算函数无需再进行数据类型间的转换，这些函数能够接受uint8和uint16数据，并返回相同格式的图像结果。
虽然在函数执行过程中元素是以双精度进行计算的，但是MATLAB工作平台并不会将图像转换为双精度类型。
代数运算的结果很容易超出数据类型允许的范围。
例如，uint8数据能够存储的最大数值是255，各种代数运算尤其是乘法运算的结果很容易超过这个数值，有时代数操作（主要是除法运算）也会产生不能用整数描述的分数结果。
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实验报告(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是非常重要的工具。
本课程设计主要目的是让学生通过实验，掌握使用MATLAB进行信号与系统分析和处理的方法和技巧。
同时，课程设计还能够加深学生对信号与系统理论知识的理解和掌握，提高其综合运用能力。
二、课程设计的内容和要求1. 实验一：信号的生成和绘制本实验主要包括以下内容：•生成几种基本信号（如正弦信号、方波信号、三角波信号等）。
•通过MATLAB绘制生成的信号，并加上合适的标注。
要求学生能够掌握信号的生成方法和MATLAB的绘图函数的使用。
2. 实验二：信号的运算与变换本实验主要包括以下内容：•对已有信号进行运算（如加、减、乘、除等）。
•对信号进行卷积、相关等线性变换操作。
•对信号进行傅里叶变换，并绘制幅度谱、相位谱等图形。
要求学生能够掌握信号的运算、变换方法和MATLAB的相应函数的使用。
3. 实验三：系统的分析和建立本实验主要包括以下内容：•对系统进行零极点分析，并绘制零极点图。
•对已有系统进行时域和频域分析（如阶跃响应、冲击响应、幅频响应等）。
要求学生能够掌握系统的分析方法和MATLAB的相应函数的使用。
4. 实验四：信号的滤波和降噪本实验主要包括以下内容：•对信号进行数字滤波（如低通滤波、高通滤波、带通滤波、带阻滤波等）。
•对信号进行去噪处理（如中值滤波、小波变换去噪等）。
要求学生能够掌握信号滤波、降噪方法和MATLAB的相应函数的使用。
三、课程设计的实施流程1.分组。
依据班级人数以及教学设备的数量，安排学生分为若干个小组，每个小组3-4人。
2.模拟分配实验。
询问小组成员们的意见，模拟分配每个小组所要完成的课程设计任务。
3.实验操作。
每个小组根据分配到的实验课程设计，使用MATLAB进行模拟操作。
4.结果展示。
每个小组进行结果展示，介绍自己的设计思路，并展示实验结果。
其他小组成员以及教师进行现场互相交流和讨论。
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.在同一个图形窗里画出实验三中第四题的三个输出，给出图题、横坐标、纵坐标名称。
matlab实验四
（理工类）实验报告书写要求实验报告原则上要求学生手写，要求书写工整。
若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。
纸张一律采用A4的纸张。
实验报告书写说明实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。
各院部可根据学科特点和实验具体要求增加项目。
填写注意事项（1）细致观察，及时、准确、如实记录。
（2）准确说明，层次清晰。
（3）尽量采用专用术语来说明事物。
（4）外文、符号、公式要准确，应使用统一规定的名词和符号。
（5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。
实验报告批改说明实验报告的批改要及时、认真、仔细，一律用红色笔批改。
实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。
实验报告装订要求实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：MATLAB的符号运算实验学时： 4同组学生姓名：实验地点： C304实验日期： 2012.5.31 6.7 实验成绩：批改教师：批改时间：一、实验目的和要求1、会用MA TLAB求电阻电路。
（节点电压法、戴维南定理等）2、会用MA TLAB求正弦稳态电路。
3、掌握自动控制系统中模型的转换方法。
4、会用MA TLAB求自动控制系统的传递函数。
二、实验仪器和设备计算机一台三、实验过程1、已知条件如下图所示。
求:RL为何值时,能获得最大功率.R L解：分析：把RL用电流源代替（Ia方向从上到下），列节点电压方程，取O为参考点。
经过电路化简，得出关系：Ia=4/3+I2。
列出节点方程：（1/6+1/3+1）U1-U2=I1 {其中I1=12/6=2(A);I2=2(A)}-U1+U2=I2-Ia运用matlab计算出U1,U2的值。
U2即是戴维南等效电路中的Uoc。
令电压源短路，电流源断路，求得ac左边等效电阻Req。
matlab实验报告(实验4)
学生实验报告开课学院及实验室: 机电学院2012年12月21日学院机电学院年级、专业、班姓名学号实验课程名称MATLAB程序设计成绩实验项目名称实验4: 数据和函数的可视化指导老师一、实验目的1、掌握MATLAB绘图的基本步骤和相关指令调用的先后顺序。
2、掌握MATLAB绘图指令的调用方法。
二、实验内容数学函数从形式上可以分为离散函数和连续函数。
MATLAB对这两种函数数据的可视化都提供了相应的指令。
仔细阅读教材【例5.1-1】的实现代码, 运行并保存结果；并改用stem函数, 画出【例5.1-1】的序列图。
仔细阅读教材【例5.1-2】的实现代码, 运行并保存结果；并分别使用描点和连折线方式, 画出连续函数y=xcosx的近似图形（采样点数自定, 要求画出的图尽量接近原连续函数的图）。
仔细阅读【例5.2-2】的实现代码, 理解plot指令画多条曲线的运用方法, 运行并保存结果；并使用plot函数和legend函数, 在同一个图形窗口上画出y=sint和y=sin(2t)在[0,2pi]区间上的图形, 并标出图例。
仔细阅读【例5.2-4】的实现代码, 理解图形标识选项的运用方法, 运行并保存结果；并修改代码, 把“sin(t)”字体改为正体, 大小改为20, “极大值”改为宋体。
阅读【例5.2-6】, 理解使用hold on指令画多幅图的方法, 运行并保存结果。
阅读【例5.2-8】, 理解使用subplot函数画多个子图的方法, 运行并保存结果。
（1）综合实验: 阅读以下关于通过绘制二阶系统阶跃响应综合演示图形标识的示例, 理解示例中所有图形标识指令的作用, 掌握各个图形标识指令的运用方法, 并在原指令上改动以实现以下功能:（2）把横坐标范围改为0至5pi, 纵坐标范围改为0至2；（3）把图中的横轴的刻度改为从0开始到4pi, 中间各点间隔为pi/2；纵轴刻度改为从0开始到1.5, 中间各点间隔为0.3；（4）把图中的α改为σ。
matlab实验四
实验四Matlab图形绘制一、实验目的：1、熟练掌握MATLAB二维曲线的绘制。
2、熟练掌握图形的修饰。
3、掌握三维图形的绘制。
二、实验内容与步骤：MATLAB的图形功能非常大，可以对二维、三维数据用图形表现，并可以对图形的线型、曲面、视角、色彩、光线等进行处理。
MATLAB的语法规则简洁，编程效率高，作为一个完整的程序语言，MATLAB也有其各种程序流程控制、文件格式和函数调用的规则，通过函数的调用就能够组成庞大的程序，完成复杂的功能。
1.绘制二维曲线（1）plot(x) %绘制以x为纵坐标的二维曲线plot(x,y) %绘制以x为横坐标、y为纵坐标的二维曲线>>x=linspace(0,2*pi,100);>>plot(sin(x))（2）用plot(x1,y1,x2,y2,...)命令绘制多条曲线>> x=0:0.1:2*pi;>> plot(1,sin(x),x,cos(x),x,3*cos(x))2.同一个窗口多个子图及图形处理绘制图形窗口分割为2列2行，窗口左上角一正弦曲线，y=sin(2*pi*t),t∈[0,2];窗口右上角画3条衰减的单边指数曲线y=e^-t,y=e^-2t和y=e^-3t,t∈[0,2];窗口左下角画一矩形脉冲信号，脉冲宽度为1，高度为2，开始时间为1；窗口右下角画一单位圆。
MATLAB允许在同一窗口绘制多个子图，使用subplot命令，各子图的顺序是先向右后向下。
①2行2列子图的第1个图;>>subplot(2,2,1)>>t1=0:0.1:2;>>y1=sin(2*pi*t1);>>plot(t1,y1);>>title(‘y=sin(2*pi*t)’) %添加标题②2行2列字图的第2个图：>>subplot(2,2,2)>>t2=0:0.1:2;>>y2=[exp(-t2);exp(-2*t2);exp(-3*t2)];>>plot(t2,y2)>>axis([0 2 -0.1 1.2]); %设定坐标范围>>title(‘y=e-t,y=e-2t,y=e-3t’)③2行2列子图的第3个图：>>subplot(2,2,3);>>t3=[0 1 1 2 2 3 4]>>y3=[0 0 2 2 0 0 0]>>plot(t3,y3, ':*r'); %点线、星号、红色线输出>>axis([0 4 -0.5 3]);>>title(‘脉冲信号’)④2行2列子图的第4个图：>>subplot（2，2，4);>>t4=0:0.1:2*pi;>>plot（sin(t4),cos(t4)，'--og'); %虚线、圆。
MATLAB第四章实验报告
MATLAB图形基础班级：11电信一班姓名：何得中学号：20111060108实验目的：1.掌握MATLAB图形绘制功能；2.学会使用MATLAB提供的绘图函数；3.了解MATLAB对图形线性、色彩、光线、视角等的制定和处理。
实验仪器：MA TLAB软件，电脑实验数据：绘制fracreturnl变量同一图轴下绘制两条曲线增加子图分割子图删除图轴后仍保留两个子图效果添加图轴信息设定曲线样式在交互式绘图工具中增加图例>> rand('state',27)>> startprice=50;>>fracreturns1=0.0015*randn(200,1)+0.0003; >> x1=[startprice;1+fracreturns1];>> prices1=cumprod(x1); >> t=(1:length(prices1))'; >> randn('state',7) >>fracreturns2=0.0015*randn(200,1)+0.0003; >> x2=[startprice;1+fracreturns2]; >> prices2=cumprod(x2); >> whos Name Size Bytes Class L 51x51 20808 double array fracreturns1 200x1 1600 double array fracreturns2 200x1 1600 double array l1 1x1 8 double array l2 1x1 8 double array prices1 201x1 1608 double array prices2 201x1 1608 double array s 1x1 8 double arraystartprice 1x1 8 double array t 201x1 1608 double array x1 201x1 1608 double array x2 201x1 1608 double array Grand total is 4010 elements using 32080 bytes>> plot (prices2, 'DisplayName', 'prices2','YDataSource', 'prices2'); figure(gcf)>> x=0:pi/1000:2*pi; >> y=exp(-0.5*x).*sin(2*pi*x+pi/4);>> plot(x,y) plot(x,y,x,y+1,x,y-1) >> B=pascal(6) B = 1 1 1 1 1 1 1 2 3 4 5 6 1 3 6 10 15 21 1 4 10 20 35 56 1 5 15 35 70 126 1 6 21 56 126 252 >> plot(B)>> x=(0:pi/100:2*pi)'; >> y1=2*exp(-0.5*x)*[1,-1]; >> y2=2*exp(-0.5*x).*sin(2*pi*x); >> x1=(0:12)/2; >> y3=2*exp(-0.5*x1).*sin(2*pi*x1);>> plot(x,y1,'g:',x,y2,'b--',x1,y3,'rp');>> x1=0:pi/100:2*pi; >> x2=0:pi/100:3*pi;>> y1=exp(-0.5*x1).*sin(2*pi*x1); >> y2=1.5*exp(-0.1*x2).*sin(x2); >> plotyy(x1,y1,x2,y2);>> x=(0:pi/100:2*pi)';y1=2*exp(-0.5*x)*[1,-1];y2=2*exp(-0.5*x).*sin(2*pi*x);x1=(0:12)/2;y3=2*exp(-0.5*x1).*sin(2*pi*x1);plot(x,y1,'g:',x,y2,'b--',x1,y3,'rp');>> title('曲线及其包络线'); %加图形标题>> xlabel('independent variable X'); %加X轴说明>> ylabel('independent variable Y'); %加y轴说明>> text(2.8,0.5,'包络线'); %在指定位置添加图形说明>> text(0.5,0.5,'曲线y');>> text(1.4,0.1,'离散数据点');>> legend('包络线','包络线','曲线y','离散数据点') %加图例>> x=(0:pi/100:2*pi)';>> y1=2*exp(-0.5*x)*[1,-1];>>y2=2*exp(-0.5*x).*sin(2*pi*x);>> plot(x,y1,'b:');>> axis([0,2*pi,-2,2]);>> hold on;>> plot(x,y2,'k');>> grid on;>> box off;>> hold off;>> x=(0:pi/100:2*pi)';>> y1=2*exp(-0.5*x)*[1,-1];>> y2=2*exp(-0.5*x).*sin(2*pi*x);>> plot(x,y1,'b:');>> hold on;>> plot(x,y2,'k');>> grid on;>> box off;>> hold off>>alpha=-0.5;>>beta=3;>>A=50;>>t=0:0.01:10;>>y=A*exp(alpha*t).*sin(beta*t);>>plot(t,y);>>title('\fontname{隶书}\fontsize{16}{隶书}\fontname{Impact}{Impact}')>>xlabel('^{上标}and_{下标}')>>ylabel('Some\bf粗体\rm and some\it{斜体}')>>txt={'y={\itAe}^{\alphax}sin(\beta\itt)',... >> ['\itA\rm','='.num2str(A)],...>> ['\alpha='num2str(alpha)],...>> ['\beta=',num2str(beta)];>> text(2,22,txt)x=linspace(0,2*pi,60);y=sin(x);z=cos(x);t=sin(x)./(cos(x)+eps);ct=cos(x)./(sin(x)+eps); subplot(2,2,1);plot(x,y);title('sin(x)');axis([0,2*pi,-1,1]); subplot(2,2,2);plot(x,z);title('cos(x)');axis([0,2*pi,-1,1]);subplot(2,2,3);plot(x,t);title('tangent(x)');axis([0,2*pi,-40,40]);subplot(2,2,4);plot(x,ct);title('cotangent(x)');axis([0,2*pi,-40,40]);x=linspace(0,2*pi,60); y=sin(x);z=cos(x);t=sin(x)./(cos(x)+eps);ct=cos(x)./(sin(x)+eps); subplot(2,2,1);stairs(x,y);title('sin(x)-1');axis([0,2*pi,-1,1]); subplot(2,1,2);stem(x,y);title('sin(x)-2');axis([0,2*pi,-1,1]); subplot(4,4,3);plot(x,y);title('sin(x)');axis([0,2*pi,-1,1]); subplot(4,4,4);plot(x,z);title('cos(x)');axis([0,2*pi,-1,1]); subplot(4,4,7);plot(x,t);title('tangent(x)');axis([0,2*pi,-40,40]);subplot(4,4,8);plot(x,ct);title('cotangent(x)');axis([0,2*pi,-40,40]);x=0:0.35:7;y=2*exp(-0.4*x);subplot(2,2,1);bar(x,y,'g');title('bar(x,y,''g'')');axis([0,7,0,2]); subplot(2,2,2);fill(x,y,'r');title('fill(x,y,''r'')');axis([0,7,0,2]); subplot(2,3,4);stairs(x,y,'b');title('stairs(x,y,''b'')');axis([0,7,0,2]);subplot(2,3,5);stem(x,y,'k');title('stem(x,y,k'')');axis([0,7,0,2]); subplot(2,3,6);area(x,y);title('area(x,y)');axis([0,7,0,2]);theta=0:0.01:2*pi;rho=sin(3*theta).*cos(3*theta);polar(theta,rho,'r');x=0:0.1:10; y=12*x.*x;subplot(2,2,1);plot(x,y);title('plot(x,y)');grid on;subplot(2,2,2);semilogx(x,y);title('semilogx(x,y)');grid on;subplot(2,2,3);semilogx(x,y);title('semipogx(x ,y)');grid on;subplot(2,2,4);loglog(x,y);title('loglog(x,y)');grid on;function y=myf(x); y=cos(tan(pi*x));fplot('myf',[-0.2,1.2],1e-4);subplot(1,2,1); pie([7,17,23,19,5]);title('饼图');legend('优秀','良好','中等','及格','不及格');subplot(1,2,2);compass([3+2i,4.5-i,-1.5+5i]);title('相量图');t=0:pi/50:2*pi;x=8*cos(t);y=4*sqrt(2)*sin(t);z=-4*sqrt(2)*si n(t);plot3(x,y,z,'p');title('Line in 3-D space');text(0,0,0,'origin');xlabel('X'),ylabel('Y'),zlabel('Z');grid;x=5:29;y=14:35; [x,y]=meshgrid(x,y); z=2*x+5*y; k=find(z==126); x(k)',y(k)' >> bans =8 13 18 23 28 ans =22 20 18 16 14x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);mesh(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zla bel('z-axis');title('mesh');x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x);surf(x,y,z);xlabel('x-asis'),ylabel('y-axis'),zlabel('z-axis');title('surf');x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);plot3(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('plot3');grid;[x,y]=meshgrid(-10:0.5:10);z=sin(sqrt(x.^3+y.^2))./sqrt(x.^2+y.^2+eps); subplot(2,2,1); meshc(x,y,z);title('meshc(x,y,x)'); subplot(2,2,2); meshz(x,y,z);title('meshz(x,y,x)'); subplot(2,2,3); surfc(x,y,z);title('surfc(x,y,z)'); subplot(2,2,4); surfl(x,y,z);title('surfl(x,y,z)');t=0:pi/20:2*pi;[x,y,z]=cylinder(2+sin(t),30); subplot(1,3,1); surf(x,y,z);subplot(1,3,2); [x,y,z]=sphere; surf(x,y,z); subplot(1,3,3); [x,y,z]=peaks(300);meshz(x,y,z);subplot(2,2,1); bar3(magic(5)); subplot(2,2,2);y=3*sin(0:pi/10:2*pi); stem3(y);subplot(2,2,3);pie3([2347,1827,2043,3025]); subplot(2,2,4);fill3(rand(3,6),rand(3,6),rand(3,6),'r');subplot(1,2,1); [x,y,z]=peaks; waterfall(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis'); title('waterfall of peaks'); subplot(1,2,2);contour3(x,y,z,12,'k');xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('contour3 of peaks');subplot(2,2,1);mesh(peaks); view(-37.5,30);title('azimuth=-37.5,elevation=30') subplot(2,2,2);mesh(peaks); view(0,90);title('azimuth=0,elevation=90') subplot(2,2,3);mesh(peaks); view(90,0);title('azimuth=90,elevation=0') subplot(2,2,4);mesh(peaks); view(-7,-10);title('azimuth=-7,elevation=-10')z=peaks(20);colormap(copper); subplot(1,3,1);surf(z);subplot(1,3,2);surf(z);shading flat;subplot(1,3,3);surf(z);shading interp;x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);[I,J]=find(z>0.2); for ii=1:length(I)z(I(ii),J(ii))=NaN; endsurf(x,y,z);。
matlab 实验报告
matlab 实验报告Matlab实验报告引言：Matlab是一种强大的数值计算和可视化软件，广泛应用于科学、工程和经济等领域。
本实验报告将介绍我在使用Matlab进行实验过程中的一些经验和结果。
实验一：矩阵运算在这个实验中，我使用Matlab进行了矩阵运算。
首先，我创建了一个3x3的矩阵A和一个3x1的矩阵B，并进行了矩阵相乘运算。
通过Matlab的矩阵乘法运算符*，我得到了一个3x1的结果矩阵C。
接着，我对矩阵C进行了转置操作，得到了一个1x3的矩阵D。
最后，我计算了矩阵C和矩阵D的点积，并将结果输出。
实验二：数据可视化在这个实验中，我使用Matlab进行了数据可视化。
我选择了一组实验数据，包括时间和温度两个变量。
首先，我将数据存储在一个矩阵中，并使用Matlab的plot函数将时间和温度之间的关系绘制成曲线图。
接着，我使用Matlab的xlabel、ylabel和title函数添加了横轴、纵轴和标题。
最后，我使用Matlab的legend函数添加了图例，以便更好地理解图表。
实验三：数值积分在这个实验中，我使用Matlab进行了数值积分。
我选择了一个函数f(x)进行积分计算。
首先，我使用Matlab的syms函数定义了符号变量x，并定义了函数f(x)。
接着，我使用Matlab的int函数对函数f(x)进行积分计算，并将结果输出。
为了验证结果的准确性，我还使用了Matlab的diff函数对积分结果进行了求导操作，并与原函数f(x)进行了比较。
实验四：信号处理在这个实验中，我使用Matlab进行了信号处理。
我选择了一个音频文件，并使用Matlab的audioread函数读取了该文件。
接着，我使用Matlab的fft函数对音频信号进行了傅里叶变换，并将结果绘制成频谱图。
为了进一步分析信号的特征，我还使用了Matlab的spectrogram函数绘制了信号的时频图。
通过对信号的频谱和时频图的观察，我可以更好地理解信号的频率和时域特性。
MATLAB实验四：MATLAB程序设计 参考答案
function y=f(x) y=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14)
河南财经政法大学数学与信息科学学院
6
实验报告
f(1)*f(2)+f(3)
function f=f(a,b,c) f=g(a)*g(b)+g(c)^2; function g=g(x) g=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14);
a=rand(1,100); [b,c]=sort(a); index=c(1); i=1; while i<=7 if i==8 disp('You Lost!'); i=i+1; else d=input('please input the number you guess:'); if d>index disp('High'); elseif d<index disp('Low'); elseif d==index disp('You won!'); i=9;
syms i j k l n i=1;k=0; n=input('Input n:'); for j=1:n i=i*j; k=k+j; end fprintf('%d！=%d \n Sum(1,...,%d)=%d\n',n,i,n,k);
4、用 while-end 循环语句求不超过 1000 的偶数之和，并求显示出最大值。 s=0,n=0; while s<=1000 n=n+2; s=s+n; if s>1000 break end
北科大matlab第4次实验报告
《数学实验》报告实验名称数学实验学院自动化专业班级姓名学号2015年4月一、 【实验目的】使用MATLAB 进行线性代数相关运算及求数值积分 二、 【实验任务】1、求下列线性方程组的通解: (1)1234123123420302340x x x x x x x x x x x ++-=⎧⎪-++=⎨⎪-+-=⎩ (2)123412341234031231/2x x x x x x x x x x x x --+=⎧⎪-+-=⎨⎪--+=-⎩2、用三种方法求积分的数值解:20sin 1cos x xdx xπ+⎰3、用多种数值方法计算定积分411sin dx xπ-⎰,,观察不同方法相应的误差.三、 【实验程序】1、 (1)A = [1 1 2 -1; -1 1 3 0; 2 -3 4 -1]; rref(A)(2)A = [1 -1 -1 1 0; 1 -1 1 -1 1; 1 -1 -2 3 -1/2]; rref(A)2、建立函数文件integration.m function y = integration(x) y = x.*sin(x)./(1+cos(x).^2); 编程如下：x = 0:pi/100:pi;y = x.*sin(x)./(1+cos(x).^2); t = length(x);s1 = sum(y(1:(t-1)))*pi/100 %矩形公式7-1 s2 = sum(y(2:t))*pi/100 %矩形公式7-2 s3 = trapz(x,y) %梯形公式7-5 s4 = quad('integration',0,pi) %辛普生公式7-63、建立函数文件integration.mfunction y = integration(x) y = 1./(1-sin(x)); 编程如下：x = 0 : pi/400 : pi/4; y = 1./(1-sin(x)); t = length(x);s1 = sum(y(1 : (t - 1)))*pi/400 %矩形公式7-1 s2 = sum(y(2 : t))*pi/400 %矩形公式7-2 s3 = trapz(x,y) %梯形法7-5 s4 = quad('integration',0,pi/4) %辛普生法7-6 vpa(sqrt(2))四、 【实验结果】 1、 （1）ans =1 0 0 -14/25 0 1 0 -1/5 0 0 1 -3/25A 的秩为3，方程组有无穷多解同解方程组：⎪⎪⎪⎩⎪⎪⎪⎨⎧===425351251434241x x x x x x取1x 4=，解得方程组的基础解系为：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=1253512514ξ所以方程组的通解为：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡12535125144321k x x x x ，其中k 为任意实数（2）ans =1 -1 0 0 1/2 0 0 1 0 1/2 0 0 0 1 0增广矩阵的秩为3，等于系数矩阵的秩，而小于未知量的个数4，所以方程组有无穷多解，原方程组对应的同解方程组为：⎪⎪⎪⎩⎪⎪⎪⎨⎧==+=021221431x x x x 可找到其中一个特解为：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=*021021η再求解对应的齐次方程组⎪⎩⎪⎨⎧===002431x x xx ，可得到一个基础解系：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=0011ξ因此，此方程组的通解为：()R k k x x x x ∈⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡021021001143212、s1 =2.4673 s2 =2.4673 s3 =2.4673 s4 =2.4674矩形公式以及梯形法计算结果都是2.4673，辛普生法计算结果是2.46743、s1 =1.4048 s2 =1.4237s3 =1.4143s4 =1.4142ans =1.4142135623730950488016887242097矩形公式7-1计算，结果是1.4048；矩形公式7-2计算，结果是1.4237；梯形法计算，结果是1.4143；辛普生法计算结果是1.4142。
MATLAB基础性实验报告4
学生实验报告一、实验目的1. 归纳和学习求解常微分方程(组)的基本原理和方法；2. 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析；3. 熟悉MATLAB软件关于微分方程求解的各种命令；4. 通过范例学习建立微分方程方面的数学模型以及求解全过程；通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。
这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。
二、实验仪器、设备或软件：电脑,MATLAB软件三、实验内容1．微分方程及方程组的解析求解法；2．微分方程及方程组的数值求解法——欧拉、欧拉改进算法；3．直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)；4．利用图形对解的特征作定性分析；5．建立微分方程方面的数学模型，并了解建立数学模型的全过程。
四、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口；2．根据微分方程求解步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．根据观察到的结果和体会写出实验报告。
五、实验要求与任务1．求微分方程的解析解，并画出它们的图形。
y '= y + 2 x, y (0) = 1, 0< x <1；解：命令为：dsolve('Dy=y+2*x','y(0)=0','x') 结果为：ans =2*exp(x) - 2*x – 2绘制图形程序：x=linspace(0,1,50);y=2*exp(x) - 2*x - 2;plot(x,y)图形如下：2．求微分方程⎪⎩⎪⎨⎧====-+]100[0)0(;0)0(01.03t uu u u u 的数值解，要求编写求解程序。
MATLAB实验报告_4
《数学实验》综合实验报告报告成员（2人）:李庆超自动化1301 41321038陈吉采矿理留1401 614020152015年 6月运用MATLAB求解电路问题电路问题主要分为稳态电路与动态电路, 这里结合电路原理中的节点电压法, 通过MATLAB求解稳态电路；结合电路原理中的状态方程, 通过MATLAB求解动态电路。
§1运用MATLAB求解稳态电路一、问题提出:求解如图所示的稳态电路, 求出各电阻、电感、电容两端的电压及流过他们的电流和有功功率和功率因数。
（w=100rad/s;Is3=1A;Is4=2A;L1=10mH;L2=50mH;R3=5Ώ;R4=2Ώ;R5=10Ώ;C6=0.1mF）模型建立:本题首先可以运用电路原理的知识, 列出节点电压方程, 将电路结构及其参数用矩阵形式表达出来, 然后运用MATLAB 的知识对其进行求解。
根据电路原理知识, 任何稳态电路的求解可以根据以下公式列节点电压方程：s n TI A U AYA ..=其中, A 为电路的关联矩阵, Y 为支路导纳矩阵, 为节点电压, 为电源列向量。
对本题而言:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=011100001010011101A ⎥⎦⎤⎢⎣⎡=65432111111jwC R R R jwL jwL diag Y⎥⎦⎤⎢⎣⎡=0000.4.3.s s s I I I因此, 在根据 , 可以运用MATLAB 求解出本题的节点电压矩阵: , 进而得到各处节点相对于参考节点的电压。
问题求解:程序代码:w=100;Is3=1;Is4=2;L1=0.01;L2=0.05;R3=5;R4=2;R5=10;C6=0.0001;A=[1 0 1 1 0 0;-1 1 0 0 0 1;0 -1 0 -1 1 0];Is=[0;0;Is3;Is4;0;0];y=[1/(j*w*L1),1/(j*w*L2),1/R3,1/R4,1/R5,1/(j*w*C6)];Y=diag(y,0)Un=inv(A*Y*A')*A*IsIR3=Un(1)/R3IR4=(Un(1)-Un(3))/R4IR5=Un(3)/R5IC6=Un(2)*j*w*C6IL1=(Un(1)-Un(2))/(j*w*L1)IL2=(Un(3)-Un(2))/(j*w*L2)求解结果:Un =0.7705 + 1.1223i0.0023 + 0.0112i-2.6991 + 0.0347iIR3 =0.1541 + 0.2245iIR4 =1.7348 + 0.5438iIR5 =-0.2699 + 0.0035iIC6 =-1.1158e-04 + 2.2793e-05iIL1 =1.1111 - 0.7682iIL2 =0.0047 + 0.5403i结论:通过MATLAB求解得：相对于参考节点4, 节点1、2、3的节点电压分别为（0.7705+j1.223）V,（0.0023+j0.0112）V, （-2.6991+j0.0347）V。
信号与系统 matlab 综合 实验 第四章 傅里叶变换
1．如图4.4所示锯齿波信号，分别取一个周期的抽样数据X1（t），0<=t<=1和五个周期的数据X（t），0<=t<5,计算其傅立叶变换X1（w）和X（w），比较有和不同并解释原因。
图4.4 练习题2图编程如下：%计算单位锯齿波和五个周期波形的傅立叶变换%数值算法用矩阵实现，大大加快了运行速度；并且直接调用“sawtooth”生成5个周期的锯齿波T1=1; %单个周期时域范围N1=10000; %时域抽样点数t1=linspace(0,T1-T1/N1,N1)'; %生成抽样时间点f1=1-2*t1; %生成抽样函数值OMG=32*pi; %频域范围K1=100; %频域抽样点数omg=linspace(-OMG/2,OMG/2-OMG/K1,K1)'; %生成抽样频率点X1=T1/N1*exp(-j*kron(omg,t1.'))*f1; %傅里叶正变换求解傅里叶系数fs1=OMG/2/pi/K1*exp(j*kron(t1,omg.'))*X1; %傅里叶逆变换还原时域函数T2=5; %五个周期时域范围N2=10000; %时域抽样点数t2=linspace(0,T2-T2/N2,N2)'; %生成抽样时间点fs2=0*t2;f2=sawtooth(t2*2*pi,0); %生成五个周期的锯齿波X2=T2/N2*exp(-j*kron(omg,t2.'))*f2; %傅里叶正变换求解傅里叶系数fs2=fs2+OMG/2/pi/K1*exp(j*kron(t2,omg.'))*X2; %傅里叶逆变换还原时域函数figure; %生成一个2*2矩阵子图subplot(2,2,1);plot(omg,abs(X1),'r'); %一个周期时的频谱图xlabel('Frequency'),ylabel('Amplitude')title('单个锯齿周期幅频特性曲线');subplot(2,2,2);plot(t1,fs1,'r'); %还原的时域函数xlabel('Time'),ylabel('Amplitude')title('Function after recovered');subplot(2,2,3);plot(omg,abs(X2),'r'); %五个周期时的频谱图xlabel('Frequency'),ylabel('Amplitude')title('五个锯齿周期幅频特性曲线');subplot(2,2,4);plot(t2,fs2,'r'); %还原的时域函数xlabel('Time'),ylabel('Function after recovered')title('Function after recovered');阅读两函数后不难发现，其实就是改变了linspace的范围。
matlab实验报告
matlab实验报告引言：Matlab（矩阵实验室）是一款功能强大的数值计算和科学计算软件，广泛应用于工程、科学和经济等领域。
本实验报告将探讨我在使用Matlab进行实验过程中的心得体会和实验结果。
实验一：图像处理在这个实验中，我使用Matlab对一张图像进行了处理，并应用了各种图像处理算法。
这包括图像增强、边缘检测和图像分割等技术。
通过Matlab的图像处理工具箱，我能够轻松调用各种算法函数，并对图像进行快速处理。
实验结果表明，Matlab图像处理工具箱提供了丰富的函数和算法，极大地方便了我们的图像处理工作。
实验二：模拟信号处理模拟信号处理是Matlab中的一个重要应用领域。
在这个实验中，我模拟了一个带噪声的正弦信号，并使用Matlab进行了噪声滤波和频谱分析。
通过使用Matlab的滤波函数，我能够有效地去除信号中的噪声，并还原出原始信号。
同时，Matlab提供了功能强大的频谱分析工具，我可以轻松地对信号的频率特性进行分析和可视化。
实验三：数据分析与统计数据分析与统计是Matlab的另一个重要应用领域。
在这个实验中，我使用Matlab对一组实验数据进行了分析和统计。
通过使用Matlab的统计函数和工具，我能够计算出数据的均值、方差、标准差等统计指标，并绘制出数据的直方图和散点图。
这些统计分析结果对我的实验研究提供了有力的支持，并帮助我更好地理解实验数据。
实验四：数值计算与优化数值计算与优化是Matlab的核心功能之一。
在这个实验中，我使用Matlab进行了一组数值计算和优化实验。
通过使用Matlab的数值计算函数和优化工具箱，我能够快速计算出复杂的数学问题，并找到最优解。
同时，在进行优化实验时，我可以设置各种约束条件和目标函数，从而得到最优解的参数值。
这些数值计算和优化工具极大地提高了我的研究效率和准确度。
结论：通过这些实验，我深刻认识到Matlab的强大功能和广泛应用领域。
无论是图像处理、信号处理、数据分析还是数值计算与优化，Matlab都提供了丰富的函数和工具，让我们能够快速高效地完成实验和研究工作。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验4、输入输出的应用
目的与要求：
掌握matlab 图形的绘制
掌握自定义函数的使用
掌握matlab 基础运用
要求：请将程序和运行结果的截屏附在题目后面。
一、编程题
1、产生一个1x10的随机矩阵，大小范围为[-5 5]，并且按照从大到小的顺序排列好！
2、设 y=cos[0.5+((3sinx)/(1+x^2))]把x=[0,2π]间分为101点，画出以x 为横坐标，y 为纵坐标的曲线；线条颜色为红色，线型为虚线。
3、在[-10，10；-10，10]范围内画出函数z 的三维网格图形。
4、将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。
①]2
2[)cos(π
π，，-∈x x （曲线图）； ②4)y 2,-4x (-242),(2222≤≤≤≤+=；y x y x f （曲面图）。
5、请修改下面的程序，让他们没有for 循环语句！修改后的程序与源程序得到相当的结果。
A=[1 2 3; 4 5 6; 7 8 9];
[r c]=size(A);
for i=1:1:r
for j=1:1:c
if (A(i,j)>8 | A(i,j)<2)
A(i,j)=0;
end
end
end
6、请自己定义函数，求出递归调用（即1*2*3*…*n）
（1）求前n项积，n为自然数。
（2）给出调用命令，如n=100时的结果。
二、填空题
1、已知a=[4 5 6],b=[7 8],[x,y]=meshgrid(a,b),则x=
y= 。
2、定义x=[0 1 2 3 4]的用冒号表达式表示为。
3、在算术、关系、逻辑三种运算符中，优先级最高。
4、设x=[4 8 1;9 7 6],[a,b]=size(x)得到的结果是a= b= 。
5、x=find(A)命令的含义是。
6、元素取整的命令有、。
7、a=[1 2 2]，求a的转置的matlab命令。
8、b=[1 4 5 2 3 7 8 11 9 13 15 6]，求所有元素和的matlab 命令。
9、syms t k的含义。
10、在线形图型格式的设置中，字符y表示什么颜色，,栅格的命令是。
11、一个字符串由多个字符组成，用符号来界定。