matlab作业5答案

MATLAB

1、 试求出下面线性微分方程的通解。

dy(t)dt55作业5

?13dy(t)dt44?64dy(t)dt33?152dy(t)dt22?176dy(t)dt?80y(t)?e?2t[sin(2t??3)?cos(3t)]&&假设上述微分方程满足已知条件，y(0)?1,y(1)?3,y(?)?2,y(0)?1,y(1)?2试

求出满足该条件的微分方程的解析解。

解： >> syms t y ;

y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=','exp(-2*t)*(sin(2*t+pi/3)+cos(3*t))'],'y(0)=1','y(1)=3','y(pi)=2','Dy(0)=1','Dy(1)=2'); vpa(y,20) ans =

.20576131687242798354e-2*exp(-2.*t)*cos(3.*t)+.15538705805619602373e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035587e-2*exp(-2.*t)*cos(2.*t)+98.

MATLAB

1、 试求出下面线性微分方程的通解。

dy(t)dt55作业5

?13dy(t)dt44?64dy(t)dt33?152dy(t)dt22?176dy(t)dt?80y(t)?e?2t[sin(2t??3)?cos(3t)]&&假设上述微分方程满足已知条件，y(0)?1,y(1)?3,y(?)?2,y(0)?1,y(1)?2试

求出满足该条件的微分方程的解析解。

解： >> syms t y ;

y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=','exp(-2*t)*(sin(2*t+pi/3)+cos(3*t))'],'y(0)=1','y(1)=3','y(pi)=2','Dy(0)=1','Dy(1)=2'); vpa(y,20) ans =

.20576131687242798354e-2*exp(-2.*t)*cos(3.*t)+.15538705805619602373e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035587e-2*exp(-2.*t)*cos(2.*t)+98.

MATLAB

作业5

1、 试求解下面微分方程的通解以及满足x(0)?1,x(?)?2,y(0)?0条件下的解析解。

????(t)?4x(t)?3y(t)?e?6tsin(4t)x(t)?5x ? ?6t??2y(t)?y(t)?4x(t)?6x(t)?ecos(4t)?解：

>> syms t x y;

>>[x,y]=dsolve('D2x+5*Dx+4*x+3*y=exp(-6*t)*sin(4*t)',

'2*Dy+y+4*Dx+6*x=exp(-6*t)*cos(4*t)','x(0)=1','x(pi)=2','y(0)=0');vpa(x,20)

ans =

.24816266536011758942e-1*exp(-6.*t)*cos(4.*t)-.16682998530132288094e-1*exp(-6.*t)*sin(4.*t)+.85861668591377882749e-1*exp(t)-.57658489325155296910e-1*exp(-5.1374586088176874243*t)+.9469805541977656

MATLAB作业1参考答案

1 、启动MATLAB 环境，并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc，试运行该语句，观察得出的结果，并利用help 命令对你不熟悉的语句进行帮助信息查询，逐条给出上述程序段与结果的解释。

【求解】在MATLAB 环境中感触如下语句，则可以看出，求解500?500随机矩阵的逆，并求出得出的逆矩阵与原矩阵的乘积，得出和单位矩阵的差，得出范数。一般来说，这样得出的逆矩阵精度可以达到10?12。

>> tic, A=rand(500); B=inv(A); norm(A*B-eye(500)), toc ans = 1.2333e-012

Elapsed time is 1. seconds.

2、试用符号元素工具箱支持的方式表达多项式f(x)?x5?3x4?4x3?2x2?3x?6，并令

x?s?1，将f(x) 替换成s 的函数。 s?1【求解】可以先定义出f 函数，则由subs() 函数将x 替换成s 的函数 >> syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1)) F =

MATLAB作业1参考答案

1 、在你的机器上安装MATLAB 语言环境，并键入demo 命令，由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。 【求解】略.

2、启动MATLAB 环境，并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc，试运行该语句，观察得出的结果，并利用help 命令对你不熟悉的语句进行帮助信息查询，逐条给出上述程序段与结果的解释。

【求解】在MATLAB 环境中感触如下语句，则可以看出，求解500?500随机矩阵的逆，并求出得出的逆矩阵与原矩阵的乘积，得出和单位矩阵的差，得出范数。一般来说，这样得出的逆矩阵精度可以达到10?12。

>> tic, A=rand(500); B=inv(A); norm(A*B-eye(500)), toc ans =

1.2333e-012

Elapsed time is 1.301000 seconds.

3、试用符号元素工具箱支持的方式表达多项式f(x)?x?3x?4x?2x?3x?6，并令

