辽工大仿真技术综合设计_Matlab

课程设计报告题目：Matlab报告班级：姓名：学号：指导教师：成绩：电子与信息工程学院信息与通信工程系1求最大值和最小值1.1题目输入10个数，求其中最大数和最小数。

要求分别用循环结构和调用MATLAB 的max函数、min函数来实现。

1.2 用到的函数input：请求用户输入max：取数组中的最大值min：取数组中最小值1.2流程图图1-1求最大值最小值流程图1.4 运行截图图1-2调用max和min函数运行截图图1-2运用循环方法运行截图2求Fibonacci数列2.1 题目求Fibonacci数列。

(1)大于4000的最小项.(2)5000之内的项数2.2什么是Fibonacci数列即斐波那契数，亦称之为斐波那契数列（意大利语：Successione di Fibonacci），又称黄金分割数列、费波那西数列、费波拿契数、费氏数列，指的是这样一个数列：0、1、1、2、3、5、8、13、21、……在数学上，斐波纳契数列以如下被以递归的方法定义：F0=0，F1=1，F n=F n-1+F n-2（n>=2，n∈N*），用文字来说，就是斐波那契数列列由0 和1 开始，之后的斐波那契数列系数就由之前的两数相加。

(百度百科)2.2 用到的函数length 取得一个对象的长度2.3 流程图图2-1 2求Fibonacci数列流程图2.4 运行截图图2-2 Fibonacci运行截图3 用求逆解法和直接解法解方程组3.1 题目解方程组Ax＝b，分别用求逆解法与直接解法求其解3.2 用到的函数inv 矩阵求逆3.3 流程图图3-1用求逆解法和直接解法解方程组流程图3.4 运行截图图3-2用求逆解法和直接解法解方程组运行截图4 求N阶方阵A的行列式4.1 题目编一个m程序，求N阶方阵A的行列式的值。

4.2 用到的函数size(A):获取矩阵的行数和列数，当只有一个输出参数时，返回一个行向量，该行向量的第一个元素是数组的行数，第二个元素是数组的列数（多维，而length 是取最大的那一维）。

4.3流程图图4-1求N阶方阵A的行列式流程图4.4运行截图图4-2求N阶方阵A的行列式运行截图5计算xi=0.2*i各点上的P(xi)5.1 题目今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。

用到的函数poly2sym:把系数数组转换为符号多项式.5.2流程图图5-1 计算xi=0.2*i各点上的P(xi)流程图5.3运行截图图5-2 计算xi=0.2*i各点上的P(xi)运行截图6一维数组倒序6.1 题目试编一个m程序，将一维数组x中的N个数按颠倒的次序重新存储。

如N=5，原来x为：x=[ 1 3 5 7 9 ]而经过颠倒处理后x中数据的次序应该为：x=[ 9 7 5 3 1 ]用到的函数rot90:旋转数组90°（倒序）.6.2流程图图6-1 一维数组倒序流程图6.3 运行截图图6-2 一维数组倒序运行截图7绘制二维曲线7.1 题目在[0 2π]范围内绘制二维曲线图y=sin(x)*cos(5x)。

用到的函数plot: 二维曲线绘图。

7.2流程图图7-1绘制二维曲线流程图7.3运行截图图7-2绘制二维曲线运行截图8名字转二进制再以脉冲绘图输出5.1 题目把自己姓名中的各个拼音字母（小写）所对应的ASCII值转换为二进制序列，并用二进制脉冲绘图表示出来。

用到的函数double:用double强制转换后得到的是ASCII值.dec2bin:dec2bin把ASCII转换成一个字符串形式表示的二进制数.str2num:string转number.stem:stem画出来是离散函数，杆图。

plot是连续函数.流程图图8-1名字转二进制再以脉冲绘图输出流程图8.2运行截图图8-2名字转二进制再以脉冲绘图输出运行截图图5-2 5计算xi=0.2*i各点上的P(xi)运行截图9GUI与SIMULINK9.1 题目设计一个第七章提到的GUI界面、第八章提到的SIMULINK模型的设计。

9.2流程图图9-1 GUI与SIMULINK流程图9.3设计思路总因为题目要求是设计需要结合GUI界面以及Simulink，所以决定尝试运用GUI来控制Simulink中的参数，实现两者的结合，首先要在Simulink中设计了一个简单的仿真电路，接着进行GUI的设计，最后是编写脚本，通过get函数获取相应参数，再运用set_param更改相应Simulink中的参数，最后用sim函数运行，再通过axes和plot讲结果显示到GUI上即完成了一个GUI结合Simulink 的简单设计。

9.4 运行截图图9-2 GUI与SIMULINK运行截图9.5 不足及改进首先是界面上，由于水平有限以及时间的仓促，没有设计出比较科学的界面，界面不够美观，同时有个很大的问题，就是脚本将scope的前十个数显示在坐标上时，可能会因为不够十个数而发生数组越界的报错，如果是大于了10个数，没有什么问题，这是不行的，在日后还要改进，另外就是结束按钮的功能，因为水平有限，研究了半天，还是没整出来，还需更加努力。

