exe0115.m

来自「Matlab程序为清华大学胡广书版信号处理的上机习题答案」· M 代码 · 共 134 行

clear;


n=[1770:1869];
x=[
    101
    82
    66
    35
    31
    7
    20
    92
    154
    125
    85
    68
    38
    23
    10
    24
    83
    132
    131
    118
    90
    67
    60
    47
    41
    21
    16
    6
    4
    7
    14
    34
    45
    43
    48
    42
    28
    10
    8
    2
    0
    1
    5
    12
    14
    35
    46
    41
    30
    24
    16
    7
    4
    2
    8
    17
    36
    50
    62
    67
    71
    48
    28
    8
    13
    57
    122
    138
    103
    86
    63
    67
    24
    11
    15
    40
    62
    98
    124
    96
    66
    64
    54
    39
    21
    7
    4
    23
    55
    94
    96
    77
    59
    44
    47
    30
    16
    7
    37
    74];
x=x';
M=30;
rx=xcorr(x,M,'biased');
subplot(221);
stem(n,x);
xlabel('n');
ylabel('x(n)');%grid on;

subplot(223);
stem((-M:M),rx);%grid on;
xlabel('m');ylabel('rx(m)');

%信号除去均值
umean=sum(x)/length(x);
xu=x/umean;
rxu=xcorr(xu,M,'biased');
subplot(222);
stem(n,xu);
xlabel('n');
ylabel('xu(n)');%grid on;

subplot(224);
stem((-M:M),rxu);%grid on;
xlabel('m');ylabel('rxu(m)');