MATLAB 练习(含代码)

MATLAB 练习

实验一常见分布的概率密度、分布函数生成

[实验目的]

1. 会利用MATLAB软件计算离散型随机变量的概率，连续型随机变量概率密度值。

2.会利用MATLAB软件计算分布函数值，或计算形如事件。

3.会求上分位点以及分布函数的反函数值。

[实验要求]

1.掌握常见分布的分布律和概率密度的产生命令，如binopdf,normpdf

2. 掌握常见分布的分布函数命令，如binocdf,normcdf

3. 掌握常见分布的分布函数反函数命令，如binoinv,norminv

[实验内容]

1 事件A在每次试验中发生的概率是0.3，计算

（1）在10次试验中A恰好发生6次的概率；

（2）在10次试验中A至多发生6次的概率.

binopdf(6,10,0.3)

>> binocdf(6,10,0.3)

2设随机变量X服从参数是3的泊松分布，求概率

poisspdf(6,3)

3设随机变量X服从区间[2,6]上的均匀分布，求

（1）X=4时的概率密度值；

（2）.

unifpdf(4,2,6)

>> unifcdf(5,2,6)

4设随机变量X服从参数是6的指数分布，求

（1）X=0,1,2,3,4,5,6时的概率密度值;

(2).

exppdf(0:6,6)

>> expcdf(5,6)

5设随机变量X服从均值是6，标准差是2的正态分布，求

(1) X=3,4,5,6，7,8,9时的概率密度值;

（2）X=3,4,5,6，7,8,9时的分布函数值；

（3）若=0.345,求x;

（4）求标准正态分布的上0.05分位数。

normpdf(3:9,6,2)

normcdf(3:9,6,2)

norminv(0.345,6,2)

norminv(0.95,0,1)

6设随机变量X服从自由度是6的t分布,求

（1）X=-3,-2,-1,0,1,2,3时的概率密度值; （2）X=-3,-2,-1,0,1,2,3时分布函数值; （3）若=0.345,求x;

（4）求t分布的上0.05分位数.

tpdf(-3:3,6)

tcdf(-3:3,6)

tinv(0.345,6)

tinv(0.95,6)

7设随机变量X服从自由度是6的

2

χ分布,求

(1) X=0,1,2,3,4,5,6时的概率密度值;

(2) X=0,1,2,3,4,5,6时的分布函数值;

(3) 若=0.345,求x;

(4) 求

2

χ分布的上0.05分位数.

chi2pdf(0:6,6)

chi2cdf(0:6,6)

chi2inv(0.345,6)

chi2inv(0.95,6)

8设随机变量X服从第一自由度是2，第,二自由度是6的F分布,求

(1) X=0,1,2,3,4,5,6时的概率密度值;

(2) X=0,1,2,3,4,5,6时的分布函数值;

(3) 若=0.345,求x;

(4) 求F分布的上0.05分位数.

fpdf(0:6,2,6)

fcdf(0:6,2,6)

finv(0.345,2,6)

finv(0.95,2,6)

实验二概率作图

[实验目的]

1.熟练掌握MATLAB软件的关于概率分布作图的基本操作

2.会进行常用的概率密度函数和分布函数的作图

3.会画出分布律图形

1.掌握MATLAB画图命令plot

2.掌握常见分布的概率密度图像和分布函数图像的画法

[实验内容]

9事件A在每次试验中发生的概率是0.3，记10次试验中A发生的次数为X.

(1)画出X的分布律图形；

x=0:10;

>> y=binopdf(x,10,0.3);

>> plot(x,y,'.')

>>

（2）画出X的分布函数图形；

x=0:0.01:10;

>> y=binocdf(x,10,0.3);

>> plot(x,y)

10设随机变量X服从参数是6的指数分布，

（1）画出X的概率密度图形

x=0:0.01:10;

>> y=exppdf(x,6);

>> plot(x,y)

（2）画出X的分布函数图形

x=-1:0.01:10;

>> y=expcdf(x,6);

>> plot(x,y)

11设随机变量X服从参数是3的泊松分布。

(1)画出X的分布律图形；

x=0:10;

>> y=poisspdf(x,3);

>> plot(x,y,'.')

（2）画出X的分布函数图形；

>> x=0:0.01:10;

>> y=poisscdf(x,3);

>> plot(x,y)

12 设随机变量X服从区间[2,6]上的均匀分布。

（1）画出X的概率密度图形

x=0:0.01:10;

>> y=unifpdf(x,2,6);

>> plot(x,y,'*')

（2）画出X的分布函数图形

>> x=0:0.01:10;

>> y=unifcdf(x,2,6);

13设随机变量X服从均值是6，标准差是2的正态分布。

（1）画出X的概率密度图形

x=-10:0.01:10;

>> y=normpdf(x,6,2);

>> plot(x,y)

