matlab例题及答案

1.编写程序：计算1/3+2/5+3/7+……+10/21

法一： s=0;

for i=1:10

s=s+i/(2*i+1); end ss =

4.4096 法二：

sum((1:10)./(3:2:21)) ans =

4.4096

2.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。

s=0;

for i=1:100

if mod(i,3)==0&&mod(i,7)==0 s=s+i; end,end ss =

210

3.画出y=n！的图（1<=n<=10）,阶乘的函数自己编写，禁用MATLAB自带的阶乘函数。

x=1:10; for i=1:10

try y(i)=y(i-1)*i; catch y(i)=1; end,end plot(x,y)

43.532.521.510.50x 10612345678910

4.一个数恰好等于它的因子之和，这个数就称为完数。例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。编程找出2000以内的所有完数。

g=[];

for n=2:2000 s=0;

for r=1:n-1

if mod(n,r)==0

字符串举例：

>> A= 'China ' '中国' ' ' ；输出带引号的汉字 >> B='我是好学生' >> c='I am fine.'

>> s3=char('s', 'y','m','b','o','l','i','c') ；用函数char生成字符串 >> double(s3') ；字符串转换为数值代码 >> abs(s3')

>> cellstr(s3) ；字符矩阵转换为字符串 >> b=num2str(a) 比较 >> b*2 和 >> ab=[A,' ',B,'.'] >> AB=['中国';'北京']

关系运算举例

>> A=[3,4,8;9,0,2;5,3,7] A = 3 4 8 9 0 2 5 3 7 >> B=[4,4,1;7,8,4;5,1,7] B =

4 4 1 7 8 4 5 1 7 >> E=(A>B) E = 0 0 1

《现代信号处理》期末考核作业

1 MATLAB仿真均值为0，方差为1的白噪声信号，信号长度N=1024，并用周期图法分别求500、1000和1500次实现的平均功率谱密度，画图。

程序代码如下：

clear; clear all; N=1024;%数据长度

Nfft=1024;?T所采用的数据长度 n=0:N-1;

wn=randn(1,N);%产生随机白噪声 subplot(2,2,1);%绘出白噪声序列 plot(n,wn); title('白噪声');

P0次实现的平均功率谱密度 s=zeros(1,N); for i=1:500

wn=randn(1,N);%产生随机白噪声

Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为db s=s+Pxx; end s=s/500;

f=(0:length(Pxx)-1)/length(Pxx);%绘出频率序列 subplot(222); plot(f,s);

xlabel('频率/Hz');ylabel('功率谱/dB'); title('500次实现的平均功率谱密度'); grid on;

00次实现的平均功率谱密度 s=

《数学实验》在线习题3

Matlab程序设计部分 一. 分

析

向

量

组

a1?[1T2a23?]?,?T[a31T?2，0],a4?[1?2?1]T，a5?[246]T的线性相关性，找出它们的最大无关组，并将其

余向理表示成最大无关组的线性组合。

解， a1=[1 2 3]';

a2=[-1 -2 0]'; a3=[0 0 1]'; a4=[1 -2 -1]'; a5=[2 4 6]'; A=[a1,a2,a3,a4,a5] ; [R,S]=rref(A) r=length(S)

R =

1.0000 0 0.3333 0 2.0000 0 1.0000 0.3333 0 0 0 0 0 1.0000 0 S =

1 2 4 r =

3

线性相关 a1,a2,a3,a4,a5 最大无关组是a1,a2,a4

其余向量的线性组合是a3=1/3a1+1/3a2 a5=2a1

二. 计算行列式

第1章

例1- 离散时间信号产生。编写程序产生下列基本脉冲

（1）单位脉冲序列：起点ns，终点ne，在n0处有一单位脉冲（ns≤n0≤ne）。

（2）单位阶跃序列：起点ns，终点ne，在n0前为0，在n0处及以后为1（ns≤n0≤ne）。 （3）实数指数序列：x3?(0.9)n （4）正弦序列：x4?sin(n)

程序如下：

ns=0;ne=15;n0=5;

n1=ns:ne;x1=[(n1-n0)==0]; %n2=ns:ne;x2=[(n2-n0)>=0]; %n3=ns:ne;x3=(0.9).^n3; %n4=ns:ne;x4=sin(n4); %subplot(2,2,1),stem(n1,x1); title('单位脉冲序列')

