代码搜索结果

% script that generates the figure 3.4 of Carrier modulation chapter snrindB=0:0.25:18; for i=1:length(snrindB), snr=10^(0.1*snrindB(i)); P2(i)=Qfunct(sqrt(2*snr)); P4(i)=2*Qfunct(sqrt(4*sn

%Samp3_14 clf;fs=100; %采样频率 Ndata=32; %数据长度 N=32; %FFT的数据长度 n=0:Ndata-1;t=n/fs; %数据对应的时间序列 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %时间域信号 y=fft(x,N); %信号的Fourier变换 mag=abs(y); %求取振幅

%Samp6_12 Order=10; %滤波器的阶数 f=0:0.1:1; %归一化频率点 m=[0 0 1 1 0 0 1 1 1 0 0]; %幅度点 [b,a]=yulewalk(Order,f,m); %设计滤波器 [h,w]=freqz(b,a,128); %计算128个点的频率特性 figure(1) plot(f,m,'b-',w/pi,abs(h)

%Samp5_13 N=5;Rp=3;f1=100;f2=500; %滤波器阶数、边界频率（Hz） w1=2*pi*f1;w2=2*pi*f2; %边界频率（rad/s） [z,p,k]=cheb1ap(N,Rp); %设计Chebyshev I型原型低通滤波器 [b,a]=zp2tf(z,p,k); %转换为传递函数形式 Wo=sqrt(w1*w2); %中心频率

%Samp2_4 n=[0:100]; %给出序号序列 x=2*sin(0.02*pi*n+pi/4); %给出值序列 stem(n,x) %绘值离散图 xlabel('n');ylabel('x(n)');title('正弦序列') %必要标记 grid on; %添加网

%Samp5_18 wp=1500*2*pi;ws=1000*2*pi;Rp=1;Rs=30; %滤波器设计参数 [N,Wn]=ellipord(wp,ws,Rp,Rs,'s'); %求滤波器的最小阶数和截止频率 [b,a]=ellip(N,Rp,Rs,Wn,'high','s'); %设计高通椭圆滤波器 w=linspace(1,3000,1000)*2*pi; %给出计算

%Samp4_1 clf; b=[1 2 1]; a=[1 0.5 0.25]; %传递函数分子和分母多项式系数 m=0:length(b)-1; n=0:length(a)-1; %序列序号 K=512; %采用512个点绘制幅频特性 k=1:K;w=pi*k/K; %转换为归一化频率，Nyquist频率对应于pi。 num=b*exp(-j*m'*w);

%Samp7_6 wp=[0.35 0.65];N=48; %通带边界频率（归一化频率）和滤波器阶数 Fs=50; b=fir1(N,wp); %设计FIR带通滤波器 figure(1) [H,f]=freqz(b,1,512,Fs); %以50Hz为采样频率求出滤波器频率响应 subplot(2,1,1),plot(f,20*log10(abs

% Samp3_17 N1=10;dt=1; %采样点数和采样间隔 n1=0:N1-1;t1=n1*dt; %时间序列 x=cos(2*pi*0.24*t1)+cos(2*pi*.26*t1); %原始信号 subplot(2,3,1);stem(t1,x); hold on; plot([0 10],[0 0]); %绘制横轴 title('信号x(n), n=0

close all;figure(1) ws=[1000 2000]*2*pi;wp=[500 2500]*2*pi;Rp=1;Rs=50; %滤波器设计参数，乘以2 转换为弧度/s [N,Wn]=cheb1ord(wp,ws,Rp,Rs,'s'); %求模拟滤波器的最小阶数和截止频率 w=linspace(1,3000,1000)*2*pi; %设置计算复数频率响应的频率点 [