function [ output_args ] = ins( input_args )
%UNTITLED1 Summary of this function goes here
%   Detailed explanation goes here

% 原始信號(hào)
N=128;
fs=1;
f1=0.1;
w=2*pi*f1/fs;
data=sin(w*[0:N-1]);

% 0 值插值
datainsert(4*N)=0;
k=1;
while k<=N
    datainsert(4*(k-1)+1)=data(k);
    k=k+1;
end
datafft=fftshift(abs(fft(datainsert)));
fx2=-0.5:1/(4*N):0.5-1/(4*N);
fx2=fx2.*(4*fs/1000);
subplot(5,2,3);
plot(fx2,datafft);
title('4 倍插值，采樣率變?yōu)?38k*4');
subplot(5,2,4);
plot(datainsert);
title('時(shí)域波形');

% 低通濾波
b=fir1(64,0.2); % 本來應(yīng)該是 0.25 的
datafilter=filter(b,1,datainsert);
datafft=fftshift(abs(fft(datafilter)));
subplot(5,2,5);
plot(fx2,datafft);
title('低通濾波后');
subplot(5,2,6);
plot(datafilter);
title('時(shí)域波形');

% 調(diào)制
w=2*pi*fs/(4*fs);
s=sin(w*[0:4*N-1]);
datas=datafilter.*s;
datafft=fftshift(abs(fft(datas)));
subplot(5,2,7);
plot(fx2,datafft);
title('調(diào)幅調(diào)制到 38khz');
subplot(5,2,8);
plot(datas);
title('時(shí)域波形');

% 插入導(dǎo)頻
w=2*pi*(fs/2)/(4*fs);
plevel=0.0003;
s=plevel*sin(w*[0:4*N-1]);
datasum=datafilter+datas+s;
datafft=fftshift(abs(fft(datasum)));
subplot(5,2,9);
plot(fx2,datafft);
title('插入導(dǎo)頻和左+右信號(hào)');
subplot(5,2,10);
plot(datasum);
title('時(shí)域波形');