enegy_sqrt.m

来自「本程序实现产生几种调制方式的眼图并且对其进行比较」· M 代码 · 共 30 行

[h2,den] = rcosine(1,100,'fir/sqrt',0.5,100);
t=-100:0.01:100;
figure(5);
stem(t,h2,'.');
title('平方根升余弦函数序列');
grid on;
b=sum(h2.*h2);


t=-6:0.05:6;
h1=sin(pi*t)./(pi*t);
%——由于计算机没有0/0的定义，所以要另外附加此点的定义
h1(121)=1;
figure(4);
stem(t,h1,'.');
title('sinc滤波器函数序列');
grid on;
c=sum(h1.*h1)*0.05;

t=-2000:0.01:2000;
h=sin(pi*t)./(pi*t).*cos(pi/2*t)./(1-t.^2);
%——由于计算机没有0/0的定义，所以要另外附加此点的定义
h((length(t)-1)/2+1)=1;
h(199901)=0;
h(200101)=0;
figure(1);
stem(t,h,'.');
title('r=0.5的升余弦滤波器函数序列');
grid on;
d=sum(h.*h)*0.01;