matlab一直系统函数画脉冲响应,Matlab 相关分析法求系统脉冲响应（三）

第三部分：相关改进

基本采用(二)中的做法，只是在y的求解摒弃了原理中介绍的方法，直接用matlab自带的lsim函数求解，使本例更有通用性。

源代码：

function mytemp

clear all

close all

clc

%创建M序列

Np=63;%循环周期

delta_T = 1;%时钟节拍

a=1;%幅度

%初始化M序列

M(1)=1;

M(2)=0;

M(3)=0;

M(4)=1;

M(5)=1;

M(6)=0;

M_XuLie(Np) = 0;

for n = 1 : Np

temp = xor(M(6), M(5));

if(temp == 0)

M_XuLie(n) = a;

else

M_XuLie(n) = -a;

end

M(6) =

M(5);

M(5) =

M(4);

M(4) =

M(3);

M(3) =

M(2);

M(2) =

M(1);

M(1) =

temp;

end

figure(4)

stairs(M_XuLie)

ylim([-2 2])

%生成M序列完毕

r = 3; % 周期数

u=repmat(M_XuLie,1,r+1);%将M序列赋给输入，作为输入信号

%第一步，从u(k)得到x(k),y(k)

K = 120;

T0 = 1; %

采样时间

T1 =