x?s?1s?15432，将f(x) 替换成s 的函数。

【求解】可以先定义出f 函数

MATLAB作业

学号：07082002 姓名：武婷

第一次作业：

1) 熟悉Matlab运行环境;简介一下你都做了什么!

首先，安装了Matlab软件，并启动该软件。

Matlab启动后，即出现Desktop操作桌面，为一个高度集成的工作界面。同时发现有四个界面：指令窗（command window),当前目录浏览器（current

directory），历史指令窗（command history）,工作空间浏览器（workspace brower）。

指令窗（command window)：在该窗内，进行了简单的算术运算，矩阵的不同输入方式的体验，既简单的函数调用等。如在命令窗口中键入一些简单的命令：rand（4，4）生成4*4其元素在0-1之间的随机矩阵，norm(x,1)求向量的1-范数。

再如输入指令：>> 7*(12^2+66)/21，按下回车键即可得到结果： ans = 5.6084e+070

当前目录浏览器：为管理方便，首先建立用户目录并将它设置成当前目录，了解到了有两种方法;交互界面设置法和指令设置法。同时基于对该软件搜索路径知识的学习利用设置路径对话框修改搜索路径，并尝试了两种修改状态：当前

1、 与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 2、 存储在工作空间中的数组能编辑吗？如何操作？

3、 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 4、 Matlab提供了两种运行方式，它们是？ 5、 在MATLAB中有几种获得帮助的途径？

6、 有几种建立矩阵的方法？ 7、 在MATLAB中如何建立矩阵??5?4793??1?，并将其赋予变量a？

8、 试用matlab语句生成一个6行5列的矩阵，要求其中的元素全为0 9、 matlab变量命名的规则是？

9、 下列符号哪些是matlab的合法变量名？

3vars、global、help、My_exp、sin、X+Y、_input、 E-4 、AbCd 、 AB_C_

10、若不想显示命令运行的结果，应怎么做？若要注释掉一行，应该怎么做？ 11、请说出下面几个命令的含义：clc，clear, which 12、数组运算和矩阵运算的运算符有什么区别？ 13、求x?4?8i???3?2i3?5i7?6i2?7i9?4i1?4i3?9i7?5i??4?4i?的共轭转置

14、请用冒号语句生成行向量，要求其内容为1到100 ?5?15、有矩阵a?3???73793?

多幅度电平信号差错率仿真设计

一、设计任务

完成多幅度电平信号差错率的仿真，运用MATLAB语言，仿真它们的误码性能，并与理论误码曲线相比较。先产生多组符号序列，再将符号序列映射到对应的幅度电平上，以四电平为例，先产生一个四组符号的序列，将该序列映射到对应的幅度电平上。随机变量信道中的加性噪声可以认为是统计独立，均值为0的高斯分量，分别作出二、四、八、十六电平仿真得到的误码率曲线与和理论上的误码率曲线。

二、物理模型

多幅度电平信号仿真系统模型如下图所示：

出错计数器Am与Am比较^均匀随机数发生器X映射到幅度电平Am+r2检测器(0,?)高斯随机数发生器Am^

图1

三、预期目标

通过以上的物理模型，作出二、四、八、十六电平信号仿真得到的误码率

曲线与和理论上的误码率曲线。 四、系统中所涉及的公式如下：

1、平均比特SNR定义为：

即可得：

1d2 二电平时平均比特SNR为：SNR?

2?25d2 四电平时平均比特SNR为：SNR? 24?21d2 八电平时平均比特SNR为：SNR?

6?2

十六电平时平均比特SNR为:

SNR?85d28?2

2、M电平PAM系统最佳

m a t l a b综合大作业(附详细

答案)

标准化文件发布号：（9312-EUATWW-MWUB-WUNN-INNUL-DQQTY-

《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：

1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)

实验结果：

A =

2)将矩阵A按列拉长得到矩阵B；

实验程序：B=A(:)

实验结果：

B =

3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的

矩阵C；

实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]

实验结果：

C =

4)寻找矩阵A中大于0的元素；]

实验程序：G=A(find(A>0))

实验结果：

G =

5)求矩阵A的转置矩阵D；

实验程序：D=A'

实验结果：

D =

6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；

实验程序：E=flipud(fliplr(A))

实验结果：

E =

7)删除矩阵A的第2列和第4列得到矩阵F；

实验程序：F=A;

F(:,[2,4])=[]

实验结果：

F =

8)求矩阵A的特征值和特征向量；

实验程序：[Av,Ad]=eig(A)

实验结果：

特征向量Av =

+ -

+

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后