《科技学院matlab实验》由会员分享，可在线阅读，更多相关《科技学院matlab实验（12页珍藏版）》请在金锄头文库上搜索。
1、华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 0 页实验实验 环境环境MATLAB 7.5实验实验 名称名称实验一：熟悉 MATLAB 环境及 Matlab 基本操作实实 验验 目目 的的1、熟悉 MATLAB 环境及各种基本 Matlab 操作实实 验验 原原 理理运算符= 赋值 ： 用来输入行向量，例如： x =1:2:20 %用户自定义增加量为 2； 分隔行 控制是否显示结果 例如： A 显示矩阵 A 的结果 %显示矩阵 A 的结果 A； 不显示矩阵 A 的结果 %不显示矩阵 A 的结果 变量 ans 自动存储未指定输出变量的表达式结果Pi 圆周率 语句 1）表达式语句
2、A=3+4+5B=abC=xor(a,b)D=sin(x)2)条件分支语句 ifend ifelseendwhileend3)转移语句 break Return 图形处理 figure %创建一个新的空图形窗口 plot（x,y） %以 x 为横坐标，y 为纵坐标绘制二维连续图形 subplot(m,n,p)%把当前窗口分成 m*n 个子窗口，选择低 p 个为当前子窗口 图形修饰 plot（x,y,color mark）%color 代表颜色，如r代表红色， g代表绿色；mark 代表线型， +代表用加号划线 title(sting)%给图形添加标题 xlabel(sting)%给 x 轴添加
3、说明 ylabel(sting)% 给 y 轴添加说明华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 1 页实实 验验 内内 容容1、将 Matlab 当前目录变更到 D:自己姓名拼音文件夹下,并建立新文件画出 y=2 * x 的曲线， 并将曲线图形保存在当前目录下。源程序：x=1:1:10;y=2*xplot(x,y,g)title(y=2x)2、使用 for 语句及 while 语句求解 1 至 100 整数中奇数的和。源程序： s=0 for i=1:2:100;s=s+i;endsum=s;sum3、用 MATLAB 程序建立九九乘法表。 源程序：for x=1:9; f
4、or y=1:x;A(x,y)=x*y; end end A=4、求连续自然数的和，当和大于等于 1000 时，求最后一个自然数以及自然数的和。 源程序：x=0; i=1; while(i) x=x+i; i=i+1; if(x=1000)y=i-1;i=0; end end j=y;j sum=x;sum5、MATLAB函数是可以递归调用的，即在函数的内部调用函数自身，试用递归方式编写n!函 数。 源程序：function y=fun(n) n=input(n=); y=1; for k=1:n y=y*k; end华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 2 页实实 验
5、验 结结 果果 及及 分分 析析1、当X=110，Y= 2、4、6 、8 、10 、12 、14 、16 、18 、20123456789102468101214161820y=2x2、s =0 sum =25003、A = 1 0 0 0 0 0 0 0 02 4 0 0 0 0 0 0 03 6 9 0 0 0 0 0 04 8 12 16 0 0 0 0 05 10 15 20 25 0 0 0 06 12 18 24 30 36 0 0 07 14 21 28 35 42 49 0 08 16 24 32 40 48 56 64 09 18 27 36 45 54 63 72 814、
6、j =45 sum =10355、在n=后边输入数字即可得结果，例如：输入4得结果24华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 3 页实验实验 名称名称实验二： Matlab 在通信原理中的应用实验实验 目的目的1、掌握通信原理课程中基本知识的 Matlab 编程、仿真方法实实 验验 原原 理理具有离散大载波的双边带调幅（AM） 设模拟基带信号为 M（t）,调幅新号为 s(t)=A+m(t)cos2ft 其中 A 是一个常数，可以将调幅新号看成一个余弦载波加抑制载波双边带调幅 信号，当 Am(t)时，称此调幅信号欠调幅；当 A 1、用 Matlab 产生一个频率为 2Hz、
7、功率为 1 的正弦信源 m(t),设载波频率为 10Hz，A 分别 等于 0.5、1、1.5，试画出： （1）AM 调制信号； （2）调制信号的功率谱密度； （3）相干解调后的信号波形 （4）研究调制深度与失真的关系。 源程序： clear all; dt=0.01; fm=2; fc=10; t=0:dt:5; mt=sqrt(2)*sin(2*pi*fm*t); A=2; s_am=(A+mt).*cos(2*pi*fc*t); figure(1) plot(t,s_am);hold on; plot(t,A+mt,r-); title(AM调制信号及其包络); xlabel(t); t1
8、=0:0.01:2; mt=sqrt(2)*sin(2*pi*fm*t1); z=(A+mt).*cos(2*pi*fc*t1) y=fft(z,1024); n=0:1023;f=100*n/1024; figure(2) plot(f,abs(y);title(FFT);华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 4 页实实 验验 结结 果果 及及 分分 析析00.511.522.533.544.55-4-3-2-101234AM以 以 以 以 以 以 以 以t0102030405060708090100020406080100120140160180200FFT华 北
9、电 力 大 学 (科 技 学 院) 实 验 报 告第 5 页实验实验 名称名称实验三 Matlab 在信号与系统中的应用实实 验验 目目 的的1、掌握信号与系统课程中基本知识的 Matlab 编程、仿真方法实实 验验 原原 理理对于连续时间信号 f(t)，可以用两个行向量 f 和 t 来表示，向量 t 是形 如 t=t1:p:t2 的命令，定义时间范围，t1 为信号起始时间，t2 为终止时间， p 为时间间隔。F 为连续时间信号 f(t)在 t 向量所定义时间点上的样值。对 于连续信号 f=sin(t)，可以用如下两个行向量 f 和 t 表示： t=0:0.1:10; f=sin(t) MAT
10、LAB提供了求连续系统冲激响应及阶跃响应，并绘制其时域波形的函数 impulse和step。描述连续系统的微分方程时，用向量a、b表示该系统，以微分 方程时间求导的降幂次序来排列，且缺项填0。例如：2y(t)+3y(t)+6y(t)=f(t)+f(t)表示该系统的对应向量为：a=2 3 6;b=1 0 1; 表示冲激响应：impulse(b,a) impulse(b,a,t) impulse(b,a,t1:p:t2) :画时域波形 y=impulse(b,a,t1:p:t2)：求出数值解 表示阶跃响应：step (b,a) step(b,a,t) step(b,a,t1:p:t2) :画时域波
11、形 y=step(b,a,t1:p:t2)：求出数值解华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 6 页实实 验验 内内 容容1、用 MATLAB 在时域中编程求解 y(t)+y(t)=f(t), f(t)= exp(-3t)(t)的冲激响应、 阶跃响应及零状态响应。在 simulink 仿真环境下，f(t)= (t),设计系统框图，分析系 统的冲激响应、阶跃响应。 clear all; b=1; a=1 1; p=0.5; %定义取样时间间隔 t=0:p:5; %定义时间范围 impulse(b,a); impulse(b,a,t); step(b,a,t); q=exp(
12、-3*t); %定义输入信号 lsim(b,a,q,t); %对系统输入信号进行仿真 Simulink仿真：1s+1 Transfer FcnStepScopedu/dtDerivative华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 7 页实实 验验 结结 果果 及及 分分 析析0500.10.20.30.40.50.60.70.80.91Impulse ResponseTime (sec)Amplitude02400.10.20.30.40.50.60.70.80.91Step ResponseTime (sec)Amplitude02400.10.20.30.40.50.
13、60.70.80.91Linear Simulation ResultsTime (sec)Amplitude函数 lism（）对系统响应进行仿真的效果取决于向量 t 的时间间隔很多密集 程度。华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 8 页实验实验 名称名称实验四 Matlab 在数字信号处理中的应用实实 验验 目目 的的1、掌握数字信号处理课程中基本知识的 Matlab 编程、仿真方法实实 验验 原原 理理在数字信号处理中，若要表示具有特定采样频率色信号，需定义时间轴 向量。例如：以采样频率 fs=1KHZ 对信号 y(t)=sin(2f1t) +2sin(2f2t)进
14、 行采样，其中，f1=50Hz，f2=120Hz,用 MATLAB 实现如下：f1=50；f2=120；fs=1000t=0:1/fs:1;y=sin(2*pi*f1*t)+ 2*sin(2*pi*f2*t) ;plot(t,y) 傅里叶变换：用 y=fft(x,N)来计算 x 的 N 点 FFT。当 x 的长度大于 N 时，截断 x ;否则补零。 计算模值（或绝对值）的函数：abs（x）当 x 为实数时，计算 x 的绝对 值；当 x 为复数时，计算 x 的模值。实实 验验 内内 容容1、对于连续信号 x(t)=1+cos(2ft)，其中 f=5kHz,分别以采样频率 fs=6 kHz 和 f
15、s=12kHz 对 其进行采样， （1）分别绘出对应的采样信号。 （2）对信号进行傅里叶变换，绘出对应的曲线。 源程序： clear all; f=5000; fs1=6000; fs2=12000; t1=0:1/fs1:0.01; t2=0:1/fs2:0.01; y1=1+2*cos(2*pi*f*t1); y2=1+2*cos(2*pi*f*t2); subplot(2,2,1); plot(t1,y1); title(以6KHz采样); xlabel(t1); subplot(2,2,2); plot(t2,y2); title(以12KHz采样); xlabel(t2);华 北 电 力 大 学 (科 技 学 院) 实 验 报 告第 9 页
