系统辨识第一次作业.zip

系统辨识课课后作业的代码，最小二乘法辨识，包括运行结果。

function [ba]=Hankel（ngkT0） 
%用Hankel矩阵法求系统传函；
%输入参数：n：系统阶数；gk：系统脉冲响应序列；T0：采样间隔
%返回参数：b: 系统脉冲传函分子系数；a:系统脉冲传函分母系数

H=[];
for i=1:n                   %构造Hankel矩阵H
    for j=1:n
        temp（j）=gk（i+j）;        
    end
    H=[H;temp];
end

g1=gk（n+2:2*n+1）*（-1）;
a=flipud（H\g1）;            %获得脉冲传函分母系数a

g2=gk（2:n+1）;
for i=1:n                 %根据a构造矩阵A
    for j=1:n
        if i            A（ij）=0;
        else if i==j
                A（ij）=1;
            else A（ij）=a（i-j）;
            end
        end
    end
end
b=A*g2;                  %获得脉冲传函分子系数b
%调整系数
b=b.‘;                  
a=[1 a.‘];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-05-09 00:38  系统辨识第一次作业\1st\
     文件         740  2009-09-17 20:40  系统辨识第一次作业\1st\Hankel.m
     文件         684  2019-05-09 00:30  系统辨识第一次作业\1st\HankelTestMain.m
     文件       40376  2019-05-09 00:38  系统辨识第一次作业\1st\impulse_si.mdl
     文件       47635  2019-05-09 00:14  系统辨识第一次作业\1st\ori_c&d_impulse_response.jpg
     文件       12147  2019-05-08 21:09  系统辨识第一次作业\1st\param.mat
     文件         328  2019-05-08 23:21  系统辨识第一次作业\1st\test.m
     文件       18378  2019-05-08 19:43  系统辨识第一次作业\1st\untitled.slx
     文件       59592  2019-05-09 00:14  系统辨识第一次作业\1st\辨识后系统脉冲响应对比.jpg
     文件       56660  2019-05-09 00:20  系统辨识第一次作业\1st\辨识后系统脉冲响应对比0.01.jpg
     文件       57052  2019-05-09 00:23  系统辨识第一次作业\1st\辨识后系统脉冲响应对比0.1.jpg
     文件       60099  2019-05-09 00:25  系统辨识第一次作业\1st\辨识后系统脉冲响应对比0.23.jpg
     文件       59096  2019-05-09 00:26  系统辨识第一次作业\1st\辨识后系统脉冲响应对比t1.3.jpg
     文件      677672  2019-05-09 00:37  系统辨识第一次作业\脉冲响应法系统辨识.pptx