南京邮电大学matlab软件设计(超详细).doc
南邮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 。
南京邮电大学软件设计实验报告
软件设计报告( 2014 / 2015 学年第二学期)课程名称软件设计指导老师赵江实习时间第十八周学生姓名学号____学院______专业软件设计课程编号：B0465011C适用专业：班级：一、所涉及的课程及知识点涉及的课程：第6学期之前的专业基础课程。
知识点：专业基础课程中所学的知识点。
二、目的与任务目的：通过软件设计，培养学生的实践能力和创新精神，加强学生对专业基础课程的理解和掌握，加强学生高级语言编程能力、应用软件以及仿真能力。
任务：选择以下任一模块进行设计：Matlab软件仿真、C语言及应用。
软件设计的内容题目1：如果给出两个矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=136782078451220124A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=087654321B ，执行下面的矩阵运算命令。
（1）B A *5+和I B A +-分别是多少（其中I 为单位矩阵） （2）B A *⋅和B A *将分别给出什么结果，它们是否相同为什么 逻辑功能程序：function [ ] = EXP1()A=[4,12,20;12,45,78;20,78,136]; B=[1,2,3;4,5,6;7,8,0]; I=eye(3); disp('A+5*B='); disp(A+5*B); disp('A-B+I=') disp(A-B+I); disp('A.*B='); disp(A.*B) disp('A*B='); disp(A*B); End 实验过程与结果打开matlab ，在命令窗口“Command Window ”中键入edit,启动程序编辑器。
输入完整程序后利用save as 储存为M 文件，文件名为EXP1。
返回主界面，在命令窗口“Command Window”中输入函数EXP1()，按下回车，得到程序运行结果如下：>> EXP1( )A+5*B=9 22 3532 70 10855 118 136A-B+I=4 10 178 41 7213 70 137A.*B=4 24 6048 225 468140 624 0A*B=192 228 84738 873 3061284 1518 528实验结果分析（1）利用MATLAB提供的disp函数既可以输出表达式、数值，也可以输出字符串，其调用方式为：disp(表达式或数值）、disp（‘待显示字符串’）；（2）在MATLAB的矩阵运算中，+、-运算符通用，表示矩阵相加、减；*与.*不同在于*表示矩阵乘法，而.*表示矩阵对应位置元素相乘，所以*要求两个矩阵的行、列数互为转置，而.*则要求两个矩阵行、列数要相同；（3）使用eye可以获得单位矩阵函数（矩阵对角线处元素为1，其余元素为0），矩阵的阶数由括号内的值决定，格式为eye(n),n为矩阵阶数。
matlab课程设计完整版
matlab课程设计完整版一、教学目标本节课的教学目标是让学生掌握MATLAB的基本语法和操作，能够利用MATLAB进行简单的数学计算和数据分析。
具体来说，知识目标包括：了解MATLAB的历史和发展，掌握MATLAB的基本语法和数据类型，熟悉MATLAB的工作环境。
技能目标包括：能够使用MATLAB进行矩阵运算，编写简单的MATLAB脚本程序，进行数学计算和数据分析。
情感态度价值观目标包括：培养学生对科学计算软件的兴趣，增强学生的动手能力和团队协作能力。
二、教学内容本节课的教学内容主要包括MATLAB的基本语法和操作。
首先，介绍MATLAB的历史和发展，使学生对MATLAB有一个整体的认识。
然后，讲解MATLAB的基本语法和数据类型，如矩阵的创建和操作，数据的输入和输出等。
接着，介绍MATLAB的工作环境，包括命令窗口、变量浏览器和脚本文件等。
最后，通过实例演示和练习，使学生能够熟练使用MATLAB进行简单的数学计算和数据分析。
三、教学方法为了达到本节课的教学目标，将采用讲授法、实践法和讨论法等多种教学方法。
首先，通过讲授法向学生介绍MATLAB的基本概念和语法。
然后，通过实践法，让学生动手操作MATLAB软件，进行实际的数学计算和数据分析。
在实践过程中，引导学生进行讨论，分享自己的心得和经验，互相学习和进步。
最后，通过讨论法，对学生的学习情况进行总结和评价，及时调整教学策略。
四、教学资源为了保证本节课的教学质量，将准备教材、多媒体资料和实验设备等多种教学资源。
教材是学生学习的基础，多媒体资料可以丰富教学手段，实验设备则是学生进行实践操作的重要工具。
此外，还将利用网络资源，如在线教程和讨论区，为学生提供更多的学习资料和实践机会。
五、教学评估本节课的教学评估将采用多元化的评价方式，以全面、客观、公正地评估学生的学习成果。
评估方式包括平时表现、作业和考试等。
平时表现主要考察学生的课堂参与度和团队合作能力，通过观察和记录学生在课堂上的表现来进行评估。
MATLAB软件的程序设计
构思程序结构
1. Read a, b, and c 2. delt1 ← b^2 - 4*a*c 3. if delt1 > 0 x（1） ← (-b + sqrt (delt1)) / (2*a) x（2） ← (-b – sqrt (delt1)) / (2*a) Write msg that equation has two distinct real roots. elseif delt1 == 0 x（1） ← -b / (2*a) Write msg that equation has two identical real roots. else real_part ← -b / (2*a) imag_part ← sqrt (abs (discriminant)) / (2*a) Write msg that equation has two complex roots. end
将子程序放 回到主程序
用MATLAB 语言实现
1/3 : plan what to do. 1/6 : write the program. 1/2 : test and debug.
合成
联合调试
测试程序运 行结果
正式测试
3.2 流程控制语句
分支语句 Branches
if 语句 switch 语句 try/catch 语句
25
铁心平均长度 + i g 铁心 绕ntinue
结束循环 continue: 直接进行下一次循环，不执行本次循环体 余下的语句
break:
27
3.4 数据导入和导出
input
file->import data
28
3.4 数据导入和导出
南邮通达《Matlab编程技术实验》期末试卷A
南京邮电大学通达学院 2015/2016学年第 二 学期《Matlab 编程技术实验》期末试卷 A本试卷共4页； 考试时间 110 分钟；专业 班级 学号 姓名1. 在桌面新建一个文件夹，命名为“你的学号+姓名”，如：学号为14110101，姓名为张三，则新建的文件夹名为：14110101张三；2. 将word 版试卷拷贝到“你的学号+姓名”的文件夹内，并改名为“你的学号+姓名.doc”，把解题代码和答案保存在此word 文档相应题目下方；3. 请把解答每道题过程中使用的Matlab 程序代码分别用一个M 文件保存，每个M 文件的文件名命名为T+题号。
如保存第一题的代码的M 文件为T1.m ，T1为该M 文件的文件名，.m 为M 文件的后缀。
最后把所有M 文件也保存在以你的学号+姓名命名的文件夹内。
一、简单计算题（60分，每小题10分）(以下各题请写出解题过程所用的代码和结果并保存)1、已知数组x=[9,27,15,0,29,7,4,17,39,117]，请找出数组x 中大于12的元素及这些元素所对应的位置。
（10分）2、已知矩阵A=[1,3,4; 2,6,7; 1,9,5]和矩阵B=[4,2,6;4,9,3;5,7,10]，请计算出A+B, A .*B, A .^3，4*B 。
（10分）3、已知矩阵C=[4,2,0;0,3,1;0,0,6],请计算出矩阵C 的特征值和特征向量。
（10分）守试则试作装 订 线 内 不 要 答 题自 觉遵 考 规 ，诚 信 考 ，绝 不 弊4、创建符号对象计算01cos 4lim 4x x y x →-=和2221lim(1)1n n z n →∞=++（10分）5、有如下联立方程组,请先根据系数矩阵与未知数的关系判断是否存在唯一确定解，并求出方程组的解。
（10分）34712457423859652108x y z w x y z w x z w x y z w +--=⎧⎪-++=-⎪⎨+-=⎪⎪-+-+=-⎩6、创建符号对象x , 求积分211x R dx e ⎛⎫=- ⎪+⎝⎭⎰。
南邮MATLAB数学实验答案(全)[精品文档]
![南邮MATLAB数学实验答案(全)[精品文档]](https://img.taocdn.com/s3/m/3d22b17eb307e87101f6963b.png)
第一次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。
补充命令vpa(x,n) 显示x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin limx mx mxx →∞- syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2) ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100xy e dxdy +⎰⎰dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.13941.4 计算4224x dx m x +⎰ syms xint(x^4/(902^2+4*x^2)) ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最高次幂为4).syms xtaylor(sqrt(902/1000+x),5,x) ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。
matlab教程(完整版)
01 MATLABChapterMATLAB简介MATLAB是一种高级编程语言和环境，主要用于数值计算、数据分析、信号处理、图像处理等多种应用领域。
MATLAB具有简单易学、高效灵活、可视化强等特点，被广泛应用于科研、工程、教育等领域。
MATLAB提供了丰富的函数库和工具箱，方便用户进行各种复杂的数学计算和数据分析。
MATLAB安装与启动MATLAB界面介绍工作空间用于显示当前定义的所有变量及其值。
命令历史记录了用户输入过的命令及其输出结果。
基本运算与数据类型02矩阵运算与数组操作Chapter01020304使用`[]`或`zeros`、`ones`等函数创建矩阵创建矩阵使用`size`函数获取矩阵大小矩阵大小通过下标访问矩阵元素，如`A(i,j)`矩阵元素访问使用`disp`或`fprintf`函数显示矩阵信息矩阵信息矩阵创建与基本操作对应元素相加，如`C = A+ B`加法运算矩阵运算对应元素相减，如`C = A-B`减法运算数与矩阵相乘，如`B = k *A`数乘运算使用单引号`'`进行转置，如`B = A'`转置运算满足乘法条件的矩阵相乘，如`C = A * B`矩阵乘法使用`inv`函数求逆矩阵，如`B = inv(A)`逆矩阵数组创建数组大小数组元素访问数组操作数组操作01020304线性方程组求解数据处理与分析特征值与特征向量图像处理矩阵与数组应用实例03数值计算与数据分析Chapter数值计算基础MATLAB基本运算数值类型与精度变量与表达式函数与脚本数据分析方法数据导入与预处理学习如何导入各种格式的数据（如Excel、CSV、TXT等），并进行数据清洗、转换等预处理操作。
数据统计描述掌握MATLAB中数据统计描述的方法，如计算均值、中位数、标准差等统计量，以及绘制直方图、箱线图等统计图表。
数据相关性分析学习如何在MATLAB中进行数据相关性分析，如计算相关系数、绘制散点图等。
matlab大学课程设计
matlab 大学课程设计一、课程目标知识目标：1. 掌握MATLAB的基本操作，包括数据类型、矩阵运算、数组操作等；2. 学习MATLAB编程技巧，能够编写简单的程序解决问题；3. 了解MATLAB在工程领域的应用，如数值计算、信号处理、图像处理等；4. 理解MATLAB与其他编程语言的区别及优势。
技能目标：1. 能够运用MATLAB进行数据可视化，绘制各种图表；2. 掌握MATLAB的文件操作，包括读写文件、导入导出数据等；3. 学会使用MATLAB进行科学计算和仿真，解决实际问题；4. 培养良好的编程习惯，编写高效、可读性强的代码。
情感态度价值观目标：1. 培养学生对MATLAB编程的兴趣，激发学习热情；2. 增强学生的团队协作能力，通过合作完成课程设计项目；3. 培养学生独立思考、解决问题的能力，提高自信心；4. 引导学生认识到编程在工程领域的应用价值，为未来职业发展奠定基础。
课程性质：本课程为大学本科阶段的一门实践性课程，旨在帮助学生掌握MATLAB编程技能，提高解决实际问题的能力。
学生特点：学生具备一定的编程基础，但对MATLAB编程尚不了解，需要通过本课程学习，逐步掌握相关技能。
教学要求：教师应注重理论与实践相结合，通过讲解、演示、练习和项目实践等多种方式，帮助学生达到课程目标。
同时，关注学生的个体差异，给予有针对性的指导，确保每个学生都能在学习中取得进步。
在教学过程中，将课程目标分解为具体的学习成果，便于后续的教学设计和评估。
二、教学内容1. MATLAB基础知识- 数据类型与变量- 矩阵与数组操作- 基本数学运算与函数- 程序控制结构2. MATLAB编程技巧- 函数与脚本- 逻辑判断与循环结构- 文件操作- 错误处理与调试3. MATLAB数据可视化- 二维图形绘制- 三维图形绘制- 图形修饰与动画制作- 数据导出与报告生成4. MATLAB应用实例- 数值计算- 信号处理- 图像处理- 优化算法5. 课程设计项目- 项目要求与分组- 项目实施与指导- 项目成果展示与评价教学内容安排与进度：第一周：MATLAB基础知识第二周：MATLAB编程技巧第三周：MATLAB数据可视化第四周：MATLAB应用实例第五周：课程设计项目启动与分组第六周：课程设计项目实施与中期检查第七周：课程设计项目成果展示与评价教学内容与教材关联性：本课程教学内容依据教材章节进行组织，确保学生能够循序渐进地掌握MATLAB编程技能。
matlab软件设计实验报告(各源程序均运行无误)
南京邮电大学通信学院软件课程设计实验报告模块名称：__ MATLAB软件设计___专业班级：_____ 通信2班______姓名：___ _ ______学号：____ __ _____实验日期： 2012 年 3月5 日实验报告日期： 2012年 3月8 日1．在时间区间 [0,10]中，绘出t e y t2cos 15.0--=曲线。
程序close all clc clear t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 运行结果：2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序：close all clc cleart=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; 然后回答以下问题： 程序：close all clc clearA=[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)V(:,1)'*null(A)U(:,1:2)=orth(A)（1）sum(diag(S)>tol) - rank(A) 是多少？ans =（2）S(1,1) - norm(A) = =0 是多少?ans =1（3）sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么?ans =（4）S(1,1)/S(3,3) - cond(A) 的结果是什么?ans =（5）S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ans =1（6）V(:,1)'*null(A) 得到什么结果?ans =-5.5511e-017（7）abs(A*null(A)) < tol 得到什么结果?ans =111（8）U(:,1:2) = = orth(A) 的运行结果是什么?ans =1 11 11 14. 求积分dt t x⎰]sin[cos ，π20≤≤x 。
matlab软件设计PPT课件
控制系统的分析和优化
控制系统稳定性分析
利用Matlab对控制系统进行稳定性 分析，判断系统是否稳定，并优化系 统参数。
控制系统性能优化
通过Matlab对控制系统进行性能优化 ，提高系统的响应速度、减小超调和 调节时间等。
控制系统的设计和实现
控制系统设计
使用Matlab进行控制系统设计，包括开环和闭环系统设计，以及PID控制、模 糊控制等先进控制策略。
强大的支持。
插值与拟合
03
Matlab中的polyfit和polyval函数可用于多项式插值和拟合，提
高数值计算的精度。
优化算法
线性规划
Matlab中的linprog函数可以用来 求解线性规划问题，通过最小化 目标函数来找到最优解。
非线性规划
Matlab中的fmincon函数可以用 来求解非线性规划问题，通过迭 代搜索找到最优解。
案例二：机器学习算法应用
机器学习算法概述
Matlab实现步骤
案例展示
机器学习是人工智能的一个重 要分支，通过从数据中学习并 自动发现模式，让机器能够做 出预测和决策。
在Matlab中，我们可以使用内 置的机器学习工具箱来实现各种 机器学习算法。例如，我们可以 使用k-最近邻（k-NN）、决策 树、支持向量机（SVM）等分 类器来对数据进行分类，或者使 用聚类算法来对数据进行聚类。
总结词
列举Matlab的应用领域。
详细描述
Matlab被广泛应用于科学计算、工程设计、数据分析和机器学习等领域。在科研领域，Matlab被用 于模拟、建模和数据分析等方面。在工程领域，Matlab被用于控制系统设计、信号处理和图像处理等 方面。此外，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教程(全)09Matlab程序设计
2020/11/8
Application of Matlab Language
14
5.5 Matlab矩阵分析与处理
5.5.1 特殊矩阵 常见的特殊矩阵有零矩阵、幺矩阵、单位矩阵等，这类特殊矩阵在应用
中具有通用性。 1、通用的特殊矩阵 常用的产生通用殊矩阵的函数有： zeros：产生全0矩阵（零矩阵）。 ones： 产生全1矩阵（幺矩阵）。 eye： 产生单位矩阵。 rand：产生0~1间均匀分布的随机矩阵。 randn：产生均值为0，方差为1的标准正态分布随机矩阵。
5.4.1 程序调试概述 一般说来，应用程序的错误有两类，一类是语法错误，另一类是运行时
的错误。语法错误，给出相应的错误信息，并标出错误在程序中的行 号。例如：输入下列程序： A = 87；
B = 9.3;
C = A+*B; 系统将给出错误信息：
??? Error: File: Untitled1.m Line: 3 Column: 7
2020/11/8
Application of Matlab Language
4
说明：
将以上函数文件以文件名fcircle.m保存，然后在命令窗口调用。
[s,p] = fcircle(10) 输出结果是： s=
314.1593 p=
62.8319 采用help命令或lookfor命令可以显示出注释说明部分的内容。 help fcircle 屏幕显示
进行存取和修改。
全局变量用global命令定义，格式为：
global 变量名
例5.13 全局变量应用示例。
先建立函数文件wadd.m，该函数将输入的参数加权相加：
function f = wadd(x,y)
南京邮电大学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 ，并对所得结果作出解释。
南邮MATLAB上机实验二
南邮MATLAB上机实验二实验二 MATLAB 的符号运算符号计算的特点：一，运算以推理解析的方式进行，因此不受计算误差积累问题困扰；二，符号计算，或给出完全正确的封闭解，或给出任意精度的数值解（当封闭解不存在时）；三，符号计算指令的调用比较简单；四，计算所需时间较长。
在MATLAB 中，符号计算虽以数值计算的补充身份出现，但涉及符号计算的指令使用、运算符操作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便捷的。
一、实验目的1、熟悉符号对象和符号表达式；2、掌握符号表达式和符号函数的操作；3、学习利用符号代数方程的求解 。
二、实验基本知识1、基本运算指令 ；2. 掌握符号运算的应用。
三、上机练习1. （1）通过符号计算求t t y sin )(=的导数dt dy 。
（2）然后根据此结果，求-=0t dt dy 和2π=t dt dy 。
>> clear allsyms ty=abs(sin(t))df=diff(y)df1=limit(df,t,0,'left')df2=subs(df,'t',sym(pi/2))y =abs(sin(t))df =sign(sin(t))*cos(t)df1 =-1df2 =2.计算二重积分⎰⎰+211222)(x dydx y x 。
>>clear allsyms x yf=x^2+y^2Fint=(int(int(f,y,1,x^2),x,1,2))f =x^2 + y^2Fint =1006/1053. 已知22y x f +=(1)求22dy f d ;(2)求⎰-11fdx .>> clear allsyms x yf=x^2+y^2;dfdy2=diff(f,y,2)dfdy2 =2>> intf=int(f,x,-1,1)intf =2*y^2 + 2/34.解方程: 06432=-+x x>> clear allsyms xeq=solve(3*x^2+4*x-6)eq =-2/3+1/3*22^(1/2)-2/3-1/3*22^(1/2)5.已知方程:62333=++y dx dy dxy d (1) 求通解;(2)求满足3)0(,2)0(,1)0(=='=''y y y 的特解.>>clear all>> y=dsolve('D3y+3*Dy+2*y-6=0','x')y =3+C1*exp((1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+C2*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*sin(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x )+C3*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*cos(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)>> y=dsolve('D3y+3*Dy+2*y-6=0','y(0)=3,D2y(0)=1,Dy(0)=2','x')y =3+(-1/12*(1+2^(1/2))^(1/3)*2^(1/2)+1/6*(1+2^(1/2))^(2/3)+1/3*(1+2^(1/2))^(1/3)-1/4*(1+2^(1/2))^(2/3)*2^(1/2))*exp((1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+(1/3*3^(1/2)*(1+2^(1/2))^(1/3)-1/12*3^(1/2)*(1+2^(1/2))^(1/3)*2^(1/2)-1/6*3^(1/2)*(1+2^(1/2))^(2/3)+1/4*3^(1/2)*(1+2^(1/2))^(2/3)*2^(1/2))*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*sin(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+(1/12*(1+2^(1/2))^(1/3)*2^(1/2)-1/6*(1+2^(1/2))^(2/3)-1/3*(1+2^(1/2))^(1/3)+1/4*(1+2^(1/2))^(2/3)*2^(1/2))*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x )*cos(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)6.求方程 221, 2x y xy +== 的解。
南邮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 。
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、求230x e x -=的所有根。
（先画图后求解）（要求贴图）>> ('(x)-3*x^2',0)=-2*(-1/6*3^(1/2))-2*(-11/6*3^(1/2))-2*(1/6*3^(1/2))3、求解下列各题：1）30sin lim x x x x->->> x;>> (((x))^3)=1/62） (10)cos ,x y e x y =求 >> x;>> ((x)*(x),10)=(-32)*(x)*(x)3）21/20(17x e dx ⎰精确到位有效数字）>> x;>> ((((x^2),0,1/2)),17)=0.544987104183622224）42254x dx x+⎰>> x;>> (x^4/(25^2))=125*(5) - 25*x + x^3/35）求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
>> t;>> ((1^2))(t);>> ()()=16）设函数(x)由方程e所确定，求y′(x)。
>> x y;*(y)(1);>> ()()=(x + (y))7）sin2xe xdx+∞-⎰>> x;>> ()*(2*x); >> (y,0)=2/58）08x =展开（最高次幂为）>> x(1);taylor(f,0,9)=- (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 19） 1sin (3)(2)x y e y =求>> x y;>> ((1));>> ((y,3),2)=-0.582610）求变上限函数x⎰对变量x的导数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计专业班级：__通信工程姓名：____ ____学号：___实验日期：2013年 6 月 17—28日实验报告日期： 2013 年 7 月 1 日一、要求练习的实验部分1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。
程序： t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果：2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。
图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。
程序：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) < 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)结果：设 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 。
程序：t=0:pi/100:2*pi; f=abs(sin(cos(t))); F=cumsum(f)*pi/100; plot(t,F);xlabel('x 轴'); ylabel('y 轴');title('描绘y=abs(sin(cos(t)))在(0,x)间积分曲线[x ∈(0,2π)]') grid on; 结果：5. 求方程2,122==+xy y x 的解。
程序： clear; clc;disp('计算方程组 x^2+y^2=1 x*y=2 的根 x y') [x,y] = solve('x^2+y^2 = 1','x*y = 2') 结果：计算方程组 x^2+y^2=1 x*y=2 的根 x y x =[ -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)] y =[ 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=ξ对它们相应的两条曲线进行醒目的文字标志。
程序：b=0.2:0.2:2.0; t=0:0.1:18;color=['r','g.','bo','cx','m+','y*','kv','r--','g:','b-.']; for i=1:10p(i)=sqrt(1-b(i).*b(i));q(i)=atan(sqrt(1-b(i).*b(i))./b(i));y=1-1./p(i).*exp(-b(i).*t).*sin(p(i).*t+q(i)); if i==1plot(t,y,color(i))text(t(30),y(30),' \leftarrow ξ=0.2','FontSize',18) hold onelse if i==10plot(t,y,color(i))text(t(50),y(50),' \leftarrow ξ=2.0','FontSize',18) hold on elseplot(t,y,color(i)) hold on end end end 结果：7.构建一个简单的全波整流模型，并用示波器分两路同时观察原信号和整流后的信号波形。
要求：信源模块被重新命名为“输入正弦波”；信宿模块被重新命名为“示波器”；连接到信宿的信号线上分别标注“原信号”和“整流信号”；在模型窗中添加注释文本。
(1)、原理图：各模块参量设置：[示波器]ParametersNumber of axes: 2其他默认[其余模块]默认选项结果(2)、原理图：各模块参数设定[Gain]Gain: -1其他默认[其余模块]默认选项结果8.利用SIMULINK及其标准模块设计一个低通滤波器，从受噪声干扰的多频率混合信号π，而)(tx中获取10Hz的信号。