（2）画出X的分布函数图形

> x=-10:0.01:10;

>> y=normcdf(x,6,2);

>> plot(x,y)

（3）在同一个坐标系中画出均值是6，标准差是1，,3的正态分布概率密度图形x=0:0.01:13;

>> y1=normpdf(x,6,1);

>> y2=normpdf(x,6,2);

>> y3=normpdf(x,6,3);

>> plot(x,y1,x,y2,x,y3)

14设随机变量X服从自由度是6的t分布

（1）画出X的概率密度图形

x=-10:0.01:10;

>> y=tpdf(x,6);

>> plot(x,y)

（2）画出X的分布函数图形

x=-10:0.01:10;

>> y=tcdf(x,6);

>> plot(x,y)

15设随机变量X服从自由度是6的

2 分布

（1）画出X的概率密度图形

x=0:0.01:10;

>> y=chi2pdf(x,6);

>> plot(x,y)

（2）画出X的分布函数图形

x=0:0.01:10;

>> y=chi2cdf(x,6);

>> plot(x,y)

16设随机变量X服从第一自由度是2，第,二自由度是6的F分布（1）画出X的概率密度图形

x=0:0.001:10;

>> y=fpdf(x,2,6);

>> plot(x,y)

（2）画出X的分布函数图形

x=0:0.001:10;

>> y=fcdf(x,2,6);

实验三数字特征

[实验目的]

1 加深对数学期望，方差的理解

2理解数学期望，方差的意义，以及具体的应用

3 加深对协方差，相关系数的理解

4 了解协方差，相关系数的具体的应用

[实验要求]

1 概率与频率的理论知识，MATLAB软件

2 协方差，相关系数的理论知识，MATLAB命令cov,corrcoef

[实验内容]

17若~B(20,0.3), 求E(X),D(X).

[M,,V]=binostat(20,0.3)

18随机变量X的概率密度为，求E(X),D(X).

Syms X

f1=x;

f2=2-x;

Ex=int(xf1,0,1)+int(xf2,1,2);

Ex2=int(x^2f1,0,1)+int(x^2f2,1,2);

Dx=Ex2-Ex^2

19设(X,Y)的概率密度为, 求EX,EY.

Syms x y

fxy=8xy;

Ex=int(int(fxyx,y,0,x),x,0,1)

Ey=int(int(fxyy,y,0,x),x,0,1)

20（续19）求cov（X,Y）

Syms x y

fxy=8xy;

Ex=;

Ey=;

Cxy=int(int(fxy(x-Ex)(y-Ey),y,0,x),x,0,1)

21 某种商品每件表面上的疵点数X服从泊松分布，平均每件上有0.8个疵点。若规定表面不超过一个疵点的为一等品，价值10元，表面疵点数大于一个不多于4个的为二等品，价值8元。表面疵点数多于4个则为废品，求产品价值的均值。

解设X表示产品表面上的疵点数，由已知，EX=0.8,且X服从泊松分布，故EX=0.8, =

设Y表示产品价值，则Y有分布律：

==0.8088

=0.1898

=0.0014

故有EY=0(元)

> pro=[ ];

>> price=[0 10 8];

>> pro(2)=poisscdf(1,0.8);

>> pro(3)=poisscdf(4,0.8)-pro(2);

>> pro(1)=1-pro(2)-pro(3)

>> Ey=pro*price'

22 设随机变量X~N(1,9),Y~N(0,16) ，且X与Y的相关系数为，令Z=X/3+Y/2.求（1）E(Z) ,D(Z) ;

(2) 求X 与Z的相关系数。

解根据题意，有E(X)=1,D(X)=9,E(y)=0 ,D(y)=16

由E(Z)=E(X/3)+E(y/2)=E(X)/3+E(Y)/2得

E(Z)=1/3+0=1/3

由，cov(X,Y)= 有

Cov(X,Y)=0.5=6

D(Z)=D(X/3+Y/2)=D(X/3)+D(Y/2)+2cov(X/3,Y/2)

