信号加减

zhangchengpeng

function [f,n]=sigadd(f1,n1,f2,n2)
n=min(min(n1),min(n2)):max(max(n1),max(n2));x1=zeros(1,length(n));x2=x1;
x1(find((n>=min(n1))&(n<=max(n1))==1))=f1;
x2(find((n>=min(n2))&(n<=max(n2))==1))=f2;
f=x1+x2;

谁能解释一下这个程序  ，详细点  主要是第三行和第四行

先thanks！

tenglang

x1(find((n>=min(n1))&(n<=max(n1))==1))=f1;
把向量f 1 赋给向量x1 中下标 从n1(min) 开始到 n1(max)对应数组值,

先help find
find 返回一个(时间)向量(tnmin: tnmax);
x1(tnmin: tnmax)=f1;

1. 
   
2. % t1, t2 两信号的测量时间段， 一个从2秒开始，一个从3秒开始
   
3. % 在时间空间上叠加两个信号
   
4. %信号的起始时间不一样
   
5. function r=Signal()
   
6. t1=2:1:10 ;
   
7. t2=3:1:20 ;
   
8. y1=sin(t1);
   
9. y2=cos(t2);
   
10. %采样得到的f1,f2
    
11. signal1=struct('x',t1,'y',y1);
    
12. signal2=struct('x',t2,'y',y2);
    
13. plot(t1,y1,'bx:');
    
14. hold on;
    
15. plot(t2,y2,'ro--');
    
16. [newsignal,newsignaltime]= sigadd(signal1.y,signal1.x,signal2.y,signal2.x);
    
17. plot(newsignaltime, newsignal,'g*-')
    
18. r=0;
    
19. end

复制代码

zhangchengpeng

回复 2 # tenglang 的帖子

谢谢了   看懂了  但还理解的不是很透  
f2的值是不是少了一个

tenglang

上面的运算结果

zhouyang664

function [y,n]=sig_add(x1,n1,x2,n2)
% Implements y(n)=x1(n)+x2(n);

%n1,n2是两个元素的向量，分别代表序列的起始/终止下标

%即：n1 = [ a , b ] ; a代表x1的起始下标，b代表x2的终止下标！
n=min(min(n1),min(n2)):max(max(n1),max(n2));%找出下标的最大范围
y1=zeros(1,length(n));%预定义，赋初值
y2=y1;
y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;%将序列x1扩展为长度为n的序列y1，下标在n1之间的赋值x1
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;%将序列x2扩展为长度为n的序列y2，下标在n2之间的赋值x2
y=y1+y2;

ChaChing

x1(find((n>=min(n1))&(n<=max(n1))==1))=f1;
好像直接使用x1( find(n>=min(n1) & n<=max(n1)) )=f1;即可

还有程序是看了懂, 只是有点好奇其应用性, 难道仅是为了省储存空间而已吗!?

zhangchengpeng

回复 4 # tenglang 的帖子

明白了  谢谢   我还想问一下 ：为什么在 r2009b版本的MATLAB里，看不了demo的动画，我也装了FLASH插件
显示是
This video, which will play in your default web browser, requires Macromedia Flash Player (version 5 or later) and an Internet connection.

tenglang

回复 6 # ChaChing 的帖子

是多了一步逻辑判断，应该是程序作者受其他编程语言的影响，程序代码编写的风格问题。

ChaChing

在找序列的起始/终止下标, 使用min(n1), max(n1)...等, 好像根本多餘! min(n1)不就是n1(1), max(n1)不就是n1(end)! 难道时间序列可以不是单调递增吗?
所以练习下, 依个人习惯并无特别优化

1. function [f,n]=sigadd(f1,n1,f2,n2)
   
2. n=min(n1(1),n2(1)):max(n1(end),n2(end)); f=zeros(1,length(n));
   
3. ppp=n1(1)<=n & n<=n1(end); f(ppp)=f1;
   
4. ppp=n2(1)<=n & n<=n2(end); f(ppp)=f(ppp)+f2;

复制代码