无理频谱笔记

代办

添加链接描述
添加链接描述
添加链接描述

1）请导出该周期信号的无理频谱的模相式。
2）若A为10v，T为0.05s，Tp为0.02s，T为0.02s，请画出该周期信号及其无理频谱。
3）请分析该周期信号及其无理频谱的奇偶性。

%原信号
Xt = [];
i = 1;
k = [];
 for t = -0.1:0.0001:0.1
         k(i)=t;
      for m = -40:1:40
           if (t<=0.05*m-0.01)&(t>=0.05*m-0.03)
             Xt(i) = 10;
             break
           else 
             Xt(i) = 0;
           end
      end
      i = i+1;
 end
plot(k,Xt)

n = -60:1:60;

A = abs(4*sinc(0.4*n));
stem(n,A)        %模

Re = 4*sinc(0.8*n);
stem(n,Re);       %实部

Im = 4*(1-cos(pi*0.8*n) )  ./   (0.8*n*pi);
stem(n,Im);       %虚部

O = pi.*(fix(0.4*n)+0.8*pi*n);
stem(n,O)       %相

用 matlab将一组离散点周期化构造

网上找到的一篇是关于正半轴的周期化，但是偶函数对于负半轴用mod 函数取余就会出错。
目前还没找到简单的方法… 自己写的太麻烦了

a = 0.02;
m = [];
i = 1;
T = 0.4;

for t = -1.6:0.01:1.6

     if t<0 
      m(i) = mod(-t,T);
     else 
       m(i) = mod(t,T);
     end
     
     if  ( (t<T)&(t>T/2) ) |  ( (t>T)&(mod(t,T)>T/2)  ) 
      m(i) = mod(t,T)-T;
     end
     
     if    (t<-T) &  (mod(-t,T)>T/2)
          m(i) = mod(-t,T)-T;
     end
     
     if (t>-T)&(t<-T/2)  
       m(i) = t+T;
     end
    
    if (t<T/2)&(t>-T/2)
      m(i) =t;
    end
     
  
  x(i)=exp(-(m(i)/a).^2);
  i = i+1;
end

t = -1.6:0.01:1.6;
stem(t,x)

x为高斯函数

建立6阶巴特沃斯抗混滤波器的传递函数模拟测试系统，并以类脉冲激励法模拟测量该抗混滤波器的传递函数
巴特沃斯介绍
matlab 函数使用
Python 函数使用