得到D(Z)=D(X)/9+D(Y)/4+2cov((X,Y)

=9/9+16/4+2=3

得cov（X,Z）=cov(X,X/3)+cov(X,Y/2)

=1/3cov(X,X)+1/2cov(X,Y)

=9/3+(-6)=0

[Ex,Dx]=normstat(1,sqrt(9));

>> [Ey,Dy]=normstat(0,sqrt(16));

>> rxy=-0.5;

>> syms x y z

>> z=x/3+y/2;

>> covxy=rxy*sqrt(Dx)*sqrt(Dy);

>> Ez=Ex/3+Ey/2

>> Dz=Dx/9+Dy/4+2*1/3*1/2*(covxy)

>> covxz=Dx/3+covxy/2

实验四统计中的样本数字特征

23随机生成4组10个整数数据，求每组数据的平均值。

> X=fix(20*rand(10,4))

>> M=mean(X)

24 随机生成服从标准正态分布的6组10个数据，求每组数据的极差、样本方差、样本标准差。

X=normrnd(0,1,10,6)

M1=range(X)

M2=var(X)

M3=std(X)

实验五两个正态总体均值差，方差比的区间估计

[实验目的]

1掌握两个正态总体均值差，方差比的区间估计方法

2会用MATLAB求两个正态总体均值差，方差比的区间估计

[实验要求]

两个正态总体的区间估计理论知识

[实验内容]

25 为比较甲乙两种型号步枪子弹的枪口速度，随机地取甲种型号子弹10发，得到枪口速度的平均值为500（m/s ）,标准差1.10（m/s ）；乙种型号子弹20发，得到枪口速度的平均值为496（m/s ）,标准差1.20（m/s ）。两个总体都近似服从正态分布，且方差相等。求两个总体均值差的置信水平为0.95的置信区间。

t=tinv(0.975,28);

>> s=sqrt((9*1.10^2+1.20^2)/28);

>> d1=(500-496)-t*s*sqrt(1/10+1/20)

>> d2=(500-496)+t*s*sqrt(1/10+1/20)

26 从甲乙两个蓄电池厂生产的产品中，分别抽取10个产品，测得它们的电容量（单位：Ah ）为：

甲厂：146 141 138 142 140 143 138 137 142 137

乙厂：141 143 139 139 140 141 138 140 142 136

若蓄电池的电容量服从正态分布，求两个工厂生产的蓄电池的电容量方差之比的置信水平为0.90的置信区间。

x1=[146 141 138 142 140 143 138 137 142 137];

>> x2=[141 143 139 139 140 141 138 140 142 136];

>> s1=var(x1);

>> s2=var(x2);

>> f1=finv(0.95,9,9);

>> f2=finv(0.05,9,9);

>> d1=s1/s1*(1/f1)

>> d2=s1/s2*(1/f2)

实验六 假设检验

[实验目的]

1 会用MATLAB 进行单个正态总体均值及方差的假设检验

2 会用MATLAB 进行两个正态总体均值差及方差比的假设检验

[实验要求]

熟悉MATLAB 进行假设检验的基本命令与操作

[实验内容]

27 某工厂生产10

Ω的电阻，根据以往生产的电阻实际情况，可以认为其电阻值服从正态分布，标准差=0.1Ω.现随机地抽取10个电阻，测得它们的电阻值为（Ω）：

9.9 ，10.1 ，10.2 ，9.7 ，9.9， 9.9 ，10 ，10.5 ，10.1 ，10.2

问我们能否认为该厂的电阻的平均值为10Ω？（取=0.1）

x=[9.9 10.1 10.2 9.7 9.9 9.9 10 10.5 10.1 10.2];

>> [h,sig,ci]=ztest(x,10,0.1,0.1)

28 用自动包装机包装大米，袋装米的质量服从正态分布。在机器正常工作情况下，其均值

为0.5kg.某天检查机器工作情况，从已经包装好的大米中随机地抽取9袋，称得质量为（单位：kg）：0.497,0.512,0.506,0.515,0.520,0.522,0.498,0.519,0.510

问机器工作是否正常？（取=0.05）

x=[0.497,0.512,0.506,0.515,0.520,0.522,0.498,0.519,0.510];

>> [h,sig,ci]=ttest(x,0.5,0.05)

29 杜鹃总是把蛋生在别的鸟巢中，现在从甲、乙两种鸟巢分别得到杜鹃蛋9只、15只，测得杜鹃蛋的长度数据（单位：mm）如下：

甲鸟巢：21.2,21.6,21.9,22.0,22.0,22.2,22.8,22.9,23.2

乙鸟巢：19.8,20.0,20.3,20.8,20.9,20.9,21.0,21.0,21.1,21.2,21.5,22.0,21.9,21.1,22.3

假设这两个样本来自同方差的正态总体，试鉴别杜鹃蛋的长度差异是由于随机因素造成的，还是与它们被发现的鸟巢不同有关？（取=0.05）

x=[21.2,21.6,21.9,22.0,22.0,22.2,22.8,22.9,23.2];

>> y=[19.8,20.0,20.3,20.8,20.9,20.9,21.0,21.0,21.1,21.2,21.5,22.0,21.9,21.1,22.3];

>> [h,sig,ci]=ttest2(x,y,0.05)

matlab源代码实例

1.硬币模拟试验 源代码： clear; clc; head_count=0; p1_hist= [0]; p2_hist= [0]; n = 1000; p1 = 0.3; p2=0.03; head = figure(1); rand('seed',sum(100*clock)); fori = 1:n tmp = rand(1); if(tmp<= p1) head_count = head_count + 1; end p1_hist (i) = head_count /i; end figure(head); subplot(2,1,1); plot(p1_hist); grid on; hold on; xlabel('重复试验次数'); ylabel('正面向上的比率'); title('p=0.3试验次数N与正面向上比率的函数图'); head_count=0; fori = 1:n tmp = rand(1); if(tmp<= p2) head_count = head_count + 1; end p2_hist (i) = head_count /i; end figure(head); subplot(2,1,2); plot(p2_hist); grid on; hold on; xlabel('重复试验次数'); ylabel('正面向上的比率'); title('p=0.03试验次数N与正面向上比率的函数图'); 实验结果：

2.不同次数的随机试验均值方差比较 源代码： clear ; clc; close; rand('seed',sum(100*clock)); Titles = ['n=5时' 'n=20时' 'n=25时' 'n=50时' 'n=100时']; Titlestr = cellstr(Titles); X_n_bar=[0]; %the samples of the X_n_bar X_n=[0]; %the samples of X_n N=[5,10,25,50,100]; j=1; num_X_n = 100; num_X_n_bar = 100; h_X_n_bar = figure(1);

matlab代码大全

MATLAB主要命令汇总 MATLAB函数参考 附录1.1 管理用命令 函数名功能描述函数名功能描述 addpath 增加一条搜索路径 rmpath 删除一条搜索路径 demo 运行Matlab演示程序 type 列出.M文件 doc 装入超文本文档 version 显示Matlab的版本号 help 启动联机帮助 what 列出当前目录下的有关文件 lasterr 显示最后一条信息 whatsnew 显示Matlab的新特性 lookfor 搜索关键词的帮助 which 造出函数与文件所在的目录 path 设置或查询Matlab路径 附录1.2管理变量与工作空间用命令 函数名功能描述函数名功能描述 clear 删除存中的变量与函数 pack 整理工作空间存 disp 显示矩阵与文本 save 将工作空间中的变量存盘 length 查询向量的维数 size 查询矩阵的维数 load 从文件中装入数据 who,whos 列出工作空间中的变量名 附录1.3文件与操作系统处理命令 函数名功能描述函数名功能描述 cd 改变当前工作目录 edit 编辑.M文件 delete 删除文件 matlabroot 获得Matlab的安装根目录 diary 将Matlab运行命令存盘 tempdir 获得系统的缓存目录 dir 列出当前目录的容 tempname 获得一个缓存(temp)文件 ! 执行操作系统命令 附录1.4窗口控制命令 函数名功能描述函数名功能描述 echo 显示文件中的Matlab中的命令 more 控制命令窗口的输出页面format 设置输出格式 附录1.5启动与退出命令 函数名功能描述函数名功能描述 matlabrc 启动主程序 quit 退出Matlab环境 startup Matlab自启动程序 附录2 运算符号与特殊字符附录 2.1运算符号与特殊字符 函数名功能描述函数名功能描述 + 加 ... 续行标志 - 减 , 分行符(该行结果不显示) * 矩阵乘 ; 分行符(该行结果显示) .* 向量乘 % 注释标志 ^ 矩阵乘方 ! 操作系统命令提示符 .^ 向量乘方矩阵转置 kron 矩阵kron积 . 向量转置 \ 矩阵左除 = 赋值运算 / 矩阵右除 == 关系运算之相等 .\ 向量左除 ~= 关系运算之不等 ./ 向量右除 < 关系运算之小于

matlab考试题及答案

%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (iA(i)) y=A(i); end end disp('最小元素是：') disp(y) % 7、某系统由四个典型环节构成，试根据四阶龙格—库塔公式，求输出量y的阶跃响应，取仿真时间Tf=10s， %步长h=。

matlab程序大全答案

频率特性类题目 1 一个系统的开环传递函数为 ，试绘制其当K=5、30时系统的开环频率特性Nyquist 图，并判断系统的稳定性。 w=linspace(0.5,5,1000)*pi; sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1) nyquist(sys1,w) title('system nyquist charts with k=5') figure(2) nyquist(sys2,w) title('system nyquist charts with k=30') 由图可知K=5时，开环Nyquist 曲线没有包围（-1，j0）点，所以系统稳定。 K=30时，开环Nyquist 曲线包围（-1，j0）点，所以系统不稳定。 2系统开环传递函数为 ，建立其零极点增益模型， 然后分别绘制当K=5、K=30时系统的开环频率特性Bode 图，并判断系统的稳定性。 sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1) [Gm1,Pm1,Wcg1,Wcp1]=margin(sys1) bode(sys1) title('system bode charts with k=5'),grid figure(2) [Gm2,Pm2,Wcg2,Wcp2]=margin(sys2) bode(sys2) title('system bode charts with k=30'),grid 因为K=5时，Wcg>Wcp,所以系统稳定。 K=10时，Wcg

