最全系统辨识源代码，包括多种最小二乘法，极大似然法，模型阶数辨识等

本人系统辨识课程的全部代码 以及报告 报告里有所有算法原理。内容如下： 第一章 最小二乘法 1 1.1 问题重述 1 1.2 最小二乘法 1 1.2.1 基本最小二乘法 1 1.2.2 不需矩阵求逆的最小二乘法 2 1.2.3 递推最小二乘法 4 1.3 辅助变量法 6 1.3.1 一次辅助变量法 6 1.3.2 递推辅助变量法 7 1.4 广义最小二乘法 9 1.4.1 一次广义最小二乘法 9 1.4.2 递推广义最小二乘法 10 1.5 夏式法 12 1.5.1 夏式偏差修正法 12 1.5.2 夏式改良法 13 1.5.3 递推夏式法 13 1.6 增广矩阵法 16 1.7 自编方法-多阶段最小二乘法 18 1.8 噪声特性分析 19 1.8.1 时域波形 20 1.8.2 均值分析 20 1.8.3 方差分析 21 1.8.4 自相关函数分析 21 1.8.5 功率谱密度分析 22 1.8.6 总结 22 第二章 极大似然法 23 第三章 方法比较 25 3.1 问题重述 25 3.2 各方法精度对比 25 3.3 各方法计算量对比 25 3.4 噪声方差的影响 26 3.5 白噪声和有色噪声对辨识的影响 27 第四章 系统模型阶次的辨识 28 4.1 问题重述 28 4.2 按残差方差定阶 28 4.2.1 按估计误差方差最小定阶 28 4.2.2 F检验法 29 4.3 按AKAIKE信息准则定阶 29 4.4 按残差白色定阶 30 4.5 噪声对定阶的影响 31 4.6 三种方法的优劣及有效性 31 附录 32

% 噪声特性分析
clear;clc;
load（‘uy3.mat‘）;
theta=[1.5;0.8;0.5;0.2];        %参数真值
n=2;                            %系统阶数
N=length（u）-n;                  %矩阵维数
Y=y（n+1:n+N1）;                 %输出矩阵
Phi=zeros（N2*n）;
for i=1:n
    Phi（:i）=-y（n-i+1:n+N-i1）;
    Phi（:n+i）=u（n-i+1:n+N-i1）;
end
Y_=Phi*theta;
e=Y_-Y;                         %残差

%波形
figure（1）
plot（1:Ne）


%均值分析
e_mean=mean（e）

%方差分析
e_var=var（e）


%自相关函数分析
C=xcorr（e‘unbiased‘）;
NC=（length（C）+1）/2;
j=1;
for i=-（NC-1）:（NC-1）
    C（j）=C（j）/abs（i）;
    j=j+1;
end
figure（2）
plot（-（NC-1）:（NC-1）C）


%功率谱密度分析
nfft=1024;
cxn=xcorr（e‘unbiased‘）;
CXk=fft（cxnnfft）;
Pxx=abs（CXk）;
figure（3）
plot（1:1024Pxx）
axis（[11024010]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-04-05 09:05  code\
     文件         776  2015-12-16 23:26  code\Analysis.m
     文件        1091  2015-12-16 16:04  code\GLS.m
     文件        1071  2015-12-16 17:57  code\Hsia1.m
     文件         973  2015-12-16 18:03  code\Hsia2.m
     文件         452  2015-12-16 03:33  code\IV.m
     文件         346  2015-12-15 22:31  code\LS.m
     文件        1001  2015-12-16 02:36  code\LS_without_inv.m
     文件        1495  2015-12-17 02:21  code\ML.m
     文件         737  2015-12-17 03:26  code\MSLS.m
     文件         702  2015-12-16 21:44  code\MSLS2.m
     文件        1842  2015-12-16 18:46  code\RAM.m
     文件        2068  2015-12-16 16:35  code\RGLS.m
     文件        2015  2015-12-16 18:22  code\RHsia.m
     文件         896  2015-12-16 03:58  code\RIV.m
     文件         608  2015-12-15 22:58  code\RLS.m
     文件         746  2015-12-17 05:06  code\SD1.m
     文件        1575  2015-12-17 05:15  code\SD2.m
     文件        2278  2015-12-17 05:28  code\SD3.m
     文件         986  2015-12-14 17:47  code\TS.m
     文件       13147  2015-11-11 23:39  code\uy1.mat
     文件       13372  2015-11-11 23:40  code\uy2.mat
     文件       13126  2015-11-11 23:40  code\uy3.mat
     文件     1232451  2015-12-18 11:21  系统辨识大作业 - 殷泽阳.pdf