南京邮电大学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语言平台与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软件设计南京邮电大学通信学院软件课程设计实验报告模块名称___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程序设计的基本概念和操作步骤。

首先，我们需要了解Matlab的基本语法和编程元素。

Matlab 程序由语句和函数构成，语句是实现具体功能的基本单位，函数则是对语句进行封装和重用的机制。

Matlab采用类似于C 语言的结构化编程风格，包括顺序、选择和循环三种基本的控制结构。

此外，Matlab还支持向量化操作和矩阵运算，可以大大简化程序的编写和运算的复杂度。

其次，我们需要了解Matlab的开发环境和常用工具。

Matlab 提供了一个交互式的命令窗口，用户可以直接在命令行输入和执行Matlab语句。

除此之外，Matlab还提供了一个集成的开发环境，可以方便地编辑、调试和运行Matlab程序。

开发环境中还包括了各种工具箱，可以扩展Matlab的功能，如信号处理、图像处理、控制系统等。

接下来，我们将介绍Matlab程序的编写和调试过程。

在编写Matlab程序时，我们首先需要确定程序的功能和数据结构，然后选择合适的函数和语句来实现。

编写完成后，我们可以在命令窗口中运行程序，观察程序的输出结果是否满足预期。

如果发现程序出现了错误或不符合预期，我们可以使用Matlab 提供的调试工具来进行调试和排错。

调试工具包括断点、单步执行、变量监视等功能，可以帮助我们分析程序的执行过程和找到问题所在。

最后，我们需要注意Matlab程序的性能和代码优化。

Matlab是一种解释性语言，相对于编译型语言来说，其执行速度较慢。

为了提高程序的性能，我们可以使用向量化操作和矩阵运算来减少循环和临时变量的使用。

此外，我们还可以使用Matlab提供的性能分析工具来找出程序的瓶颈，并进行相应的算法优化和代码改进。

综上所述，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解决实际问题的能力。

具体的教学目标如下：1.知识目标：（1）理解Matlab的基本语法和编程方法。

（2）熟悉Matlab的数据类型、运算符和控制结构。

（3）掌握Matlab函数的定义和调用方法。

（4）了解Matlab的绘图功能和图形用户界面设计。

2.技能目标：（1）能够运用Matlab编写简单的程序，解决数学、工程和科学计算问题。

（2）具备阅读和分析Matlab程序的能力。

（3）能够独立完成Matlab程序的调试和优化。

（4）学会使用Matlab进行数据分析和可视化。

3.情感态度价值观目标：（1）培养学生对计算机编程的兴趣和热情。

（2）增强学生运用计算机技术解决实际问题的意识。

（3）培养学生团队合作精神和自主学习能力。

二、教学内容本课程的教学内容主要包括以下几个部分：1.Matlab概述：介绍Matlab的发展历程、功能特点和应用领域。

2.Matlab基本语法：包括数据类型、运算符、控制结构等。

3.Matlab函数：介绍常用数学函数、绘图函数和图形用户界面设计函数。

4.Matlab编程实践：通过实例讲解如何运用Matlab解决实际问题。

5.Matlab与其他软件的接口：介绍Matlab与其他软件的联合使用方法。

三、教学方法为了提高教学效果，本课程将采用多种教学方法相结合的方式进行教学：1.讲授法：用于讲解Matlab的基本语法和函数。

2.案例分析法：通过分析实际案例，使学生掌握Matlab在工程和科学计算中的应用。

3.实验法：安排上机实验，让学生动手实践，巩固所学知识。

4.讨论法：学生进行小组讨论，培养团队合作精神和解决问题的能力。

四、教学资源为了支持本课程的教学，我们将准备以下教学资源：1.教材：选用权威、实用的Matlab教材作为主要教学资料。

2.参考书：提供相关的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 。

软件设计报告姓名：。

学号：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程序设计及应用课程设计一、课程设计背景和目的MATLAB是一种广泛应用于各个领域的计算机语言，是许多科学和工程领域的标准工具之一，其强大的数学计算能力和图像处理能力使其在科学计算、数据分析、仿真和控制系统设计等方面拥有广泛的应用。

