Matlab实验报告
Matlab实验报告
四院四队3班
马鑫201204013002
马也201204014005
冯帆201204014023
吴冠楠20120415002
解旭辉201204015007
赖恪201204015019
第一专题作业1、完成课程GUI设计。
第二专题作业
1、 设矩阵A 和B 满足关系式AB=A+2B ，已知423A 110123?? ?= ? ?-??
求矩阵B 。
A=[4 2 3;1 1 0;-1 2 3];
C=2*eye(3,3);
B=(A-C)\A
B =
3.0000 -8.0000 -6.0000
2.0000 -9.0000 -6.0000
-2.0000 12.0000
9.0000
2、求解线性方程组
123
123
123
2x x3x5 3x x5x5 4x x x9
-+=
?
?
+-=
?
?-+=
?
。
A=[2 -1 3;3 1 -5;4 -1 1];B=[5;5;9]; X=A\B
X =
2
-1
3、 求极限1/ln 0lim (tan )x x x +→。
syms x
f=(tan(x))^(1/log(x));
ls=limit(f,x,0,'right')
ls = exp(1)
4、求积分
0x
?
f=(x^(1/2))/(x+1)^2; s=int(f,x,0,+inf)
s =
pi/2
5、 求一阶微分方程的通解y'tan cos y x x +=。
dsolve('Dy+y*tan(x)=cos(x)')
ans = (cos(x) - C2*exp(-t*tan(x)))/tan(x)
第三专题作业
1、 曲线绘制
画出衰减振荡曲线t/2y sin 2t e -=?及其包络线t/2t/201y e ,e y -==（虚线）,t 的取值范围为
[0,6]π,加上适当标注。
2、 曲面绘制
画出22z /(y )x =+所表示的三维曲面，用暖色调，x 、y 的取值范围为[-8,8]。
3、特殊图形绘制
为作业1中衰减振荡曲线绘制直方图，直方图数目为20个。
4、街头骗局揭秘-摸球游戏
一袋装有16个大小、形状相同的球，其中8个红色、8个白色。游戏者从中一次摸出8个，
论值和模拟值。（nchoosek、randsample函数）。
function [theory,simulation]=h34(x)
%计算理论值
p=nchoosek(16,8);
p1=2*nchoosek(8,0)*nchoosek(8,8)/p;
p2=2*nchoosek(8,1)*nchoosek(8,7)/p;
p3=2*nchoosek(8,2)*nchoosek(8,6)/p;
p4=2*nchoosek(8,3)*nchoosek(8,5)/p;
p5=nchoosek(8,4)*nchoosek(8,4)/p;
pp=[p1 p2 p3 p4 p5];
c=[10 1 0.5 0.2 -3];
theory=pp*c';
%模拟摸球
c=0;
b=0;
for i=1:x
R= randsample('01',8,true,[0.5 0.5]);
a=sum(R)-384;
switch a
case 0
b=10;
case 1
b=1;
case 2
b=0.5;
case 3
b=0.2;
case 4
b=-3;
case 5
b=0.2;
case 6
b=0.5;
case 7
b=1;
case 8
b=10;
end
c=c+b;
end
simulation=2*c/x;
%命令窗口执行x=input('输入摸球次数
x=:');[theory,simulation]=h34(x);disp(theory);disp(simulation);
第四专题作业
1、 假定两个输入信号分别是12(t)2cos2100t (t)sin 2500t f f ππ==仿真
12y()(t)(t)t f f =?在仿真速率3000Hz 时的输出波形，
并存入工作空间变量sinprod 中。仿真时间1s 。
2、产生一个均值为2、方差为3的高斯信号，用统计模块测试该信号的直流分量、交流功
率、信号中值，并画出分布的归一化直方图。设信号采样率为1000Hz，仿真时间为10s。
注：histogram，mean，var模块选择running模式直方图输出到m文件中用bar函数绘制，统计分段为100，信号中值用buffer接median模块实现。
第五专题作业
1、同第一专题所列。
时域程序
ts=0:0.001:1;
x=sin(2*pi*50.*ts)+sin(2*pi*100.*ts)+sin(2*pi*150.*ts); subplot(2,1,1)
plot(ts,x);
频域程序
f=-500:1:500;
subplot(2,1,2)
plot(f,abs(fftshift(fft(x))));
2、 根据线性调频信号表达式
22
t t s()cos jsin 22t αα=- 画出脉冲宽度为100微秒、调频频带为1MHz 的线性调频信号功率谱密度。
说明：2B /T απ=为线性调频斜率，B 为调频带宽，T 为信号持续时间。
t=0:0.01:100;
T=100e-006;
B=1e006;
a=2*pi*B/T;
s=cos(a*t.^2/2)-i*sin(a*t.^2/2);
windows=boxcar(length(s));
periodogram(s,windows,512,1000)
3、信号的持续时间为10，且该信号是频率分别为单位功率47Hz、219Hz正弦信号及白噪声
之和X()cos(247t)cos(2219t)()
t n t
ππ
=?+?+其中0≤t ≤10，以每秒1000Hz速率
对X(t)进行抽样，求功率谱密度。
数据在matlab/toolbox/sl3d/vrealm/program/sounds
[y,fs,bits]=wavread('bird');%读取wav格式的文件sound(y,fs,bits);%将数列的数据通过声卡转化为声音Y=fft(y,4096);%fft函数的一种输入输出形式subplot(211);
plot(y);
title('声音信号的波形');
subplot(212)
plot(abs(Y));
title('声音信号的频谱')
4、读入火车声音数据文件train.mat，采样速率8192，分析频谱load(‘train’)。
数据在matlab/toolbox/sl3d/vrealm/program/sounds