当前位置：文档之家 › Matlab 实验报告

Matlab 实验报告

  • 格式：docx
  • 大小：358.05 KB
  • 文档页数：22

四、程序截图
1）
2）
3）
3）
5）
五、实验结果截图
1）
2）
3）
4）
[1]r=0.2
时
[2]r=0.5
r=0.8
时
[3]r=1.2
时
5)椭圆面积 p=4.8442
六、实验结语
通过这次实验，使我掌握了 matlab 这个强大的数学工具的 基本操作，也使我认识到，科学的进步可以使我们以更加准确的 工具来模拟难以实际实验的一些数学实验现象。 动态的变化过程， 形象地刻画了事物发展的渐变规律，混沌现象更使我体会到“量 变达到质变”的数学表现形式。 最后，谢谢老师的无私教导，指导我成功完成实验！
Matlab 实验报告
班级软件 12 学号 2111601027 姓名李旭锋 日期 2012有 5 个小实验，其要求如下： 1） 在 同 一 坐 标 系 下 画 出 函 数 y=sinx, y=cosx,
y=0.2e^(0.1*x)sin(0.5x) 和 y=0.2e^(0.1*x)cos(0.5x) 在 区 间 [0,2pi]的曲线图，并对该图进行修饰； 2） 每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进 行统计,统计内容包括:全班人数,总得分,平均得分,不及格的 人数及 90 分(包括 90 分)以上的人数.请编制程序解决这一 问题,并自给一组数据验证程序的正确性； 3） 编制程序：要求任给一个自然数 n，写出 1 到 n 之间所有 能够表示成连续若干个自然数的和的数及其所有表示方法. （如 6=1+2+3; 15=1+2+3+4+5 或 15=4+5+6 或 15=7+8 等等） ； 4） Feigenbaum 曾对超越函数 y=rsin（������x） （r 为非负实数）进 行了分岔与混沌的研究， 试利用迭代格式 x （k+1） =rsin （������x （k ） ） ，做出相应的 Feigenbaum 图； 5） 计算椭圆
x=0.1; y=[]; r=3.2; hold on axis([0 100 0 1]) for i=1:100 x=r*sin(pi*x); y=[y,x]; plot(i,x,'k.','markersize',10) fprintf('x(%d)=%.10f\n',i,x);pause(0.005); end t=1:100; plot(t,y,'k-'); grid 5）程序 f=inline('sqrt(3.*sin(t).*sin(t)+1)'); t=(0:0.1:1)*pi; y=f(t); p=trapz(t,y); fprintf('p=%.5f\n',p)
2）程序 clear;clc; count1=0; count2=0; n=input('please inpue n='); a=input('please input the first mark:') if(a<60) count1=count1+1; else if(a>=90) count2=count2+1; end end t=a; for i=1:n-1 p=input('please input the next mark:'); t=t+p; if(p<60) count1=count1+1; else if(p>=90) count2=count2+1; end end
2）用 count1 表示不及格的人数，用 count2 表示 90 分以上（包 括 90 分）的人数，n 表示总人数，ave 则表示平均成绩。每输入 一个成绩，用 if 判断其所属范围，如小于 60 则 count1 自增 1， 如大于等于 90 则 count2 自增 1，与此同时，用 t 累计 n 个学生 的总成绩，然后 ave 的值则用 t 与 n 的比值确定。结果输出不及 格人数 count1 和不低于 90 分的人数 count2、 以及平均成绩 ave； 3）用 input 命令输入欲求解的数，用 floor 函数降低算法的复杂 度，以加快计算的速度。采用 for 循环实现对此数的和式因子求 解，其因子的和可用式 k*n+k*(k-1)/2 来表示，用 if 判断若因子 和等于所求数，则输出因子的连加形式，并且控制输出格式为每 行输出一种因子表示； 4）按比例系数的情况分三种情况试验： [1]取比例系数 r 为 0.1 到 0.3 之间的任意一个值， 按迭代 式 x（k+1）=rsin（������x（k） ）迭代 100 步，产生迭代序列{x（k）}， =0,1，„„，100 并绘图； [2]再取比例系数为 0.4 到 0.8 之间的任意一个数值， 并运 行上述程序； [3]再取比例系数 r 为 0.9 到 3.2 之间的任意一个值运行上 述程序进行试验，其实呈现一种紊乱状态； 5） 首先进行数学建模。 由高等数学知识可知， 用代换 x=cox （t） ， y=sin(t),0<=t<= ������ ， 则 所 求 问 题 刻 画 为 求 一 定 积 分 f=2
������ 0
3 sin x ∗ sin⁡ (x) + 1的值。用 inline 建立内联函数，然后
用命令 trapz 进行积分。
三、程序设计
其程序列单如下： 1）程序 x=(0:0.1:2)*pi; y=sin(x); z=cos(x); s=0.2.*exp(0.1.*x).*sin(0.5.*x); t=0.2.*exp(0.1.*x).*cos(0.5.*x); plot(x,y,'b-',x,z,'g-','linewidth',3,'markersize',15) text(2.5,0.7,'sin(x)') text(1.6,-0.6,'cos(x)') hold on plot(x,s,'m-',x,t,'r-','linewidth',3,'markersize',15) text(4.0,0.3,'s(x)') text(4.0,-0.1,'t(x)') xlabel('Variable \it{x}') ylabel('Variable \it{y}') title('four figures are here,thank you') plot([0,2*pi],[0,0],'r-.') grid hold off
������ ^2 4
+y^2=1 的周长，使结果保留五位有效数字。
二、问题的分析
1） 使用命令 hole on 实现在同一窗口中多次绘制图形， 用 hold off 取消；使用 grid 命令对图形性窗口加坐标网格；使用 linewidth 和 markersize 命令实现图形中线宽和点型大小的设置；使用 title,xlable,ylable,text 等命令实现对图形的文字说明；
ave=t/n; end fprintf('n=%d count1=%d count2=%d
ave=%.4f',n,count1,count2,ave); 3）程序 Clear;clc; f=input('please input your number=') n=0; for i=1:f-1 n=n+1; for k=1:f/n if f==k*n+k*(k-1)/2 fprintf('num=%.0f',n) for i=1:k-1 fprintf('+%.0f',n+i) end fprintf ('\n') end end end 4）程序 clc;clf;

合集下载

Matlab实验报告(三)-MATLAB绘图

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; ％关闭所有的图形视窗。

x=1:10；y=rand （size(x ）); bar(x,y ）; ％绘制直方图.123456789100.51Bar(）函数还有barh （）和errorbar （）两种形式，barh()用来绘制水平方向的直方图，其参数与bar()相同，当知道资料的误差值时，可用errorbar （)绘制出误差范围,其一般语法形式为：errorbar （x,y,l,u)其中x,y 是其绘制曲线的坐标,l ，u 是曲线误差的最小值和最大值，制图时，l 向量在曲线下方,u 向量在曲线上方。

基于matlab的实验报告

基于matlab的实验报告

基于matlab的实验报告实验报告：基于MATLAB 的实验一、实验目的通过使用MATLAB 软件，掌握如何进行数据分析、图像处理、算法实现等一系列实验操作，提高实验者的实践能力和动手能力。

二、实验原理MATLAB 是一种在科学计算和技术开发领域广泛应用的计算机软件。

它能进行矩阵计算、绘制函数和数据图像、实现算法以及进行数据分析等。

通过掌握MATLAB 的使用，能够快速、高效地解决各种科学和工程问题。

三、实验内容1. 数据分析：使用MATLAB 的数据分析工具进行数据的导入、处理和分析。

2. 图像处理：利用MATLAB 的图像处理工具包对图像进行滤波、增强、分割等操作。

3. 算法实现：使用MATLAB 实现常用的算法，如排序、搜索、图像压缩等。

四、实验步骤1. 数据分析：（1）使用MATLAB 的读取数据函数将数据导入MATLAB 环境中。

（2）利用MATLAB 的数据处理函数进行数据清洗和预处理。

（3）使用MATLAB 的统计工具进行数据分析，如求平均值、标准差等。

（4）利用MATLAB 的绘图函数将分析结果可视化。

2. 图像处理：（1）使用MATLAB 的读取图像函数将图像导入MATLAB 环境中。

（2）利用MATLAB 的图像处理工具包进行滤波操作，如均值滤波、中值滤波等。

（3）使用MATLAB 的图像增强函数对图像进行锐化、变换等操作。

（4）利用MATLAB 的图像分割算法对图像进行分割。

3. 算法实现：（1）使用MATLAB 编写排序算法，如冒泡排序、快速排序等。

（2）使用MATLAB 编写搜索算法，如二分查找、线性搜索等。

（3）使用MATLAB 实现图像压缩算法，如离散余弦变换（DCT）。

五、实验结果实验中，我们使用MATLAB 完成了数据分析、图像处理和算法实现的一系列实验操作。

通过数据分析，我们成功导入了数据并对其进行了清洗和预处理，最后得到了数据的统计结果。

在图像处理方面，我们对图像进行了滤波、增强和分割等操作，最终得到了处理后的图像。

MATLAB实验报告(1-4)

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信号与系统实验报告19472[五篇范文]

MATLAB信号与系统实验报告19472[五篇范文]

MATLAB信号与系统实验报告19472[五篇范文]第一篇：MATLAB信号与系统实验报告19472信号与系统实验陈诉（5）MATLAB 综合实验项目二连续系统的频域阐发目的：周期信号输入连续系统的响应可用傅里叶级数阐发。

由于盘算历程啰嗦，最适适用MATLAB 盘算。

通过编程实现对输入信号、输出信号的频谱和时域响应的盘算，认识盘算机在系统阐发中的作用。

任务：线性连续系统的系统函数为11)(+=ωωjj H，输入信号为周期矩形波如图 1 所示，用MATLAB 阐发系统的输入频谱、输出频谱以及系统的时域响应。

-3-2-1 0 1 2 300.511.52Time(sec)图 1要领：1、确定周期信号 f(t)的频谱nF&。

基波频率Ω。

2、确定系统函数 )(Ω jn H。

3、盘算输出信号的频谱n nF jn H Y&&)(Ω=4、系统的时域响应∑∞-∞=Ω=nt jnn eY t y&)（MATLAB 盘算为y=Y_n*exp(j*w0*n“*t);要求（画出 3 幅图）：1、在一幅图中画输入信号f(t)和输入信号幅度频谱|F(jω)|。

