trianglewave.m

来自「Galileo BOC 码产生捕获相关仿真结果」· M 代码 · 共 21 行

clear;

N=1024;
n=1:1:N;
y=zeros(1,N);

for i=1:N,
    y(i)=7*((-1)^fix(n(i)/14))+((-1)^(fix(n(i)/14)+1))*(n(i)-14*fix(n(i)/14));
end

yy=zeros(1,32);
for i=1:32,
    yy(i)=y(i);
end

freq_y=fft(y);
mag_y=abs(freq_y);
subplot(2,1,1);
plot(n,y)
subplot(2,1,2);
plot(n,mag_y)