原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
1、
%---filter求卷积，B(Z)/A(Z)=H(Z),已知B(Z)和A(Z)，求y(n)=x(n)*h(n)-----
clear;
x=ones(100);
t=1:100;
b=[.001836,.007344,.011016,.007374,.001836];
a=[1,-3.0544,3.8291,-2.2925,.55075];
%
y=filter(b,a,x);
% 求所给系统的输出，本例实际上是求所给系统的阶跃响应；
plot(t,x,r.,t,y,k-);grid on;
ylabel(x(n) and y(n))
xlabel(n)
1、
%---filter求卷积，B(Z)/A(Z)=H(Z),已知B(Z)和A(Z)，求y(n)=x(n)*h(n)-----
clear;
x=ones(100);
t=1:100;
b=[.001836,.007344,.011016,.007374,.001836];
a=[1,-3.0544,3.8291,-2.2925,.55075];
%
y=filter(b,a,x);
% 求所给系统的输出，本例实际上是求所给系统的阶跃响应；
plot(t,x,r.,t,y,k-);grid on;
ylabel(x(n) and y(n))
xlabel(n)
第一章 产生信号，求卷积和自相关函数
1、
%信号产生
n=0:100;
%工频
f0=50;A=220;fs=400;
x1=A*sin(2*pi*f0*n/fs);
subplot(321);plot(n,x1);xlabel(n);ylabel(x1(n)) ;grid on;
%率减正弦
f0=2;A=2;alf=0.5;fs=16;
x2=A*exp(-alf*n/fs).*sin(2*pi*f0*n/fs);
subplot(323);plot(n,x2);xlabel(n);ylabel(x2(n)) ;grid on;
%谐波信号
f0=5;A1=1.0;A2=0.5;A3=0.2;fs=100;
x3=A1*sin(2*pi*f0*n/fs)+A2*sin(2*pi*2*f0*n/fs)+A3*sin(2*pi*3*f0*n/fs);
subplot(322);plot(n,x3);xlabel(n);ylabel(x3(n)) ;grid on;
%哈明窗
f0=10;fs=1000;
x4=0.54-0.46*cos(2*pi*f0*n/fs);
subplot(324);plot(n,x4);xlabel(n);ylabel(x4(n)) ;grid on;
%采样
n=-50:50;
f0=10;fs=400;
w=2*pi*f0*n/fs;
x5=sinc(w);
subplot(325);plot(n,x5);xlabel(n);ylabel(x5(n)) ;grid on;
2、
% 产生均匀分布的白噪信号，使均值为0，功率为p
%-----------------------------------------------------------------
clear;
p=0.01;
N=50000;
u=rand(1,N);
u=u-mean(u);
a=sqrt(12*p);
u1=u*a;
power_u1=dot(u1,u1)/N
subplot(211)
plot(u1(1:200));grid on;
ylabel(u(n))
xlabel(n)
3、
% 产生高斯分布的白噪信号，使功率为p，并观察数据分布的直方图
%-----------------------------------------------------------------
clear;
p=0.1;
N=500000;
u=randn(1,N);
a=sqrt(p);
u=u*a;
power_u=var(u);
subplot(211)
plot(u(1:200));grid on;
ylabel(u(n));
xlabel(n)
subplot(212)
hist(u,50);grid on;
ylabel(histogram of u(n));
4、
% 产生一 sinc 函数；
%-----------------------------------------------------------------
clear;
n=200;
stept=4*pi/n;
t=-2*pi:stept:2*pi;
y=sinc(t);
plot(t,y,t,zeros(siz
文档评论（0）