10 拉普拉斯变换求信号f(t)=t*(e^-3)*u(t)的拉普拉斯变换.用到的函数laplace:进行拉普拉斯变换.10.2 流程图图10-1拉普拉斯变换流程图10.3运行截图图10-2拉普拉斯变换运行截图11 拉普拉斯逆变换求F(s)=(s^3+5s^2+9s+7)/(s^2+3s+2)的拉普拉斯反变换。

用到的函数ilaplace:进行拉普拉斯反变换.11.2流程图图11-1拉普拉斯反变换流程图11.3运行截图图11-2拉普拉斯反变换运行截图12 傅里叶变换12.1 题目求信号f(t)=u(t+1)-u(t-1)的傅里叶变换用到的函数fourier:傅里叶变换函数.12.2 流程图图12-1傅里叶变换流程图12.3运行截图图12-2傅里叶变换运行截图13 实验心得通过这次学习，我对matlab这个软件有了更加熟悉的认识，特别是学到了交互式的GUI编程，感觉很方便，也很容易上手。

对于matlab语言，感觉和其他语言比较起来，也是十分类似，在有了其他几门编程语言的基础上，基本上没有什么难度，并且matlab内置了许多函数，在解决某些问题方面，只需要一个函数就可以搞定了，很方便。

14 附件14.1 求最大值最小值代码for i=1:10a(i)=input('please input a number：'); endmaxa=max(a);mina=min(a);maxamina%运行后会在命令窗出现for i=1:10a(i)=input('请输入10个数字：'); endmaxa = a(1);mina = a(1);for i=2:10;if maxa<a(i)maxa=a(i);endif mina>a(i)mina=a(i);endendmaxamina14.2 质求Fibonccif(1) = 1;f(2) = 1;for i=3:30;%猜f(i) = f(i-1)+f(i-2);endfx=fy= length(x);for i=1:yif x(i)>4000&x(i)<5000b = x(i);disp('min='),disp(b(1))disp('i='),disp(i)%肯定只有一个了endend14.3用求逆解法和直接解法解方程组源代码A=[2 -1 0 0 0;-1 2 -1 0 0;0 -1 2 -1 0;0 0 -1 2 -1;0 0 0 -1 2];b=[1;0;0;0;0];x=inv(A)*bx=A\b14.4 位数求和源代码function s=mydet(A)n=size(A);if (n(1)~=n(2))%行不等于列，也就不是方阵return;endn=n(1);if (n==1)%只有一行的也不用算了s=A;return;ends=0;for i=1:ns=s-(-1)^i*A(1,i)*mydet(A(2:n,[1:i-1 i+1:n]));%1:i-1 表示一个行向量: 1,2,3,...,i-1%i+1:n 表示另一个行向量：i+1, i+2, ..., n%逗号用于按行连接上述两个行向量，得到1,2,3,...,i-1, i+1, i+2, ..., n%A(2:n,[1:i-1 i+1:n]）取这一列的所有元素%行列式的值：现将左下部分等价转换为零，然后对角线上的各数的积即为其值。

end14.5 计算xi=0.2*i各点上的P(xi)源代码p1=[1.0 0.0 0.0 -2.0 1.0];p2=[0.0 0.0 1.0 4.0 -0.5];p1x=poly2sym(p1);%把系数数组转换为符号多项式。

p2x=poly2sym(p2);%第二个参数没有就是默认用x。

disp(p1x);disp(p2x);p=p1x+p2xdisp(p);14.6 一位数组倒序源代码x=[ 1 3 5 7 9 ]xx=rot90(rot90(x))14.7 绘制二维曲线源代码x = 0:pi/10:2*pi; %单调递增的数列{0 pi/10 pi/10*2 ... 2*pi} pi就是圆周率y = sin(x).*cos(5*x);%y中的每一个值都是通过x中对应位置的值算出来的；plot(x,y)14.8 名字转二进制再以脉冲绘图输出源代码name='lugang';% b=dec2bin(int8(s))b=double(name)%用double强制转换后得到的是ASCII值c=dec2bin(b)%dec2bin把ASCII转换成一个字符串形式表示的二进制数c=c'%转置%转置格式不一样：A'中的A是[m x n]的矩阵；A.'中的A是[n x 1]的向量stem(str2num(c(:)))%string转number%stem画出来是离散函数，杆图。

plot是连续函数14.9 GUI与SIMULINKfunction varargout = LG(varargin)%LG M-file for LG.fig% LG, by itself, creates a new LG or raises the existing% singleton*.% H = LG returns the handle to a new LG or the handle to% the existing singleton*.% LG('Property','Value',...) creates a new LG using the% given property value pairs. Unrecognized properties are passed% varargin to LG_OpeningFcn. This calling syntax produces a% warning when there is an existing singleton*.%% LG('CALLBACK') and LG('CALLBACK',hObject,...) call the% local function named CALLBACK in LG.M with the given input% arguments.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help LG% Last Modified by GUIDE v2.5 25-Jun-2016 19:43:29% Begin initialization code - DO NOT EDITgui_Singleton = 1;%你第一次运行gui的时候，没有关闭%第二次再运行的时候，matlab不会新开一个gui,而是把第一次那个gui拿到前面。