实验四平稳时间序列模型预测.docxVIP

随机信号分析实验报告第 随机信号分析实验报告 第 PAGE 2 页 共 2 页 实验四 平稳时间序列模型预测 一、实验目的 1、掌握平稳时间序列分析模型的分析方法和步骤 2、会求平稳时间序列的自相关函数和偏相关函数 3、掌握模型类别和阶数的确定 二、实验设备 计算机、Matlab 软件 三、实验内容与步骤 ?? numberZi已知平稳时间序列{ k }一个长为 50 number Zi 1-10 289 285 289 286 288 287 288 292 291 291 11-20 292 296 297 301 304 304 303 307 299 296 21-30 293 301 293 301 295 284 286 286 287 284 31-40 282 278 281 278 277 279 278 270 268 272 41-50 273 279 279 280 275 271 277 278 279 285 51-60 301 295 281 278 278 270 286 288 279 279 每个同学以自己的学号为起点，循环计数 50 重新排序，如：学号为 3 的学生样本数据为：Z3,Z4……Z50,Z1,Z2,编程计算，并打印下列： t1、W t  ? Z ? Z t 2、W 的样本自相关函数?? 2、 t k 3、利用递推公式计算样本的偏相关系数 4、 ?? 和?? 曲线 k kk 5、确定模型的类别和阶数 四、实验原理 平稳时间序列的模型估计与预测原理 ? W W ? W W ??? W W 1 ?n?k ? ? 1 1?k 2 2?k n?k n ? W W 样本自协方差函数： k ? ?? n n j ?1 j j ?k ? 样本自相关函数： k ? k ?? 0 样本偏相关函数 ??? ? 11 ? ?? 1 ?? ?? ?k ? ? ?? ?k ? ? ??1 ?? ? k ?1 k ?1 ? ?? ? ? k ?1 ? ? k ?1 ??1 ? k j ?? ? ? ? k j ? j ?1 j ?1 ???? ?  k ?1 j ? ?? j ? ?? ?? k ?1 k ?1  k k ?( j ?1) , j ? 1,2,?, k k kk3、利用 ?? 与?? k kk 五、实验报告要求 1、写出详细的计算步骤及设计原理； 2、按实验内容的要求打印图形； 3、附上程序和必要的注解。 六．实验过程 function y = experiment4 close all;clc; % r = [];p1 = [];p = []; % Fai = [];FAI = []; %学号 21 z1 = [293 301 293 301 295 284 286 286 287 284]; z2 = [282 278 281 278 277 279 278 270 268 272]; z3 = [273 279 279 280 275 271 277 278 279 285]; z4 = [301 295 281 278 278 270 286 288 279 279]; z5 = [289 285 289 286 288 287 288 292 291 291]; z6 = [292 296 297 301 304 304 303 307 299 296]; Z = [z1 z2 z3 z4 z5 z6]; W = Z - mean(Z); figure(1), subplot(211),plot(Z);grid on; subplot(212),plot(W);grid on; N = length(W); %利用公式来求样本的自协方差函数，取K60/4 K = 15; for k = 1:K sum = 0; for i = 1:(N-k) sum = sum + W(i)*W(i+k); end %55 end r(k) = sum/N; sum = 0; for i = 1:N sum = sum + W(i)*W(i); end r0 = sum/N;% 样本方差 p1 = r/r0; p = [1 p1]; %样本相关系数 %利用递推法求偏相关函数 Fai(1,1) = p1(1); %利用公式 1 for k = 1:K - 1 sum1 = 0; sum2 = 0; for j = 1:k sum1 = sum1 + p1(k + 1)*Fai(k,j); sum2 = sum2 + p1(j)*Fai(k,j); end Fai(k + 1,k + 1) = (p1(k + 1) - sum1)/(1