Question

求一段MATLAB代码的意思啊，急求，有知道的大侠给翻译下吧

%¨x(t),100Hz,1000
N=10;
Fs=100;
Ts=1;
dt=Ts/Fs;
a=randint(1,N,[1 10]);
x=zeros(1,N*Ts/dt);
for i=1:length(x)
x(i)=a(ceil(i/Ts*dt));
end
ft=1024;
xw=fft(x,ft);f=0:1/ft:1-1/ft;
%
hw1=exp(-j*f*2*pi);
yw1=hw1.*xw;
yt1=ifft(yw1,ft);
figure(1);
subplot(2,1,1);
plot(abs(hw1));title('');
subplot(2,1,2);
plot(angle(hw1));title('');
figure(2);
subplot(2,1,1);plot(x);title('');
axis([1 1200 0 20]);
subplot(2,1,2);plot(real(yt1));title('');
axis([1 1200 0 20]);
%
hw2=zeros(1,ft);
hw2=sin(f*pi)./(f*pi);
hw2=hw2+sin((1-f)*2*pi)./((1-f)*2*pi);
yw2=hw2.*hw1.*xw;
yw2(1)=0;
figure(3);
subplot(2,1,1);
plot(abs(hw2));title('');
subplot(2,1,2);
plot(angle(hw2));title('');
yt2=ifft(yw2,ft);
figure(4);
subplot(2,1,1);plot(x);title('');
axis([1 1200 0 20]);
subplot(2,1,2);plot(real(yt2));title('');
%
hw3=zeros(1,ft);
hw3(1:ft/2)=exp(-j*(2*pi*f(1:ft/2)-pi));
hw3(ft/2+1:ft)=exp(-j*(2*pi*f(ft/2+1:ft)+pi));
yw3=hw3.*xw;
yt3=ifft(yw3,ft);
figure(5);
subplot(2,1,1);
plot(abs(hw3));title('');
subplot(2,1,2);
plot(angle(hw3));title('');
figure(6);
subplot(2,1,1);plot(x);title('');
axis([1 1200 -20 20]);
subplot(2,1,2);plot(real(yt3));title('');
axis([1 1200 -20 20]);

Answer 1

%¨x(t),100Hz,1000
N=10;
Fs=100;
Ts=1;
dt=Ts/Fs;
%% 随机产生N个1-10的整数
a=randint(1,N,[1 10]);
%% 用a初始化N*Ts/dt个数据 即1000个
x=zeros(1,N*Ts/dt);
for i=1:length(x)
x(i)=a(ceil(i/Ts*dt));
end
%% 对X进行离散傅立叶变换
ft=1024;
xw=fft(x,ft);
%% xw点乘exp(-j*f*2*pi)后进行反离散傅立叶变换
f=0:1/ft:1-1/ft;
hw1=exp(-j*f*2*pi);
yw1=hw1.*xw;
yt1=ifft(yw1,ft);
%% 显示部分
figure(1);
% 第一幅图
subplot(2,1,1);
plot(abs(hw1));title('hw1的绝对值 ');
subplot(2,1,2);
plot(angle(hw1));title('hw1的相位角');
figure(2);
% 第二幅图
subplot(2,1,1);plot(x);title('x');
axis([1 1200 0 20]);
subplot(2,1,2);plot(real(yt1));title('yt1的实部');
axis([1 1200 0 20]);
%
hw2=zeros(1,ft);% 初始化hw2为0
hw2=sin(f*pi)./(f*pi);%hw2=sin(x)/x
hw2=hw2+sin((1-f)*2*pi)./((1-f)*2*pi);%hw2=sin(x)/x+sin((1-f)*2*pi)./((1-f)*2*pi)
yw2=hw2.*hw1.*xw;%yw2=hw2点乘hw1点乘xw
yw2(1)=0;%把第一个数据改为0
% 第三幅图
figure(3);
subplot(2,1,1);
plot(abs(hw2));title('hw2的绝对值');
subplot(2,1,2);
plot(angle(hw2));title('hw2的相位角');
yt2=ifft(yw2,ft);
figure(4);
subplot(2,1,1);plot(x);title('x');
axis([1 1200 0 20]);
subplot(2,1,2);plot(real(yt2));title('yt2的实部');
%
hw3=zeros(1,ft);
hw3(1:ft/2)=exp(-j*(2*pi*f(1:ft/2)-pi));
hw3(ft/2+1:ft)=exp(-j*(2*pi*f(ft/2+1:ft)+pi));
yw3=hw3.*xw;
yt3=ifft(yw3,ft);
figure(5);
subplot(2,1,1);
plot(abs(hw3));title('hw3的绝对值');
subplot(2,1,2);
plot(angle(hw3));title('hw3的相位角');
figure(6);
subplot(2,1,1);plot(x);title('x');
axis([1 1200 -20 20]);
subplot(2,1,2);plot(real(yt3));title('yt3的实部');
axis([1 1200 -20 20]);

Answer 2

知道是解决什么问题的代码吗？看起来有点像采样，频谱分析之类的