matlab语音识别系统(源代码)最新版

matlab语音识别系统(源代码)最新版

目录 一、设计任务及要求 (1) 二、语音识别的简单介绍 2.1语者识别的概念 (2) 2.2特征参数的提取 (3) 2.3用矢量量化聚类法生成码本 (3) 2.4VQ的说话人识别 (4) 三、算法程序分析 3.1函数关系 (4) 3.2代码说明 (5) 3.2.1函数mfcc (5) 3.2.2函数disteu (5) 3.2.3函数vqlbg (6) 3.2.4函数test (6) 3.2.5函数testDB (7) 3.2.6 函数train (8) 3.2.7函数melfb (8) 四、演示分析 (9) 五、心得体会 (11) 附：GUI程序代码 (12)

一、设计任务及要求 用MATLAB实现简单的语音识别功能； 具体设计要求如下： 用MATLAB实现简单的数字1~9的语音识别功能。 二、语音识别的简单介绍 基于VQ的说话人识别系统，矢量量化起着双重作用。在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。在识别(匹配)阶段，我们用VQ方法计算平均失真测度(本系统在计算距离d时，采用欧氏距离测度)，从而判断说话人是谁。 语音识别系统结构框图如图1所示。 图1 语音识别系统结构框图 2.1语者识别的概念 语者识别就是根据说话人的语音信号来判别说话人的身份。语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可能。用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘；语音信号的采集方便，系统设备成本低；利用电话网络还可实现远程客户服务等。因此，近几年来，说话人识别越来越多的受到人们的重视。与其他生物识别技术如指纹识别、手形识别等相比较，说话人识别不仅使用方便，而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中，是唯一可以用作远程验证的识别技术。因此，说话人识别的应用前景非常广泛：今天，说话人识别技术已经关系到多学科的研究领域，不同领域中的进步都对说话人识别的发展做出了贡献。说话人识别技术是集声学、语言学、计算机、信息处理和人工智能等诸多领域的一项综合技术，应用需求将十分广阔。在吃力语音信号的时候如何提取信号中关键的成分尤为重要。语音信号的特征参数的好坏直接导致了辨别的准确性。

