【 MATLAB 】数字信号处理中的几个常用序列产生的MATLAB代码

举报
李锐博恩 发表于 2021/07/15 05:43:32 2021/07/15
【摘要】 【 MATLAB 】常用的离散时间序列的 Matlab 产生 上篇这篇博文介绍了例如单位样值序列、单位阶跃序列产生的脚本以及函数,这篇博文利用到的函数,将直接贴过来: 单位样值序列: function [x,n]=delta(n0,n1,n2);% generate x(n) = delta(n - n0); n1 <= n <= n2%________...

【 MATLAB 】常用的离散时间序列的 Matlab 产生

上篇这篇博文介绍了例如单位样值序列、单位阶跃序列产生的脚本以及函数,这篇博文利用到的函数,将直接贴过来:

单位样值序列:

\delta (n-n_0)=\left\{\begin{matrix} 1 & n=n_0\\0 & n\neq n_0 \end{matrix}\right.


  
  1. function [x,n]=delta(n0,n1,n2);
  2. % generate x(n) = delta(n - n0); n1 <= n <= n2
  3. %_____________________________________________
  4. %[x,n] = delta(n0, n1, n2);
  5. %
  6. n = [n1:n2];
  7. x = [(n-n0) == 0];

单位阶跃序列:

u(n-n_0)=\left\{\begin{matrix} 1, &n\geq n_0 \\0, & n\leq n_0 \end{matrix}\right.


  
  1. function [x,n]=stepseq(n0,n1,n2);
  2. % generate x(n) = u(n - n0); n1 <= n <= n2
  3. %_____________________________________________
  4. %[x,n] = stepseq(n0, n1, n2);
  5. %
  6. n = [n1:n2];
  7. x = [(n-n0) >= 0];


这篇博文给出几个序列,请使用脚本实现:

例1:

x(n)=2\delta(n+2)-\delta(n-4),-5\leq n\leq 5


  
  1. clc
  2. clear
  3. close all
  4. n = [-5:5];
  5. x = 2*delta(-2,-5,5)-delta(4,-5,5);
  6. stem(n,x);
  7. xlabel('n');ylabel('x(n)');
  8. ylim([-2,3]);

运行即可得到结果:


例2:

x(n)=n[u(n)-u(n-10)]+10e^{-0.3(n-10)}[u(n-10)-u(n-20)],0 \leq n \leq 20


  
  1. clc
  2. clear
  3. close all
  4. n = [0:20];
  5. x1 = n.*(stepseq(0,0,20) - stepseq(10,0,20));
  6. x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,0,20) - stepseq(20,0,20));
  7. x = x1+x2;
  8. stem(n,x);
  9. xlabel('n');ylabel('x(n)');
  10. ylim([-1,10]);


例3:

x(n)=cos(0.04\pi n)+0.2w(n),0 \leq n \leq 50

其中,w(n)是均值为0,方差为1的高斯随机序列。


  
  1. clc
  2. clear
  3. close all
  4. n = [0:50];
  5. x = cos(0.04 * pi * n) + 0.2*randn(size(n));
  6. stem(n,x);
  7. xlabel('n');ylabel('x(n)');


例4:

\tilde{x}=\left \{ ...,5,4,3,2,1,5,4,3,2,1,5,4,3,2,1,... \right \};-10 \leq n \leq 9


  
  1. clc
  2. clear
  3. close all
  4. n = [-10:9];
  5. x = [5,4,3,2,1];
  6. xtilde = x' * ones(1,4);
  7. xtilde = xtilde(:)';
  8. stem(n,xtilde);
  9. xlabel('n');ylabel('xtilde(n)');


南方姑娘

听到这首歌,我遇到的姑娘也是一个南方姑娘。


 

文章来源: reborn.blog.csdn.net,作者:李锐博恩,版权归原作者所有,如需转载,请联系作者。

原文链接:reborn.blog.csdn.net/article/details/83246159

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。