南京邮电大学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编程技能。
- 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的信号。