利用MATLAB进行时域分析

自动控制原理与系统课程实验报告 实验题目：利用MATLAB进行时域分析 班级：机电1131班姓名：刘润学号：38号 一、实验目的及内容 时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。 二、实验设备 三、实验原理 典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律： 已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn）， 其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。

四、实验步骤编出程序如下图： 五、实验结果画出图表如下图：

六、结果分析 （1）当ξ=0（无阻尼）（零阻尼）时： 无阻尼时的阶跃响应为等幅振荡曲线。如图ξ=0曲线。 （2）当0<ξ<1（欠阻尼）时： 对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。如图ξ=0.4曲线。 (3)当ξ=1（临界阻尼）时： 临界阻尼时的阶跃响应为单调上升曲线。如图ξ=1曲线。 （4）当ξ>1（过阻尼）时： 过阻尼时的阶跃响应也为单调上升曲线。不过其上升的斜率较临界阻尼更慢。如图ξ=1.6曲线 七、教师评语

基于MATLAB的潮流计算源程序代码(优.选)

%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;

matlab代码大全教学文案

m a t l a b代码大全

MATLAB主要命令汇总 MATLAB函数参考 附录1.1 管理用命令 函数名功能描述函数名功能描述 addpath 增加一条搜索路径 rmpath 删除一条搜索路径 demo 运行Matlab演示程序 type 列出.M文件 doc 装入超文本文档 version 显示Matlab的版本号 help 启动联机帮助 what 列出当前目录下的有关文件 lasterr 显示最后一条信息 whatsnew 显示Matlab的新特性 lookfor 搜索关键词的帮助 which 造出函数与文件所在的目录 path 设置或查询Matlab路径 附录1.2管理变量与工作空间用命令 函数名功能描述函数名功能描述 clear 删除内存中的变量与函数 pack 整理工作空间内存 disp 显示矩阵与文本 save 将工作空间中的变量存盘 length 查询向量的维数 size 查询矩阵的维数 load 从文件中装入数据 who,whos 列出工作空间中的变量名 附录1.3文件与操作系统处理命令 函数名功能描述函数名功能描述 cd 改变当前工作目录 edit 编辑.M文件 delete 删除文件 matlabroot 获得Matlab的安装根目录 diary 将Matlab运行命令存盘 tempdir 获得系统的缓存目录 dir 列出当前目录的内容 tempname 获得一个缓存(temp)文件 ! 执行操作系统命令 附录1.4窗口控制命令 函数名功能描述函数名功能描述 echo 显示文件中的Matlab中的命令 more 控制命令窗口的输出页面format 设置输出格式 附录1.5启动与退出命令 函数名功能描述函数名功能描述 matlabrc 启动主程序 quit 退出Matlab环境 startup Matlab自启动程序 附录2 运算符号与特殊字符附录

