通信系统计算机仿真上机实验报告,计算机仿真实验报告实验..docx

《计算机仿真》上机实验报告

姓名：

学号：2012104021

专业：测控

班级：

12 级

实验一常微分方程的求解及系统数学模型的转换

一．实验目的

通过实验熟悉计算机仿真中常用到的 Matlab 指令的使用方法，掌握常微分方程求解指令和模型表示及转换指令， 为进一步从事有关仿真设计和研究工作打下基础。

二. 实验设备

个人计算机， Matlab 软件。

三. 实验准备

预习本实验有关内容(如教材第 2、3、5 章中的相应指令说明和例题) ，编写本次仿真练习题的相应程序。

四. 实验内容

Matlab中常微分方程求解指令的使用

题目一：请用 MATLAB的ODE45算法分别求解下列二个方程。 要求 :1. 编写出 Matlab 仿真程序； 2. 画出方程解的图形并对图形进行简要分析； 3. 分析下列二个方程的关系。

1．2．

1.function fun=funl(t,x)

fun=-x^2;

[t,x]=ode45('fun1',[0,20],[1]);

figure(1);plot(t,x);

grid

2.function fun=fun2(t,x)

fun=x^2;

[t,x]=ode45('fun2',[0,20],[-1]);

figure(2);plot(t,x);

grid

题目二：下面方程组用在人口动力学中，可以表达为单一化的捕食者- 被捕食者

模式(例如，狐狸和兔子)。其中x1 表示被捕食者，x2 表示捕食者。如果被捕

食者有无限的食物，并且不会出现捕食者。于是有x1'x1 ，则这个式子是以指数

形式增长的。 大量的被捕食者将会使捕食者的数量增长； 同样，越来越少的捕食者会使被捕食者的数量增长。 而且，人口数量也会增长。 请分别调用 ODE45、ODE23算法求解下面方程组。 要求编写出 Matlab 仿真程序、画出方程组解的图形并对图形进行分析和比较。

1.ODE45

function fun=fun3(t,x)

fun=[x(1)-0.1*x(1)*x(2)+0.01*t;-x(2)+0.02*x(1)*x(2)+0.04*t];

[t,x]=ode45('fun3',[0,20],[30,20]);

plot(t,x);

Grid

2.ODE23

function fun=fun3(t,x)

fun=[x(1)-0.1*x(1)*x(2)+0.01*t;-x(2)+0.02*x(1)*x(2)+0.04*t];

[t,x]=ode23('fun3',[0,20],[30,20]);

plot(t,x);

Grid

Matlab 中模型表示及模型转换指令的使用

题目三：若给定系统的的传递函数为

6s3

12s2

6s

10

G (s)

4

2s3

3s2

1s 1

s

请用 MATLAB编程求解其系统的极零点模型。

num=[6 12 6 10];

den=[1 2 3 1 1];

sys=tf(num,den);

H=zpk(sys)

题目四：习题 2.4

系统传递函数为

G(s)

s2

4s

5

求其对角标准型实现。

(s 1)(s

2

5s 6)

num=[1 4 5];

den=[1 6 11 6];

sys=tf(num,den);

canon(sys,'modal')

题目五：习题 5.8

D ( z)

z

0.03z 0.02 ，采样周期 Ts=0.02s, 现需要在 T=0.1s

已知某离散系统

z2

下作数字仿真，求该系统的数字仿真模型。

D=zpk(0,[0.1 0.2],1,0.1)

D2=d2d(D,0.02)

五．总结与体会

初步掌握了 Matlab 常用指令的使用方法，为下一步学习打下了基础。

实验二Matlab 优化工具箱的使用

一．实验目的

通过上机操作熟悉 Matlab 优化工具箱的主要功能及其使用方法，掌握优化工具箱中常用函数的功能和语法， 并利用其进行极值运算、 求解线性和非线性问题等，为进一步的仿真设计和研究打下基础。

二. 实验设备

个人计算机， Matlab 软件。

三. 实验准备

预习本实验有关内容 (如教材第 6 章中的相应指令说明和例题) ，编写本次仿真练习题的相应程序。

四. 实验内容

应用 Matlab 优化工具箱求解优化问题

例题 6.6~6.10，选做 2 题，要求自行修改方程系数，并比较运行结果。

2x1x2e x1

例 6.6

例如求解非线性方程：

x1 2x2

e

x2

初始值为 x0 = [-5 -5] 。

原程序

function F=myfun(x)

F=[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];

x0=[-5,-5];

op