南邮matlab软件设计实验报告
南邮MATLAB软件设计
通信与信息工程学院2015 / 2016 学年第 2 学期软件设计实验报告模块名称MATLAB软件设计专业通信工程学生班级B130102学生学号B13010229学生姓名郭智超指导教师孟田甜报告内容一、要求练习的实验部分1.题目：在时间区间 [0,10]中，绘出t e y t2cos 15.0--=曲线。
程序代码：t=0:0.02:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time');ylabel('amplitude'); title('figure 1'); grid实验结果：24681000.511.5timea m p l i t u d efigure 12.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序代码：t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点）y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*z; plot(t,y,':r');hold on; %保持当前图形 plot(t,z,'b-') xlabel('t'); ylabel('z=f(t)');title('逐段解析函数');legend('y=sin(t)','z=f(t)',4); %确定向量长度 axis([0 10 -1.2 1.2]); %限定坐标轴范围 hold off; grid实验结果：24681001tz =f (t )逐段解析函数3.题目： 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？（2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ?（4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ?程序：A=[1 2 3;4 5 6;7 8 9][U,S,V]=svd(A);tol=S(1,1)*3*epssum(diag(S)>tol) - rank(A) S(1,1) - norm(A)sqrt(sum(diag(S*S))) - norm(A,'fro') S(1,1)/S(3,3) - cond(A)S(1,1)*S(2,2)*S(3,3) - det(A) < tol V(:,1)'*null(A)abs(A*null(A)) < tol U(:,1:2) == orth(A)运行结果：（1） sum(diag(S)>tol) - rank(A) 的结果是0 （2） S(1,1) - norm(A) = 0的结果是0（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015 （4） S(1,1)/S(3,3) - cond(A) 的结果是 -8（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的结果是1 （6） V(:,1)'*null(A) 的结果是 0（7） abs(A*null(A)) < tol 的结果是1；1；1（8） U(:,1:2) = = orth(A) 的结果是 [1 1 ；1 1 ；1 1] 4.题目：求积分dt t x⎰0]sin[cos ，π20≤≤x 。
（最新版）MATLAB实验报告
（最新版）MATLAB实验报告实验一典型环节的MATLAB仿真一、实验目的1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型，进行仿真和调试。
1．运行MATLAB软件，在命令窗口栏“>>”提示符下键入simulink命令，按Enter 键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK仿真环境下。
2．选择File菜单下New下的Model命令，新建一个simulink 仿真环境常规模板。
3．在simulink仿真环境下，创建所需要的系统。
以图1-2所示的系统为例，说明基本设计步骤如下：1）进入线性系统模块库，构建传递函数。
点击simulink下的“Continuous”，再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2）改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标，即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数，数字之间用空格隔开；设置完成后，选择OK，即完成该模块的设置。
3）建立其它传递函数模块。
按照上述方法，在不同的simulink 的模块库中，建立系统所需的传递函数模块。
例：比例环节用“Math”右边窗口“Gain”的图标。
4）选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”，将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口，形成一个阶跃函数输入模块。
5）选择输出方式。
用鼠标点击simulink下的“Sinks”，就进入输出方式模块库，通常选用“Scope”的示波器图标，将其用左键拖至新建的“untitled”窗口。
程序设计实验报告(matlab)
程序设计实验报告(matlab)实验一: 程序设计基础实验目的：初步掌握机器人编程语言Matlab。
实验内容：运用Matlab进行简单的程序设计。
实验方法：基于Matlab环境下的简单程序设计。
实验结果：成功掌握简单的程序设计和Matlab基本编程语法。
实验二：多项式拟合与插值实验目的：学习多项式拟合和插值的方法，并能进行相关计算。
实验内容：在Matlab环境下进行多项式拟合和插值的计算。
实验方法：结合Matlab的插值工具箱，进行相关的计算。
实验结果：深入理解多项式拟合和插值的实现原理，成功掌握Matlab的插值工具箱。
实验三：最小二乘法实验目的：了解最小二乘法的基本原理和算法，并能够通过Matlab进行计算。
实验内容：利用Matlab进行最小二乘法计算。
实验方法：基于Matlab的线性代数计算库，进行最小二乘法的计算。
实验结果：成功掌握最小二乘法的计算方法，并了解其在实际应用中的作用。
实验六：常微分方程实验目的：了解ODE的基本概念和解法，并通过Matlab进行计算。
实验内容：利用Matlab求解ODE的一阶微分方程组、变系数ODE、高阶ODE等问题。
实验方法：基于Matlab的ODE工具箱，进行ODE求解。
实验结果：深入理解ODE的基本概念和解法，掌握多种ODE求解方法，熟练掌握Matlab的ODE求解工具箱的使用方法。
总结在Matlab环境下进行程序设计实验，使我对Matlab有了更深刻的认识和了解，也使我对计算机科学在实践中的应用有了更加深入的了解。
通过这些实验的学习，我能够灵活应用Matlab进行各种计算和数值分析，同时也能够深入理解相关的数学原理和算法。
这些知识和技能对我未来的学习和工作都将有着重要的帮助。
实验二MATLAB程序设计含实验报告
实验二MATLAB 程序设计一、 实验目的1.掌握利用if 语句实现选择结构的方法。
2.掌握利用switch 语句实现多分支选择结构的方法。
3.掌握利用for 语句实现循环结构的方法。
4.掌握利用while 语句实现循环结构的方法。
5.掌握MATLAB 函数的编写及调试方法。
二、 实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。
点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、 实验内容1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。
并输入几组典型值加以检验。
（提示：提示输入使用input 函数）2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。
其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。
要求：（1）用switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
（提示：注意单元矩阵的用法）3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。
重复此过程，最终得到的结果为1。
如：2?13?10?5?16?8?4?2?16?3?10?5?16?8?4?2?1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。
请为关键的Matlab 语句填写上相关注释，说明其含义或功能。
南邮MATLAB上机实验
实验一：MATLAB语言平台与SIMULINK工具箱一、目的与任务1．学习了解MATLAB语言环境；2．练习MATLAB命令；4．掌握SIMULINK工具箱的应用。
二、内容与要求1．实验内容（1）学习了解MATLAB语言环境；（2）练习MATLAB命令的基本操作；（3）掌握SIMULINK仿真环境的使用方法;(4) 掌握线性系统仿真常用基本模块的用法2．要求按照实验文档的要求与步骤完成实验，撰写实验报告。
三、实验基本知识：1.熟悉MATLAB环境：MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.了解SIMULINK模块库中各子模块基本功能四、上机练习1.学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）>> help eyeEYE Identity matrix.EYE(N) is the N-by-N identity matrix.EYE(M,N) or EYE([M,N]) is an M-by-N matrix with 1's onthe diagonal and zeros elsewhere.EYE(SIZE(A)) is the same size as A.EYE with no arguments is the scalar 1.EYE(M,N,CLASSNAME) or EYE([M,N],CLASSNAME) is an M-by-N matrix with 1'sof class CLASSNAME on the diagonal and zeros elsewhere. Example:x = eye(2,3,'int8');See also speye, ones, zeros, rand, randn.Reference page in Help browserdoc eye>> D=eye(3)D =1 0 00 1 00 0 12．任意创建3个变量，分别为数值、矩阵、字符串,保存到mydat.mat文件，再把my_data.mat文件中的变量读取到MatLab工作空间内；查找已创建变量的信息，删除变量；学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
matlab实验报告(实验1)
开课学院及实验室：机电学院 2012年11 月16日
（2）在E盘新建一个文件夹并把其加入搜索路径。
截屏相关过程结果并保存。
（3）逐步按照以下步骤进行并记录各个步骤的结果：
1）在指令窗中输入clear指令；
5）再次在指令窗中输入clear指令后，输入
8）把cow.m暂时移动到另一个无关的文件夹，如“我的文档”。
（2）保存全部变量为数据文件Mydata.mat；在对应文件夹中找到该文件，截屏后保存结果；
（4）删除全部内存变量；
Editor/Debugger和脚本编写初步
把以下指令编写为M脚本文件，并运行：
t=0:pi/50:4*pi;
y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid
什么？此方法也是查询函数注释的有效方法，请牢记。
南邮MATLAB软件设计
通信与信息工程学院2015 / 2016 学年第 2 学期软件设计实验报告模块名称MATLAB软件设计专业通信工程学生班级B130102学生学号B13010229学生姓名郭智超指导教师孟田甜报告内容一、要求练习的实验部分1.题目：在时间区间 [0,10]中，绘出t e y t2cos 15.0--=曲线。
程序代码：t=0:0.02:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time');ylabel('amplitude'); title('figure 1'); grid实验结果：24681000.511.5timea m p l i t u d efigure 12.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序代码：t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点）y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*z; plot(t,y,':r');hold on; %保持当前图形 plot(t,z,'b-') xlabel('t'); ylabel('z=f(t)');title('逐段解析函数');legend('y=sin(t)','z=f(t)',4); %确定向量长度 axis([0 10 -1.2 1.2]); %限定坐标轴范围 hold off; grid实验结果：24681001tz =f (t )逐段解析函数3.题目： 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？（2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ?（4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ?程序：A=[1 2 3;4 5 6;7 8 9][U,S,V]=svd(A);tol=S(1,1)*3*epssum(diag(S)>tol) - rank(A) S(1,1) - norm(A)sqrt(sum(diag(S*S))) - norm(A,'fro') S(1,1)/S(3,3) - cond(A)S(1,1)*S(2,2)*S(3,3) - det(A) < tol V(:,1)'*null(A)abs(A*null(A)) < tol U(:,1:2) == orth(A)运行结果：（1） sum(diag(S)>tol) - rank(A) 的结果是0 （2） S(1,1) - norm(A) = 0的结果是0（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015 （4） S(1,1)/S(3,3) - cond(A) 的结果是 -8（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的结果是1 （6） V(:,1)'*null(A) 的结果是 0（7） abs(A*null(A)) < tol 的结果是1；1；1（8） U(:,1:2) = = orth(A) 的结果是 [1 1 ；1 1 ；1 1] 4.题目：求积分dt t x⎰0]sin[cos ，π20≤≤x 。
matlab程序设计实验报告
matlab程序设计实验报告《MATLAB程序设计实验报告》摘要：本实验报告旨在介绍MATLAB程序设计的基本原理和实践操作，通过实验演示和分析，展示了MATLAB在工程领域的应用和重要性。
本报告详细介绍了MATLAB程序设计的基本语法和常用函数，以及如何利用MATLAB进行数据处理、图像处理、信号处理等工程应用。
通过本报告的学习，读者将能够掌握MATLAB程序设计的基本技能，为工程实践提供有力的支持。
1. 引言MATLAB是一种用于算法开发、数据可视化、数据分析和数值计算的高级技术计算语言和交互式环境。
它具有强大的数学计算功能和丰富的绘图工具，广泛应用于工程、科学和金融等领域。
本实验报告将介绍MATLAB程序设计的基本原理和实践操作，帮助读者快速掌握MATLAB的基本技能。
2. 实验目的本实验的主要目的是让读者了解MATLAB程序设计的基本语法和常用函数，掌握MATLAB在工程领域的应用和重要性。
通过实验演示和分析，展示MATLAB 在数据处理、图像处理、信号处理等方面的应用。
3. 实验内容（1）MATLAB程序设计的基本语法和常用函数（2）利用MATLAB进行数据处理的实验演示（3）利用MATLAB进行图像处理的实验演示（4）利用MATLAB进行信号处理的实验演示4. 实验步骤（1）学习MATLAB程序设计的基本语法和常用函数（2）编写MATLAB程序，实现数据处理、图像处理、信号处理等功能（3）进行实验演示和分析，展示MATLAB在工程领域的应用和重要性5. 实验结果与分析通过本实验的学习，读者将能够掌握MATLAB程序设计的基本技能，包括数据处理、图像处理、信号处理等方面的应用。
通过实验演示和分析，读者将了解MATLAB在工程领域的重要性，为工程实践提供有力的支持。
6. 结论MATLAB程序设计是一种强大的工程工具，具有广泛的应用前景。
通过本实验报告的学习，读者将能够掌握MATLAB程序设计的基本技能，为工程实践提供有力的支持。
MATLAB程序设计实验报告
MATLAB 程序设计实验报告一、实验目的1. 通过实验熟悉MATLAB 仿真软件的使用方法；2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征；3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征；4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。
二、实验设备1. 计算机：2. MATLAB R2007a 仿真软件三、实验原理对系统的时域分析信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。
（1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为)()(1t f t f =)(2t f +；)()(1t f t f =)(2t f *。
（2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。
连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。
（3）信号的平移、翻转和尺度变换信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。
信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。
)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。
南京邮电大学matlab软件设计.doc
南京邮电大学matlab软件设计南京邮电大学通信学院软件课程设计实验报告模块名称___MATLAB软件设计专业班级__通信工程_17班姓名____张步涛__ ____ 学号___ _10001018 实验日期2013年6 月1728日实验报告日期2013 年7 月 1 日一、要求练习的实验部分1．在时间区间[0,10]中，绘出曲线。
程序t00.110; y1-exp-0.5*t.*cos2*t; plott,y, r- ; shg 结果2. 写出生成如图E2-1所示波形的MATLAB脚本M文件。
图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。
程序tlinspace0,3*pi,500; ysint; asinpi/3; zy0.*y; zya.*aytol - rankA 是多少（2）S1,1 - normA 0 是多少（3）sqrtsumdiagS*S - normA, fro 的结果是什么（4）S1,1/S3,3 - condA 的结果是什么（5）S1,1*S2,2*S3,3 - detA tol - rankA，求a asumdiagStol - rankA disp 设bS1,1 - normA，求 b bS1,1 - normA disp 设csqrtsumdiagS*S - normA, fro ，求c csqrtsumdiagS*S - normA, fro disp 设dS1,1/S3,3 - condA，求 d dS1,1/S3,3 - condA disp 设eS1,1*S2,2*S3,3 - detA tol - rankA，求a a 0 设bS1,1 - normA，求b b 0 设csqrtsumdiagS*S - normA, fro ，求c c 3.5527e-015 设dS1,1/S3,3 - condA，求 d d -8 设eS1,1*S2,2*S3,3 - detA tol - rankA 的结果是0；（2）S1,1 - normA 0 的结果是0；（3）sqrtsumdiagS*S - normA, fro 的结果是3.5527e-015；（4）S1,1/S3,3 - condA 的结果是-8；（5）S1,1*S2,2*S3,3 - detA wavread Error using wavread Data compression format IMA ADPCM is not supported. 我在查阅了很多资料，在网上也查阅相关信息，花费了大量时间也没找出结果，最后发现在WAV格式的语音文件有两种格式，即PCM格式和IMA ADPCM格式，而在MATLAB中用wavread函数进行语音处理时，并不能直接处理IMA ADPCM格式的语音信号，经过格式转换之后（选择PCM格式），我们运行出了正确的结果。
南邮Matlab实验报告三参考模板
Matlab上机实验报告实验名称：MATLAB的数值运算班级：自动化二班学号： B11050216姓名：李鹏飞南京邮电大学2013年4月9日一、实验目的1．学习MATLAB 的基本矩阵运算； 2．学习MATLAB 的点运算；3．学习复杂运算。
二、实验基本知识1.基本矩阵运算；2.多项式运算；3.代数方程组；4.数值微积分。
三、上机练习1. 给a,b,c 赋如下数据:]6,46,23,4,2,6,3,8,0,1[,3568382412487,278744125431-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=c b a1)求a+b,a*b,a.*b,a/b,a./b,a^2,a.^2,length(c)的结果.2)求c 中所有元素的平均值、最大值.3)求d=b(2:3,[1,3]).程序代码：（1）>> a=[1,3,4;5,12,44;7,8,27]; >> b=[-7,8,4;12,24,38;68,-5,3]; >> c=[1,0,8,3,6,2,-4,23,46,6]; >> d=a+b d =-6 11 8 17 36 82 75 3 30>> e=a*be =301 60 1303101 108 6081883 113 413 >> f=a.*bf =-7 24 1660 288 1672476 -40 81 >> g=a/bg =0.0966 0.0945 0.0080-3.6125 1.5838 -0.5778-1.9917 0.9414 -0.2682>> h=a./bh =-0.1429 0.3750 1.00000.4167 0.5000 1.15790.1029 -1.6000 9.0000>> i=a^2i =44 71 244373 511 1736236 333 1109 >> j=a.^2j =1 9 1625 144 1936 49 64 729 >> k=length(c) k =10（2）>> l=max(c) l =46>> m=mean(c) m =9.1000（3）>> d=b(2:3,[1,3]) d =12 3868 32.求解方程⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001987654321X ，并对所得结果作出解释。
实验3 Matlab程序设计1实验报告
实验3 Matlab程序设计1实验报告
实验3 Matlab程序设计1自查报告。
在本次实验中，我学习了如何使用Matlab进行程序设计，并完
成了相应的实验任务。
在实验过程中，我遇到了一些困难，但通过
查阅资料和与同学讨论，最终顺利完成了实验。
首先，我学习了Matlab的基本语法和常用函数，包括变量的定义、数组的操作、条件语句和循环结构等。
在实验中，我成功地运
用了这些知识，编写了一些简单的程序来实现特定的功能。
其次，我学习了Matlab的绘图功能，包括如何绘制二维和三维
图形，以及如何对图形进行美化和标注。
通过实验，我掌握了
Matlab中绘图函数的使用方法，并成功地绘制了一些图形来展示实
验结果。
在实验过程中，我也遇到了一些问题。
例如，在编写程序时，
我经常会忘记Matlab的语法规则，导致程序出现错误。
此外，我在
绘图时也遇到了一些困难，比如不知道如何设置图形的颜色和线型。
但通过查阅Matlab的官方文档和向同学请教，我逐渐解决了这些问
题。
总的来说，本次实验让我对Matlab的程序设计和绘图功能有了更深入的了解，也提高了我的编程能力和解决问题的能力。
通过不断地练习和实践，我相信我会在Matlab编程方面取得更大的进步。
南邮MATLAB软件设计
通信与信息工程学院2015 / 2016 学年第 2 学期软件设计实验报告模块名称MATLAB软件设计专业通信工程学生班级B130102学生学号B13010229学生姓名郭智超指导教师孟田甜报告内容一、要求练习的实验部分1.题目：在时间区间 [0,10]中，绘出t e y t2cos 15.0--=曲线。
程序代码：t=0:0.02:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time');ylabel('amplitude'); title('figure 1'); grid实验结果：24681000.511.5timea m p l i t u d efigure 12.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序代码：t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点）y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*z; plot(t,y,':r');hold on; %保持当前图形 plot(t,z,'b-') xlabel('t'); ylabel('z=f(t)');title('逐段解析函数');legend('y=sin(t)','z=f(t)',4); %确定向量长度 axis([0 10 -1.2 1.2]); %限定坐标轴范围 hold off; grid实验结果：24681001tz =f (t )逐段解析函数3.题目： 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？（2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ?（4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ?程序：A=[1 2 3;4 5 6;7 8 9][U,S,V]=svd(A);tol=S(1,1)*3*epssum(diag(S)>tol) - rank(A) S(1,1) - norm(A)sqrt(sum(diag(S*S))) - norm(A,'fro') S(1,1)/S(3,3) - cond(A)S(1,1)*S(2,2)*S(3,3) - det(A) < tol V(:,1)'*null(A)abs(A*null(A)) < tol U(:,1:2) == orth(A)运行结果：（1） sum(diag(S)>tol) - rank(A) 的结果是0 （2） S(1,1) - norm(A) = 0的结果是0（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015 （4） S(1,1)/S(3,3) - cond(A) 的结果是 -8（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的结果是1 （6） V(:,1)'*null(A) 的结果是 0（7） abs(A*null(A)) < tol 的结果是1；1；1（8） U(:,1:2) = = orth(A) 的结果是 [1 1 ；1 1 ；1 1] 4.题目：求积分dt t x⎰0]sin[cos ，π20≤≤x 。
南邮软件课程设计MATLAB软件设计
南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计专业班级：__通信工程_姓名：______ ____学号：___ _09实验日期：2012 年 6 月 11—21日实验报告日期： 2012 年 6 月 26 日索引目录一、要求练习的实验部分 (1)1、练习题1 (1)2、练习题2 (1)3、练习题3 (2)4、练习题4 (5)5、练习题5 (6)6、练习题6 (7)7、练习题7 (8)8、练习题8 (11)9、练习题9 (13)10、练习题10 (14)11、练习题11 (17)12、练习题12 (20)二、模拟数字电路仿真实验 (21)(一)数字逻辑电路基础 (21)(二)组合逻辑电路仿真 (22)(三)时序逻辑电路仿真 (26)三、数字信号处理仿真实验 (28)（一）、利用Kaiser窗函数设计的FIR低通滤波器进行数字滤波 (28)（二）、DSB-SC信号的生成与解调 (35)四、通信系统实验 (40)模拟调制---DSB调制与解调 (40)一、要求练习的实验部分1．在时间区间 [0,10]中，绘出t e y t2cos 15.0--=曲线。
程序：（zoushuCode1.m ） t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg结果：(图图1.fig)2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序：（zoushuCode2.m ） t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*z; plot(t,y,':r'); hold on; plot(t,z,'b-')xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off 结果：（图图2.fig ）3. 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1）sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8）U(:,1:2) = = orth(A) 的运行结果是什么 ?程序：(zoushuCode3.m)clear; clc;disp('设 A=[1,2,3;4,5,6;7,8,9]，得')A=[1,2,3;4,5,6;7,8,9] %一个3*3矩阵disp('设 [U,S,V]=svd(A)，得')[U,S,V]=svd(A) %奇异值分解(SVD)USV:U m*m酉矩阵 V n*n酉矩阵 S 对角阵disp('设 tol=S(1,1)*3*eps，得')tol=S(1,1)*3*epsdisp('设 a=sum(diag(S)>tol) - rank(A)，求 a')a=sum(diag(S)>tol) - rank(A)disp('设 b=S(1,1) - norm(A)，求 b')b=S(1,1) - norm(A)disp('设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 c')c=sqrt(sum(diag(S*S))) - norm(A,'fro')disp('设 d=S(1,1)/S(3,3) - cond(A)，求 d')d=S(1,1)/S(3,3) - cond(A)disp('设 e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol ，求 e')e=S(1,1)*S(2,2)*S(3,3) - det(A) < toldisp('设 f=V(:,1) *null(A)，求 f')f=V(:,1)'*null(A)disp('设 g=abs(A*null(A)) < tol，求 g')g=abs(A*null(A)) < toldisp('设 h=U(:,1:2) == orth(A)，求 h')h=U(:,1:2) == orth(A)结果：（结果3.txt）设 A=[1,2,3;4,5,6;7,8,9]，得A =1 2 34 5 67 8 9设 [U,S,V]=svd(A)，得U =-0.2148 0.8872 0.4082-0.5206 0.2496 -0.8165-0.8263 -0.3879 0.4082S =16.8481 0 00 1.0684 00 0 0.0000V =-0.4797 -0.7767 -0.4082-0.5724 -0.0757 0.8165-0.6651 0.6253 -0.4082设 tol=S(1,1)*3*eps，得tol =1.1223e-014设 a=sum(diag(S)>tol) - rank(A)，求 aa =设 b=S(1,1) - norm(A)，求 bb =设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 cc =3.5527e-015设 d=S(1,1)/S(3,3) - cond(A)，求 dd =-8设 e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol ，求 ee =1设 f=V(:,1) *null(A)，求 f f =设 g=abs(A*null(A)) < tol ，求 g g =1 1 1设 h=U(:,1:2) == orth(A)，求 h h =1 1 1 1 1 1 结果分析：由上面的结果得出如下的题目答案：（1）sum(diag(S)>tol) - rank(A) 的结果是0； （2）S(1,1) - norm(A) = 0 的结果是0；（3）sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015； （4）S(1,1)/S(3,3) - cond(A) 的结果是-8；（5）S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是1； （6）V(:,1)'*null(A)的结果是0； （7）abs(A*null(A)) < tol 的结果是1 1 1 ；（8）U(:,1:2)==orth(A)的运行结果是1 1 1 1 1 1 4. 求积分dt t x⎰]sin[cos ，π20≤≤x 。
matlab实验一实验报告
matlab实验一实验报告实验一：Matlab实验报告引言：Matlab是一种强大的数学软件工具，广泛应用于科学计算、数据分析和工程设计等领域。
本实验旨在通过使用Matlab解决实际问题，探索其功能和应用。
一、实验目的本次实验的主要目的是熟悉Matlab的基本操作和常用函数，了解其在科学计算中的应用。
二、实验内容1. 数值计算在Matlab中，我们可以进行各种数值计算，包括基本的加减乘除运算，以及更复杂的矩阵运算和方程求解。
通过编写相应的代码，我们可以实现这些功能。
例如，我们可以使用Matlab计算两个矩阵的乘积，并输出结果。
代码如下：```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A * B;disp(C);```2. 数据可视化Matlab还提供了强大的数据可视化功能，可以将数据以图表的形式展示出来，更直观地观察数据的规律和趋势。
例如，我们可以使用Matlab绘制一个简单的折线图，来展示某个物体在不同时间下的位置变化。
代码如下：```matlabt = 0:0.1:10;x = sin(t);plot(t, x);xlabel('Time');ylabel('Position');title('Position vs. Time');```3. 图像处理Matlab还可以进行图像处理，包括图像的读取、处理和保存等操作。
我们可以通过Matlab对图像进行增强、滤波、分割等处理，以及进行图像的压缩和重建。
例如，我们可以使用Matlab读取一张图片，并对其进行灰度化处理。
代码如下：```matlabimg = imread('image.jpg');gray_img = rgb2gray(img);imshow(gray_img);```三、实验结果与分析在本次实验中，我们成功完成了数值计算、数据可视化和图像处理等任务。
matlab程序设计实验报告
matlab程序设计实验报告Matlab程序设计实验报告引言：Matlab（Matrix Laboratory）是一种强大的高级编程语言和环境，广泛应用于科学计算、数据分析和工程设计等领域。
本实验报告旨在介绍我在Matlab程序设计实验中的学习和实践经验。
一、Matlab基础知识1.1 Matlab的安装与配置在实验开始前，我们首先需要安装Matlab并进行相应的配置。
Matlab的安装过程相对简单，只需按照官方指引进行操作即可。
配置方面，我们可以设置工作目录、界面风格、字体大小等，以提高工作效率。
1.2 Matlab的基本语法Matlab的语法类似于其他编程语言，但也有一些特殊之处。
例如，Matlab中的变量名不区分大小写，函数名则区分大小写。
此外，Matlab还具有丰富的数学函数库，可以方便地进行各种数值计算。
二、Matlab程序设计实践2.1 数值计算Matlab以其强大的数值计算能力而闻名，我们可以使用Matlab进行各种数学运算和数值计算。
例如，我们可以使用Matlab求解线性方程组、计算矩阵的特征值和特征向量等。
2.2 图像处理Matlab提供了丰富的图像处理函数，可以对图像进行各种操作和处理。
例如，我们可以使用Matlab读取图像文件、调整图像的亮度和对比度、进行图像滤波等。
此外，Matlab还支持图像的显示和保存，方便我们进行结果的展示和分析。
2.3 数据可视化Matlab提供了强大的数据可视化功能，可以将数据以图表的形式直观地展示出来。
我们可以使用Matlab绘制各种类型的图表，如折线图、散点图、柱状图等。
此外，Matlab还支持对图表的样式、标签、标题等进行自定义，以满足不同的需求。
三、实验心得与体会通过这次Matlab程序设计实验，我深刻体会到了Matlab在科学计算和工程设计中的重要性。
Matlab不仅提供了丰富的数学函数库和工具箱，还具备直观的图形界面和友好的交互环境，使得我们能够快速、高效地进行各种计算和分析。
2016南邮matlab软件设计报告
软件设计报告姓名：。
学号：B130207212016年6月26日～2016年7月3日目录第一章软件设计介绍 (1)第二章软件开发平台简介 (2)第三章软件设计练习题 (3)3.1 (3)3.2 (3)3.3 (4)3.4 (9)3.5 (10)3.6 (11)3.7 (12)3.8 (13)第四章设计题第一题 (14)4.1设计电路图 (14)4.2设置电路参数 (15)4.3程序代码 (16)4.4运行结果与分析 (17)第五章设计题第二题 (20)5.1迭代法介绍 (20)5.2简单迭代法的函数M文件 (21)5.3简单迭代法的运行M文件 (21)5.4简单迭代法的结果 (22)5.5超松弛迭代法的函数M文件 (22)5.6超松弛迭代法的运行M文件 (23)5.7超松弛迭代法的结果 (23)第六章设计题第三题 (25)6.1简单迭代法的函数M文件 (25)6.2简单迭代法的运行M文件 (26)6.3简单迭代法的结果 (26)6.4超松弛迭代法的函数M文件 (26)6.5超松弛迭代法的运行M文件 (27)6.6超松弛迭代法的结果 (28)第七章设计题第四题 (29)7.1设计思想 (29)7.2程序代码 (30)7.3运行结果 (30)第八章软件设计总结 (32)第一章软件设计介绍一、软件设计的性质、目的本课程是一门技能训练课，为在校学生的必修课程。
通过设计培养学生良好的编程习惯，使学生掌握编程的基本步骤，规范学生的编写程序的格式，锻炼学生编写程序的能力，为将来编写较大程序的打好基础。
二、软件设计的任务通过集中时间上机实验，使学生充分发挥已有的编程能力，熟悉编译环境，选择若干课题，进行审题、分析、设计和调试，并根据要求写出实验报告。
三、实验内容软件设计课题共分A 、B两档。
A档为练习题，第一周完成；B档为选做题，第二周完成。
四、考核办法考核采取程序验收和面试的方法进行。
其中，程序验收及面试占50%，实验报告占50%。
南京邮电大学matlab软件设计(超详细).doc
南京邮电大学matlab软件设计(超详细).doc南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计专业班级：__通信工程姓名：____ ____学号：___实验日期：2013年 6 月17—28日实验报告日期： 2013 年 7 月 1 日一、要求练习的实验部分1．在时间区间[0,10]中，绘出t15.0-=曲线。
-ey t2cos程序：t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t);plot(t,y,'r-');shg结果：2.写出生成如图E2-1所示波形的MATLAB脚本M文件。
图中虚线为正弦波，要求它的负3处被削顶。
半波被置零，且在2程序：t=linspace(0,3*pi,500);y=sin(t);a=sin(pi/3);z=(y>=0).*y;z=(y>=a).*a+(y<a).*z;plot(t,y,':r');hold on;plot(t,z,'b-')xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off结果：3. 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行[U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题：（1） sum(diag(S)>tol) - rank(A) 是多少 ？（2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro')的结果是什么 ?（4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5）S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？（6）V(:,1)'*null(A) 得到什么结果? （7）abs(A*null(A)) < tol 得到什么结果?（8）U(:,1:2) = = orth(A) 的运行结果是什么?程序：clear;clc;disp('设 A=[1,2,3;4,5,6;7,8,9]，得')A=[1,2,3;4,5,6;7,8,9] %一个3*3矩阵disp('设 [U,S,V]=svd(A)，得')[U,S,V]=svd(A) %奇异值分解(SVD)USV:U m*m酉矩阵 V n*n酉矩阵 S 对角阵disp('设 tol=S(1,1)*3*eps，得')tol=S(1,1)*3*epsdisp('设 a=sum(diag(S)>tol) - rank(A)，求a')a=sum(diag(S)>tol) - rank(A)disp('设 b=S(1,1) - norm(A)，求 b')b=S(1,1) - norm(A)disp('设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 c')c=sqrt(sum(diag(S*S))) - norm(A,'fro') disp('设 d=S(1,1)/S(3,3) - cond(A)，求 d') d=S(1,1)/S(3,3) - cond(A)disp('设 e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol ，求 e')e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol disp('设 f=V(:,1) *null(A)，求 f')f=V(:,1)'*null(A)disp('设 g=abs(A*null(A)) < tol，求 g')g=abs(A*null(A)) < toldisp('设 h=U(:,1:2) == orth(A)，求 h')h=U(:,1:2) == orth(A)结果：设 A=[1,2,3;4,5,6;7,8,9]，得A =1 2 34 5 67 8 9设 [U,S,V]=svd(A)，得U =-0.2148 0.8872 0.4082 -0.5206 0.2496 -0.8165 -0.8263 -0.3879 0.4082S =16.8481 0 0 0 1.0684 0 0 0 0.0000V =-0.4797 -0.7767 -0.4082-0.5724 -0.0757 0.8165-0.6651 0.6253 -0.4082设 tol=S(1,1)*3*eps，得tol =1.1223e-014设 a=sum(diag(S)>tol) - rank(A)，求 aa =设 b=S(1,1) - norm(A)，求 bb =设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 cc =3.5527e-015设 d=S(1,1)/S(3,3) - cond(A)，求 dd =-8设 e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol ，求 ee =1设 f=V(:,1) *null(A)，求 ff =设 g=abs(A*null(A)) < tol，求 gg =111设 h=U(:,1:2) == orth(A)，求 hh =1 11 11 1结果分析：由上面的结果得出如下的题目答案：（1）sum(diag(S)>tol) - rank(A) 的结果是0；（2）S(1,1) - norm(A) = 0 的结果是0；（3）sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015；（4）S(1,1)/S(3,3) - cond(A) 的结果是-8；（5）S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是1；（6）V(:,1)'*null(A)的结果是0；（7）abs(A*null(A)) < tol的结果是111 ；（8）U(:,1:2)==orth(A)的运行结果是1 11 11 1 4.求积分dt tx⎰0]sin[cos，π20≤≤x。
《MATLAB软件设计》实验报告
《MATLAB 软件设计》实验报告学院：&&&&&&&学院 学号：！！！！！ 姓名：****1、写代码实现或填空（5分）：（1）将矩阵123456789876⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦赋给变量a ；a=[1,2,3,4;,5,6,7,8;9,8,7,6;]（2）将a 的第二行赋给变量b ；b=a(2,:)（3）将a 的第三列赋给变量c ；c=a(:,3)（4）a(2,3)=([7])；（5）a(1:2,2:3)=(⎥⎦⎤⎢⎣⎡7362)；2、（35分）写代码段示意：从文本读入数据，用所读入数据绘制二维曲线图，二维伪彩色图，给所做图形添加标注，图例，色标，文本等，并用所学句柄绘图内容对图形进行精细修改。
x=load('as.txt');z=rand(12,7);y=[1:12]pcolor(x,y,z);Shading flat;shading interp;xlabel('x');ylabel('y');legend('sth');h6=colorbar;title('伪色彩图');h1=xlabel('x');set(h1,'string','x 轴','color',[1 0 1],'fontsize',11);h2=ylabel('y');set(h2,'string','y 轴','color',[1 0 1],'fontsize',11);set(gcf,'color',[0.4 0.3 0.7]);set(gca,'box','off');set(gca,'xlim',[0 50],'ylim',[0 10]);set(gca,'yscale','log','xscale','linear');set(gca,'xtick ',0:20:50);H=uicontrol(gcf,'Style','pushbutton','Position',[10 5 100 25],'string','关闭','CallBack','close','min',-90,'max',90,'value',10)3、（5分）结合代码段x=(1:10)/10*2*pi; y=sin(x); h1=plot(x,y);xlabel('x(km)');的执行结果填空或回答问题：（1） 获得h1的父对象句柄的指令为（ f=get(h,’parent ’)），此父对象为（ f =151.0012 ）；（2）利用句柄h1修改曲线的线型为虚线、宽度为20、颜色为红色；set(h1,'color',[1 0 0]);set(h1,'LineWidth',[20]);set(h1,'LineStyle','--');（3）获得X 轴标注的句柄h2；h2=xlabel('x(km)')（4）利用句柄h2修改X 轴的标注为distance(km)，字体颜色为红色；set(h2,'String','distance(km)');set(h2,'color',[1 0 0]);4、（5分）排列下列对象的层次关系：menu, pushbutton, axes, line, figure, root menuUicontrol ——pushbuttonRoot ——figure -Axes ——line5、（15分）编写代码段示意分支结构和循环结构的用法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
通信与信息工程学院2014 /2015 学年第2学期软件设计实验报告模块名称MATLAB软件设计专业通信工程学生班级通工2班学生学号学生姓名指导教师孟田甜报告内容一、要求练习的实验部分1． 题目：在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。
实现方法：plot 画图 程序：clccleart=0:0.01:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); xlabel('t'); ylabel('y');title('绘制曲线'); 结果及结果分析：1234567891000.511.5ty绘制曲线对源程序的改进：1.加了坐标轴标示，标题2.将步长细化2． 题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
图 E2-1 题15所要表示的波形（逐段解析函数）实现方法：plot 画图 程序 t=0:0.01*pi:10;y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*z; plot(t,y,':r'); hold on ;plot(t,z,'b-') xlabel('t');ylabel('z=f(t)');title('逐段解析函数');legend('y=sin(t)','z=f(t)',4); axis([0 10 -1.2 1.2]); hold off ; 结果及结果分析12345678910-1-0.8-0.6-0.4-0.200.20.40.60.81tz =f (t )逐段解析函数对源程序的改进：1.使用了常规的t 生成方法。
2.因为原程序图形上下都顶到了，不好看。
所以加上规定坐标轴值范围的语句3． 题目： 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题：（1） sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ?（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ?（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ?实现方法：matlab 矩阵操作 程序： clcclearA=[1 2 3;4 5 6;7 8 9]; [U,S,V]=svd(A);tol=S(1,1)*3*eps;sum(diag(S)>tol)-rank(A) S(1,1)-norm(A)sqrt(sum(diag(S*S)))-norm(A,'fro') S(1,1)/S(3,3)-cond(A)S(1,1)*S(2,2)*S(3,3)-det(A)<tol V(:,1)'*null(A) abs(A*null(A))<tol U(:,1:2)==orth(A)结果及结果分析：（1） sum(diag(S)>tol) - rank(A) 的结果是0 （2） S(1,1) - norm(A) = 0的结果是0（3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是0 （4） S(1,1)/S(3,3) - cond(A) 的结果是0（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的结果是1 （6） V(:,1)'*null(A) 的结果是-1.3124e-016 （7） abs(A*null(A)) < tol 的结果是1；1；1（8） U(:,1:2) = = orth(A) 的结果是1，1 ；1，1 ；1，14． 题目： 求积分dt t x⎰]sin[cos ，π20≤≤x 。
提示：abs,sin,cos,cumsum,disp,plot,grid实现方法：使用int 函数程序：F1=vpa(int(abs(sin(cos(t))),t,0,2*pi))结果及结果分析：F1 =3.57297496390010467337484720250395． 题目： 求方程2,122==+xy y x 的解。
实现方法:使用solve 函数 程序：syms xsyms yA=solve('x*x+y*y=1','x*y=2',x,y); Ax=A.x Ay=A.y结果及结果分析： Ax =-1/2*(1/2*5^(1/2)+1/2*i*3^(1/2))^3+1/4*5^(1/2)+1/4*i*3^(1/2) -1/2*(1/2*5^(1/2)-1/2*i*3^(1/2))^3+1/4*5^(1/2)-1/4*i*3^(1/2) -1/2*(-1/2*5^(1/2)+1/2*i*3^(1/2))^3-1/4*5^(1/2)+1/4*i*3^(1/2) -1/2*(-1/2*5^(1/2)-1/2*i*3^(1/2))^3-1/4*5^(1/2)-1/4*i*3^(1/2) Ay =1/2*5^(1/2)+1/2*i*3^(1/2) 1/2*5^(1/2)-1/2*i*3^(1/2) -1/2*5^(1/2)+1/2*i*3^(1/2) -1/2*5^(1/2)-1/2*i*3^(1/2)6． 题目： 在某激励条件下，二阶系统归一化响应可表示为)sin(11)(θββξ+-=-t e t y t ，其中ξ为阻尼系数，21ξβ-=，)1(2ξξθ-=arctg 。
请用不同的颜色或线型，在同一张图上，绘制0.2:2.0:2.0=ξ取值下系统在]18,0[∈t 区间内的响应曲线，并要求用2.0=ξ和0.2=ξ对它们相应的两条曲线进行醒目的文字标志。
实现方法: 使用了 text,for,if 函数 程序 clc clear syms k ; t=0:0.1:18; hold on ;%kz={'-r' '.g' 'b' 'c' 'm' 'y' 'k' 'w' ':r' ':g'}; kz='r g :b c m y k w r g r:'; i=0;for k=0.2:0.2:2 i=i+1; clear y ; a=sqrt(1-k^2); b=atan(sqrt(1-k^2)/k); if k==1y=(exp(t)-t-1)./exp(t); plot(t,y,kz(i)) elsey=1-1/a*exp(-k*t).*sin(a*t+b); plot(t,y,kz(i)) end% plot(t,y,kz(i)) if k==0.2text(8,1.2,'阻尼系数是0.2')end if k==2text(14,0.9,'阻尼系数是2') end endtitle('二阶系统归一化响应'); xlabel('t'); ylabel('y(t)'); hold off程序2（求极限）：syms klimit('1-1/sqrt(1-k^2)*exp(-k*t)*sin(sqrt(1-k^2)*t+atan(sqrt(1-k^2)/k))',k,1)结果及结果分析：24681012141618-0.200.20.40.60.811.21.41.6二阶系统归一化响应ty (t )7． 题目： 构建一个简单的全波整流模型，并用示波器分两路同时观察原信号和整流后的信号波形。
要求：信源模块被重新命名为“输入正弦波”；信宿模块被重新命名为“示波器”；连接到信宿的信号线上分别标注“原信号”和“整流信号”；在模型窗中添加注释文本。
实现方法: 用sine wave 产生正弦波，用abs 进行整流，用Scope 进行信号的观察。
程序:结果及结果分析结果分析：较好的实现了全波整流。
8． 题目： 利用 SIMULINK 及其标准模块设计一个低通滤波器，从受噪声干扰的多频率混合信号)(t x 中获取10Hz的信号。
在此)()1002cos()102sin()(t n t t t x +⋅⋅+⋅⋅=ππ，而)2.0,0(~)(N t n实现方法: 用Sine Wave 产生正弦波，用Random Number 产生随机信号，用sum 将信号相加，用Analog Filter Design 设计低通滤波器，用Mux 模块将信号合并在一起,用Scope 观察信号； 程序:结果及结果分析结果分析：LPC 使用Butterworth 方法，实现了不错的低通滤波效果。
9． 题目：已知系统的状态方程为：⎩⎨⎧=--=1222211)1(y y y y y y ，其中25.0)0(,25.0)0(21==y y ，请构建该系统的仿真模型，并用XY Graph 模块观察21,y y 相轨迹。
实现方法：利用乘法器，加法器，积分器搭建模型，利用XY Graph 模块观察21,y y 相轨程序：结果及结果分析：10．题目：CD74HC00仿真（数字电路）实现方法: 这里就是四个与非门的实现，利用Logical Operator模块做出四个与非门；用两个Pulse Generator分别产生“1100”和“0101”的序列码，检验系统对于“00”、“01”、“10”、“11”四种输入的输出。
程序:结果：结果分析:前两路分别为产生“1100”和“0101”的脉冲发生器的波形，后四路波形则分别为四个与非门的输出；可以看到，四路输出信号与输入信号均符合与非逻辑的关系。
11．题目： 2FSK信号的生成与接收滤波要求：1)产生二进制[0，1]随机信号波，信号波特率为1000B2)对此信号进行2FSK调制，f1=1500Hz, f2=3050Hz3)加入高斯白噪声4)在接收端进行带通滤波5)抽样速率为20KHz.6)观察滤波前后信号的波形和频谱。