matlab课程设计题目

课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图

形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。

MATLAB简单程序大全

MATLAB简单程序大全 求特征值特征向量 A=[2 3 4;1 5 9;8 5 2] det(A) A' rank(A) inv(A) rref(A) eig(A)%求特征值和特征向量 卫星运行问题 h=200,H=51000,R=6378; a=(h+H+2*R)/2; c=(H-h)/2; b=(a^2-c^2)^(1/2); e=c/a; f=sqrt(1-exp(2).*cos(t)^2); l=int(f,t,0,pi/2) L=4*a.*l 动态玫瑰线 n=3;N=10000; theta=2*pi*(0:N)/N; r=cos(n*theta); x=r.*cos(theta); y=r.*sin(theta); comet(x,y) 二重积分 syms x y f=x^2*sin(y); int(int(f,x,0,1),y,0,pi) ezmesh(f,[0,1,0,pi]) 函数画图 syms x;f=exp(-0.2*x)*sin(0.5*x); ezplot(f,[0,8*pi])

玫瑰线 theta=0:0.01:2*pi; r=cos(3*theta); polar(theta,r,'r') 求x^2+y^2=1和x^2+z^2=1所围成的体积 syms x y z R r=1; Z=sqrt(1-x^2); y0=Z; V=8*int(int(Z,y,0,y0),x,0,1) 求导数及图像 f='1/(5+4*cos(x))'; subplot(1,2,1);ezplot(f) f1=diff(f) subplot(1,2,2);ezplot(f1) 绕x轴旋转 t=(0:20)*pi/10; r=exp(-.2*t).*sin(.5*t); theta=t; x=t'*ones(size(t)); y=r'*cos(theta); z=r'*sin(theta); mesh(x,y,z) colormap([0 0 0]) 某年是否闰年 year=input('input year:='); n1=year/4; n2=year/100; n3=year/400; if n1==fix(n1)&n2~=fix(n2) disp('是闰年') elseif n1==fix(n1)&n3==fix(n3) disp('是闰年') else

matlab的编码大全

附录Matlab源程序 附录A 信息熵 % 函数说明：% % H=entropy(P,r) 为信息熵函数% % P为信源的概率矢量, r为进制数% % H为信息熵% %****************************** % function H=entropy(P,r) if (length(find(P<=0))~=0) error('Not a prob.vector,negative component'); % 判断是否符合概率分布条件end if (abs(sum(P)-1)>10e-10) error('Not a prob.vector,component do not add up to 1'); end H=(sum(-P.*log2(P)))/(log2(r)+eps); 附录B 离散无记忆信道容量的迭代计算 % 信道容量C的迭代算法% % 函数说明：% % [CC,Paa]=ChannelCap(P,k) 为信道容量函数% % 变量说明：% % P：输入的正向转移概率矩阵，k：迭代计算精度% % CC：最佳信道容量，Paa：最佳输入概率矩阵% % Pa：初始输入概率矩阵，Pba：正向转移概率矩阵% % Pb：输出概率矩阵，Pab：反向转移概率矩阵% % C：初始信道容量，r：输入符号数，s：输出符号数% %************************************************** % function [CC,Paa]=ChannelCap(P,k) % 提示错误信息 if (length(find(P<0)) ~=0) error('Not a prob.vector,negative component'); % 判断是否符合概率分布条件end

MATLAB习题参考答案

第1章 MATLAB简介 1、MA TLAB的主要特点有：①语言简洁，编程效率高。②人机界面友善，交互性好。③绘图功能强大，便于数据可视化。④学科众多、领域广泛的MATLAB工具箱。⑤源程序的开放性。 MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。 2、填空题 ⑴命令窗口、命令历史窗口、当前目录窗口 ⑵查阅、保存、编辑 ⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。 3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。 4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。） 6+7+8+9 1+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。） 第2章矩阵与数值数组 1、填空题： ⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。 ⑵全下标、单下标。 2、阅读程序题： （本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果） ⑴ans = 2 3 7 Sa = 10 20 30 A = 1 20 5 30 9 10 4 6 8 10 ⑵ A = 1 3 5 7 9 11 13 15

2 4 6 8 10 12 14 16 ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16 A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 0 ⑶ ans = -1 -4 6 4 ans = 3 0 5 -2 ans = 3 6 9 12 ans = 3 6 9 12 3、A=magic(4)； L=A<10 L = 0 1 1 0 1 0 0 1 1 1 1 0 1 0 0 1 4、省略。 第3章字符串、元胞和构架数组 1. 直接创建法：S1=['Where there is life,' 'there is hope. '] %注意第2行要加入6个空格使其与第一 行字符数相等，否则报错。 S2=str2mat('Where there is life,','there is hope.') S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行 字符的数目 2. 填空题：A（2，3）；A{2，3} 3. 阅读程序题： ⑴subch = ABc12

