实验二 基于MATLAB的序列线性卷积实验

实验二基于MATLAB的序列线性卷积实验

一、实验目的：

1.学会使用MATLAB；

2.通过实验了解MATLAB如何进行数字信号的线性卷积、相关函数的计算；

3.掌握MATLAB的编程方法；

二、实验内容：

1.用MATLAB编程计算两个序列的线性卷积，画出其波形。

代码：

clear;

N = 5;

M = 6;

L = N+M-1;

x = [1,2,3,4,5];

h = [6,2,3,6,4,2];

y = conv(x,h);

nx = 0:N-1;

nh = 0:M-1;

ny = 0:L-1;

subplot(231);

stem(nx,x,'.k');xlabel('n');ylabel('x(n)');grid on;

subplot(232);

stem(nh,h,'.k');xlabel('n');ylabel('h(n)');grid on;

subplot(233);

stem(ny,y,'.k');xlabel('n');ylabel('y(n)');grid on;

运行结果：

2.求两个序列的互相关函数，或一个序列的自相关函数，画出其波形。

代码：

clear;

N = 500;

p1 = 1;

p2 = 0.1;

f = 1/8;

Mlag = 50;

u = randn(1,N);

n = [0:N-1];

s = sin(2*pi*f*n);

%混有告示白噪音的正弦型号的自相关

u1 = u*sqrt(p1);

x1 = u1(1:N)+s;

rx1 = xcorr(x1,Mlag,'biased');

subplot(221);

plot(x1(1:Mlag));

xlabel('n');

ylabel('x1(n)');grid on;

subplot(223);

plot((-Mlag:Mlag),rx1);grid on;

xlabel('m');ylabel('rx1(m)');

%高斯白噪声功率由原来的p1减少为p2,再观察混合型号的自相关

u2 = u*sqrt(p2);

x2 = u2(1:N)+s;

rx2 = xcorr(x2,Mlag,'biased');

subplot(222);

plot(x2(1:Mlag));

xlabel('n');ylabel('x2(n');grid on;

subplot(224);

plot((-Mlag:Mlag),rx2);

grid on;

xlabel('m');

ylabel('rx2(m)');

运行结果：