1.m

通原实验课上的作业

close all	%关闭所有的窗口
clear all	%清除所有的变量
clc   %清空命令窗口
k=input('取样点数=2^k,k=[11]');  %默认取样点数=2^11
if isempty(k),k=11;
end;
N=2^k;
fs=14;   %单位为（KHz）
dt=1/fs;
T=N*1/fs;
Bs=fs/2;
t=linspace(0,T,N);
f=linspace(-Bs,Bs,N);
fm=1;
m=sin(2*pi*fm*t);
c=sin(2*pi*6*fm*t);
s=m.*c;
S=t2f(s,dt);   %傅立叶变换
figure(1)	%开窗口1
%设置窗口位置和大小
set(1,'position',[10,50,300,200])
plot(t,s,'r');
grid on		%显示网格
%设置坐标范围
axis([0,50*T/length(t),min(s),max(s)])
xlabel('t(ms)')	%给X轴加标签
ylabel('s(t)')	%给Y轴加标签
%画频谱曲线
figure(2)
set(2,'position',[350,50,300,200])
plot(f,abs(S),'g')
grid on
axis([-Bs-3,Bs+3,min(abs(S)),max(abs(S))])
xlabel('f(kHz)')
ylabel('S(f)')