matlab快捷键大全

1.使用“Tab”键可以在Matlab进行命令输入时补全变量名或者命令名； 2.使用“Ctrl+C”组合键能够强制从运行的或者进入死循环的Matlab程序中退出； 3.使用“上”“下”方向键能够调用Matlab的历史命令； 4.Matlab的变量命名规则注意： a.区分大小写； b.必须以字母开头； c.中间不能有空格、标点符号等； 5.在一条命令或语句的末尾加上一个分号，则Matlab就不会在屏幕上显示这条命令计算的结果； 6.在Matlab中，“...”(3个句点)称为“续行符”，表示同一语句的延续输入； 注意：只有当续行符出现在变量名和运算符之间，才能起到语句延续的作用。出现在一个变量名中间时，是不能实现语句延续的。换句话说，变量名不能分割成两行书写。 7.在进行三角函数运算时，Matlab使用的是“弧度制”而不是“角度制”； 8.选择需要注释的语句块以后，使用“Ctrl+R”可以进行Matlab语句注释；使用“Ctrl+T”可以进行语句块的解注释； 9.Matlab数组的下标是从“1”开始的，这点与C语言不同； 1 0.使用空格和逗号隔开的元素指定的是同一行的元素，使用分号或者回车分开的元素指定的是不同行的元素。 Q1：matlab有没有监视内存的方法? A：用函数whos。 Q2：如何解决matlab7.0命令窗口跳出一大堆java错误...

A：换matlab7的sp2。 Q3：自从安装matlab，一开机就在进程里有matlab。 能不能开机的时候进程就不运行matlab？ A：开始-->控制面板-->管理工具-->服务 把MATLAB Server的属性改成“手动”就行了。 Q4：退出matlab7程序运行的快捷键。 A：ctrl+q Q5：matlab7远程控制是否有限制？ A：不能远程控制,不过你可以先在你的remote机器上打开,然后就可以用了。 Q6：Matlab占用资源太多怎么办？ A：用matlab-nojvm启动（如果不需要图形界面）。 Q7：怎样给matlab添加新的toolbox？ A：在matlab的文件菜单里边添加路径，选set path。 Q8：请问matlab有没有命令可以调出历史输入啊？ A：command history窗口。 Q9：matlab7.0不能在64位的cpu下运行？ A：matlab应该是依赖于自己的虚拟机的 但是好像这个虚拟机是在IA32里面作出来的，所以，应该找个带64位的java虚拟机替换原来的，不过不一定能行or so,记不清了） Q10：matlab有没有注释一段的功能？ A：选中一段代码，ctrl r就是区段注释 选中一段代码，ctrl t取消区段注释

最常用的matlab图像处理的源代码

最常用的一些图像处理Matlab源代 码 #1：数字图像矩阵数据的显示及其傅立叶变换 #2：二维离散余弦变换的图像压缩 #3：采用灰度变换的方法增强图像的对比度 #4：直方图均匀化 #5：模拟图像受高斯白噪声和椒盐噪声的影响 #6：采用二维中值滤波函数medfilt2对受椒盐噪声干扰的图像滤波 #7：采用MATLAB中的函数filter2对受噪声干扰的图像进行均值滤波 #8：图像的自适应魏纳滤波 #9：运用5种不同的梯度增强法进行图像锐化 #10：图像的高通滤波和掩模处理 #11：利用巴特沃斯（Butterworth）低通滤波器对受噪声干扰的图像进行平滑处理 #12：利用巴特沃斯（Butterworth）高通滤波器对受噪声干扰的图像进行平滑处理 1.数字图像矩阵数据的显示及其傅立叶变换 f=zeros(30,30); f(5:24,13:17)=1; imshow(f, 'notruesize'); F=fft2(f,256,256); % 快速傅立叶变换算法只能处矩阵维数为2的幂次，f矩阵不 % 是，通过对f矩阵进行零填充来调整 F2=fftshift(F); % 一般在计算图形函数的傅立叶变换时，坐标原点在 % 函数图形的中心位置处，而计算机在对图像执行傅立叶变换 % 时是以图像的左上角为坐标原点。所以使用函数fftshift进 %行修正，使变换后的直流分量位于图形的中心； figure,imshow(log(abs(F2)),[-1 5],'notruesize');

2 二维离散余弦变换的图像压缩I=imread('cameraman.tif'); % MATLAB自带的图像imshow(I); clear;close all I=imread('cameraman.tif'); imshow(I); I=im2double(I); T=dctmtx(8); B=blkproc(I,[8 8], 'P1*x*P2',T,T'); Mask=[1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; B2=blkproc(B,[8 8],'P1.*x',Mask); % 此处为点乘(.*) I2=blkproc(B2,[8 8], 'P1*x*P2',T',T); figure,imshow(I2); % 重建后的图像 3.采用灰度变换的方法增强图像的对比度I=imread('rice.tif'); imshow(I); figure,imhist(I); J=imadjust(I,[0.15 0.9], [0 1]); figure,imshow(J); figure,imhist(J);

