例如 第三版数字信号处理P51 -1.14习题时域离散信号的相关性研究x(n)=Asin(ωn)+u(n),其中ω=π/16,u(n)是白噪声,现要求
⑴、产生均值为0,功率P=0.1的均匀分布白噪声u(n),求u(n)自相关函数ru(m)
⑵、使x(n)的信噪比10dB决定A的数值并画出x(n)的图形及其自相关函数的图形
(1)
1 新建一个matlab脚本---edit 文件名.m
2 然后再.m后缀的文件中编写代码
3 在脚本中调试:直接打文件名+回车键就可以调试了,
4 过一会会出现如下最右边的图形
edit jaoany.m
p=0.1;N=50000;
u0=rand(1,N);u0=u0-mean(u0);
a=sqrt(p/0.083);u=u0*a;
power_u=dot(u,u)/N
plot(u(1:100));grid on;
求自相关函数
>> ru=xcorr(u,100,\'unbiased\');plot(ru) %求自相关函数
⑵、使x(n)的信噪比10dB决定A的数值并画出x(n)的图形及其自相关函数的图形
运行可以直接在>>l后面敲下代码,不过错了不好编辑,建议还是建立一个.m的文档来编写
n=0:49999;
y=0.5*pi*sin(0.0625*pi*n);
x=y+u;
plot(x(1:100));grid on;
求自相关函数>> rx=xcorr(x,100,\'unbiased\');plot(rx)%求自相关函数