subplot(2,2,2),stem(n2,x2); title('单位阶跃序列')

subplot(2,2,3),stem(n3,x3); title('实数指数序列')

subplot(2,2,4),stem(n4,x4); title('正弦序列')

运行结果如图1-所示。

《数学实验》在线习题3

Matlab程序设计部分 一. 分

析

向

量

组

a1?[1T2a23?]?,?T[a31T?2，0],a4?[1?2?1]T，a5?[246]T的线性相关性，找出它们的最大无关组，并将其

余向理表示成最大无关组的线性组合。

解， a1=[1 2 3]';

a2=[-1 -2 0]'; a3=[0 0 1]'; a4=[1 -2 -1]'; a5=[2 4 6]'; A=[a1,a2,a3,a4,a5] ; [R,S]=rref(A) r=length(S)

R =

1.0000 0 0.3333 0 2.0000 0 1.0000 0.3333 0 0 0 0 0 1.0000 0 S =

1 2 4 r =

3

线性相关 a1,a2,a3,a4,a5 最大无关组是a1,a2,a4

其余向量的线性组合是a3=1/3a1+1/3a2 a5=2a1

二. 计算行列式

实验一、MATLAB基本操作

一、实验目的

2．学习使用图形函数计算器命令funtool及其环境。 3. 学习使用help命令进行帮助

4. 掌握向量与矩阵的创建以及矩阵的基本操作 5. 掌握数组与矩阵的概念 二、 实验内容

熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；掌握数组与矩阵的概念；学会使用help命令进行帮助；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool；

1. 命令窗 口的简单使用 （1） 简单矩阵的输入 （自由创建）

x=[1 3 5;2 4 6]

x =

1 3 5

2 4 6

（2） 求[12+2×(7-4)]÷32的算术运算结果，总结算术运算符先级

[12+2*(7-4)]/3^2

ans =

2

2. 有关向量、矩阵或数组的一些运算 （1） 设A=15；B=20；求C=A+B与c=a+b的区别 A=15;B=20; >> C=A+B C = 35 >> c=a+b

??? Undefined f

一元函数微分学

实验1 一元函数的图形（基础实验）

实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab作平面曲线图性的方法与技巧. 初等函数的图形

2 作出函数y?tanx和y?cotx的图形观察其周期性和变化趋势. 解：程序代码： >> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps);

plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象：

tan(x)50403020100-10-20-30-40-500123456

程序代码：

>> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps);

plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象：

cot(x)50403020100-10-20-30-40-500123456

8

4在区间[?1,1]画出函数y?sin1x的图形. 解：程序代码：

>> x=

一元函数微分学

实验1 一元函数的图形（基础实验）

实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab作平面曲线图性的方法与技巧.

初等函数的图形

2 作出函数y?tanx和y?cotx的图形观察其周期性和变化趋势. 解：程序代码：

>> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps);

plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象：

tan(x)50403020100-10-20-30-40-500123456程序代码：

>> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps);

plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象：

cot(x)50403020100-10-20-30-40-500123456

8

4在区间[?1,1]画出函数y?sin1x的图形. 解：程序代码：

>> x=lin

含两个递进例题,有详细解释,以及matlab程序

设高度的测量误差是均值为0、方差为1的高斯白噪声随机序列，该物体的初始高度h0和速度V0也是高斯分布的随机变量，且

Eh0 1900m h0 1000

。试求该物体高度 ,P var 0

EV0 10m/s V0 02

和速度随时间变化的最优估计。（g 9.80m/s2） 解：

1. 令X(k)

h(k)

t=1 R（k）=1 Q(k)=0 v(k)

根据离散时间卡尔曼滤波公式，则有： X(k 1) (k 1,k)X(k) U(k) Y(k 1) H(k 1)X(k 1) V(k 1)

0.5gt2 1 t

k＝) (k 1, U(k)＝ ＝ 10 H(k 1) 1 gt

1900

滤波初值：X(0|0) EX(0)

10

^

P(0|0) var[X(0)] P0

100

2

一步预测：X(k 1|k) (k 1,k)X(k|k) U(k) P(k 1|k) (k 1,k)P(k|k) (k 1,k)

滤波增益：K(k 1) P(k 1|