matlab 指令大全

分享 我的分享 当前分享 返回分享首页? 分享 matlab命令，应该很全了！来源：李家叶的日志 matlab命令 一、常用对象操作：除了一般windows窗口的常用功能键外。 1、!dir 可以查看当前工作目录的文件。!dir& 可以在dos状态下查看。 2、who 可以查看当前工作空间变量名，whos 可以查看变量名细节。 3、功能键： 功能键快捷键说明 方向上键Ctrl+P 返回前一行输入 方向下键Ctrl+N 返回下一行输入 方向左键Ctrl+B 光标向后移一个字符 方向右键Ctrl+F 光标向前移一个字符 Ctrl+方向右键Ctrl+R 光标向右移一个字符 Ctrl+方向左键Ctrl+L 光标向左移一个字符 home Ctrl+A 光标移到行首 End Ctrl+E 光标移到行尾 Esc Ctrl+U 清除一行 Del Ctrl+D 清除光标所在的字符 Backspace Ctrl+H 删除光标前一个字符 Ctrl+K 删除到行尾 Ctrl+C 中断正在执行的命令 4、clc可以命令窗口显示的内容，但并不清除工作空间。 二、函数及运算 1、运算符： ＋：加，－：减，*：乘，/：除，\：左除^：幂，‘：复数的共轭转置，（）：制定运算顺序。 2、常用函数表： sin( ) 正弦（变量为弧度） Cot( ) 余切（变量为弧度） sind( ) 正弦（变量为度数） Cotd( ) 余切（变量为度数） asin( ) 反正弦（返回弧度） acot( ) 反余切（返回弧度） Asind( ) 反正弦（返回度数） acotd( ) 反余切（返回度数） cos( ) 余弦（变量为弧度） exp( ) 指数 cosd( ) 余弦（变量为度数）

matlab经典习题及解答

第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 页脚内容1

命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 在MATLAB中有多种获得帮助的途径： （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。 第2章MATLAB矩阵运算基础 页脚内容2

BP神经网络matlab源程序代码

close all clear echo on clc % NEWFF——生成一个新的前向神经网络 % TRAIN——对 BP 神经网络进行训练 % SIM——对 BP 神经网络进行仿真 % 定义训练样本 % P为输入矢量 P=[0.7317 0.6790 0.5710 0.5673 0.5948;0.6790 0.5710 0.5673 0.5948 0.6292; ... 0.5710 0.5673 0.5948 0.6292 0.6488;0.5673 0.5948 0.6292 0.6488 0.6130; ... 0.5948 0.6292 0.6488 0.6130 0.5654; 0.6292 0.6488 0.6130 0.5654 0.5567; ... 0.6488 0.6130 0.5654 0.5567 0.5673;0.6130 0.5654 0.5567 0.5673 0.5976; ... 0.5654 0.5567 0.5673 0.5976 0.6269;0.5567 0.5673 0.5976 0.6269 0.6274; ... 0.5673 0.5976 0.6269 0.6274 0.6301;0.5976 0.6269 0.6274 0.6301 0.5803; ... 0.6269 0.6274 0.6301 0.5803 0.6668;0.6274 0.6301 0.5803 0.6668 0.6896; ... 0.6301 0.5803 0.6668 0.6896 0.7497]; % T为目标矢量 T=[0.6292 0.6488 0.6130 0.5654 0.5567 0.5673 0.5976 ... 0.6269 0.6274 0.6301 0.5803 0.6668 0.6896 0.7497 0.8094]; % Ptest为测试输入矢量 Ptest=[0.5803 0.6668 0.6896 0.7497 0.8094;0.6668 0.6896 0.7497 0.8094 0.8722; ... 0.6896 0.7497 0.8094 0.8722 0.9096]; % Ttest为测试目标矢量 Ttest=[0.8722 0.9096 1.0000]; % 创建一个新的前向神经网络 net=newff(minmax(P'),[12,1],{'logsig','purelin'},'traingdm'); % 设置训练参数 net.trainParam.show = 50; net.trainParam.lr = 0.05; net.trainParam.mc = 0.9; net.trainParam.epochs = 5000; net.trainParam.goal = 0.001; % 调用TRAINGDM算法训练 BP 网络 [net,tr]=train(net,P',T); % 对BP网络进行仿真 A=sim(net,P'); figure; plot((1993:2007),T,'-*',(1993:2007),A,'-o'); title('网络的实际输出和仿真输出结果，*为真实值，o为预测值'); xlabel('年份'); ylabel('客运量'); % 对BP网络进行测试 A1=sim(net,Ptest');