用两个子图画出。

2、画出系统函数的幅度频谱|H(jω)|。

3、在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(jω)|。

用两个子图画出。

解:(1)阐发盘算：输入信号的频谱为(n)输入信号最小周期为=2，脉冲宽度，基波频率Ω=2π/ =π，所以(n)系统函数为因此输出信号的频谱为系统响应为(2)步伐：t=linspace(-3,3,300);tau_T=1/4;%n0=-20;n1=20;n=n0:n1;%盘算谐波次数20F_n=tau_T*Sa(tau_T*pi*n);f=2*(rectpuls(t+1.75,0.5)+rectpuls(t-0.25,0.5)+rectpuls(t-2.25,0.5));figure(1),subplot(2,1,1),line(t,f,”linewidth“,2);%输入信号的波形 axis([-3,3,-0.1,2.1]);grid onxlabel(”Time(sec)“,”fontsize“,8),title(”输入信号“,”fontweight“,”bold“)%设定字体巨细，文本字符的粗细text(-0.4,0.8,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“);%输入信号的幅度频谱xlabel(”n“,”fontsize“,8),title(”输入信号的幅度频谱“,”fontweight“,”bold“)text(-4.0,0.2,”|Fn|“)H_n=1./(i*n*pi+1);figure(2),stem(n,abs(H_n),”.“);%系统函数的幅度频谱xlabel(”n“,”fontsize“,8),title(”系统函数的幅度频谱“,”fontweight“,”bold“)text(-2.5,0.5,”|Hn|“)Y_n=H_n.*F_n;y=Y_n*exp(i*pi*n”*t);figure(3),subplot(2,1,1),line(t,y,“linewidth”,2);%输出信号的波形 axis([-3,3,0,0.5]);grid onxlabel(“Time(sec)”,“fontsize”,8),title(“输出信号”,“fontweight”,“bold”)text(-0.4,0.3,“y(t)”)subplot(2,1,2),stem(n,abs(Y_n),“.”);%输出信号的幅度频谱xlabel(“n”,“fontsize”,8),title(“输出信号的幅度频谱”,“fontweight”,“bold”)text(-4.0,0.2,“|Yn|”)(3)波形：-3-2-1 0 1 2 300.511.52Time(sec)输入信号f(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输入信号的幅度频谱|Fn|-20-15-10-5 0 5 10 15 2000.10.20.30.40.50.60.70.80.91n系统函数的幅度频谱|Hn|-3-2-1 0 1 2 300.10.20.30.4Time(sec)输出信号y(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输出信号的幅度频谱|Yn| 项目三连续系统的复频域阐发目的：周期信号输入连续系统的响应也可用拉氏变更阐发。