本课程设计旨在通过学生的实践操作，运用MATLAB软件进行程序开发和相关应用，掌握MATLAB的基本操作和应用技能，培养学生的实践能力，提高MATLAB应用能力和解决问题的能力。

二、课程设计任务1. 熟悉MATLAB基本操作在学习前，学生需要对MATLAB软件的基本操作有所了解。

在课程设计开始之前，需要学生自学MATLAB的基本操作。

2. 设计并实现简单的程序在学习MATLAB的基本操作后，学生需要设计并实现一个简单的程序，例如实现一个线性函数的计算。

当程序设计完成后，需要生成程序的图像，并生成程序的代码。

3. 设计并实现高级程序在掌握基本程序设计技能的基础上，学生需要设计并实现一个较为复杂的程序，例如实现一个数字信号的滤波。

当程序设计完成后，需要生成程序的图像和代码，并进行程序的性能评估。

4. 应用MATLAB进行信号处理在以上任务完成后，学生需要根据之前学习的知识和技能，结合MATLAB进行信号处理的实践应用。

例如根据余弦函数的变化规律设计程序，完成余弦信号的生成和分析。

三、课程设计实施过程1. 阶段一：自学MATLAB基本操作在开课前，要求学生阅读相关学习资料，自学MATLAB的基本操作方法和语法。

2. 阶段二：设计并实现简单的程序在第一阶段完成后，学生需要设计并实现一个简单的程序。

教师和助教将针对学生的程序进行代码审查和基本语法的指导，让学生更好地理解和掌握MATLAB语言的特点和使用方法。

3. 阶段三：设计并实现高级程序在简单程序实现完成后，学生需要进一步设计并实现一个较为复杂的程序。

教师和助教将在本阶段对学生的代码进行查错和程序性能的分析，帮助学生更加清晰地理解程序框架和算法原理。

南京邮电大学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。

注意:在下面的题目中m 为你的学号的后4位第一次练习题1． 求032=-x e x 的所有根。

（先画图后求解）2． 求下列方程的根。

1)0155=++x x 2） 至少三个根）(021sin =-x x3） 所有根0cos sin 2=-x x x3． 求解下列各题： 1）3sin lim x mx mx x->- 2） )10(,cos y x e y x求=3）21/20(17mxe dx⎰精确到位有效数字） 4）424x dx m x+⎰ 5） 081000mx x +=将在展开（最高次幂为）6） 1sin (3)()xy e y m =求( 精确到17位有效数字)4． 1)求矩阵21102041100A m ⎛⎫ ⎪- ⎪= ⎪ ⎪ ⎪-⎝⎭ 的逆矩阵1-A 及特征值和特征向量。

2)求点(1,1,4)到直线l: (x-3)/-1 =y/0=(z+1)/2的距离。

5． 已知,21)(222)(σμσπ--=x ex f 分别在下列条件下画出)(x f 的图形：）；（在同一坐标系上作图，，＝时＝、）；（在同一坐标系上作图，－，＝时、421,0)2(110,1)1(σμμσ=、6． 画 下列函数的图形：（1）202004cos sin ≤≤≤≤⎪⎪⎩⎪⎪⎨⎧===u t t z tu y t u x(2) 30,30)sin(≤≤≤≤=y x xy z（3）sin (3cos )02cos (3cos )02sin x t u t y t u u z u ππ=+⎧≤≤⎪=+⎨≤≤⎪=⎩（第6题只要写出程序）. 第二次练习题1、 设11()/23n nnm x x x x +⎧=+⎪⎨⎪=⎩，数列}{n x 是否收敛？若收敛，其值为多少？精确到6位有效数字。

2、设 ,131211p p p n nx ++++=Λ }{n x 是否收敛？若收敛，其值为多少？精确到17位有效数字（提示：当n x 与1n x +的前17位有效数字一致时终止计算）注：学号为单号的取7=p ，学号为双号的取.8=p 书上习题：(实验四)1，2，4，7（1），8，12（改为：对例2，取 120,55,25,5.4=a 观察图形有什么变化．），13，14 。

南邮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 +== 的解。