西南交通大学研究生一年级现代信号处理作业

西南交通大学研究生一年级现代信号处理作业,代码均为本人所写，保证可以运行

%pq=3peqe=5M=30;
N = 500;%样本点数
B = 1;%分母项，确定为AR模型
A =[1-1.30.86 -0.676];%分子项，可根据它求出实际的功率谱
yangben=filter（B A0.1*randn（1N））;%滤波产生样本

r=zeros（1N）;%为总体最小二乘演化公式z，详见式3.4.45
pp=0;%有效秩初始化
for p=1:35
    for i=1:N-p
        r（p）=r（p）+yangben（i）*yangben（i+p）;%混入噪声时的自相关函数
    end
    r（p）=1/（N-p）*r（p）;
end
for l=1:30
    re（l:）=[r（l+5） r（l+4） r（l+3） r（l+2） r（l+1） r（l）];%生成矩阵RE即为3.4.43
end
[USV] =svd（re）;%奇异值分解
for i=1:6
    s=S/S（11）;%归一化奇异值
end
for i=1:6
    if s（ii）>0.05
        pp=pp+1;%确定有效秩
    end
end
pp
ss1=zeros（pp+1）;
for tt=1:6-pp
    for j=1:pp
        ss1=ss1+S（jj）^2*V（tt:tt+ppj）*V（tt:tt+ppj）‘;%式3.4.56，求出矩阵SP
    end
end
a=inv（ss1）;
b=a（1:pp+11）/a（11）%式3.4.59，求出未知参数的最小二乘估计
% bfi=b/b（1）
[h1 w1] = freqz（B A 100）;
[h2 w2] = freqz（1 b‘ 100）;
plot（w1 abs（h1） ‘r‘  w2 abs（h2））;%画图，比较最小二乘估计与实际模型的差别。
legend（‘估计值‘‘实际值‘2）
% zplane（BA）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1052  2018-01-04 21:46  信号处理代码\arma.m

     文件        845  2018-01-04 21:46  信号处理代码\bandp.m

     文件       2041  2018-01-04 21:46  信号处理代码\fmlin.m

     文件        283  2018-01-04 21:46  信号处理代码\kalman_filter.m

     文件       2409  2018-01-04 21:47  信号处理代码\mfiles\altes.m

     文件       2884  2018-01-04 21:47  信号处理代码\mfiles\ambifunb.m

     文件       5455  2018-01-04 21:47  信号处理代码\mfiles\ambifuwb.m

     文件       1620  2018-01-04 21:47  信号处理代码\mfiles\amexpo1s.m

     文件       1628  2018-01-04 21:47  信号处理代码\mfiles\amexpo2s.m

     文件       1654  2018-01-04 21:47  信号处理代码\mfiles\amgauss.m

     文件       1601  2018-01-04 21:47  信号处理代码\mfiles\amrect.m

     文件       1634  2018-01-04 21:47  信号处理代码\mfiles\amtriang.m

     文件       1971  2018-01-04 21:47  信号处理代码\mfiles\anaask.m

     文件       2060  2018-01-04 21:47  信号处理代码\mfiles\anabpsk.m

     文件       2100  2018-01-04 21:48  信号处理代码\mfiles\anafsk.m

     文件       1397  2018-01-04 21:48  信号处理代码\mfiles\anapulse.m

     文件       2217  2018-01-04 21:48  信号处理代码\mfiles\anaqpsk.m

     文件       1878  2018-01-04 21:48  信号处理代码\mfiles\anasing.m

     文件       1390  2018-01-04 21:48  信号处理代码\mfiles\anastep.m

     文件       5976  2018-01-04 21:48  信号处理代码\mfiles\atoms.m

     文件       8126  2018-01-04 21:48  信号处理代码\mfiles\Contents.m

     文件       3805  2018-01-04 21:48  信号处理代码\mfiles\contwtgn.m

     文件       2688  2018-01-04 21:48  信号处理代码\mfiles\contwtgnmir.m

     文件       2242  2018-01-04 21:48  信号处理代码\mfiles\correlmx.m

     文件       5573  2018-01-04 21:55  信号处理代码\mfiles\CVS\Entries

     文件         12  2018-01-04 21:55  信号处理代码\mfiles\CVS\Repository

     文件         57  2018-01-04 21:55  信号处理代码\mfiles\CVS\Root

     文件       1478  2018-01-04 21:48  信号处理代码\mfiles\d2statio.m

     文件       1392  2018-01-04 21:48  信号处理代码\mfiles\disprog.m

     文件       1134  2018-01-04 21:48  信号处理代码\mfiles\divider.m

............此处省略122个文件信息