資源簡介
clear all;
close all;
fd=1e3;
fs=fd*10;
delay=5;
%分別求四種情況下的升余弦函數
num1=rcosine(fd,fs,'fir/normal',0,delay);
num2=rcosine(fd,fs,'fir/normal',0.5,delay);
num3=rcosine(fd,fs,'fir/normal',0.75,delay);
num4=rcosine(fd,fs,'fir/normal',1,delay);
%分別求四種情況下的升余弦函數的幅度響應
Hwa=abs(fft(num1,1000));
Hwb=abs(fft(num2,1000));
Hwc=abs(fft(num3,1000));
Hwd=abs(fft(num4,1000));
angH1=angle(fft(num1,1000));
angH2=angle(fft(num2,1000));
angH3=angle(fft(num3,1000));
angH4=angle(fft(num4,1000));
t=0:1/fs:1/fs*(length(num1)-1); %分別畫出四種情況下的升余弦對應的時域波形
figure(1);
代碼片段和文件信息
clear?all;
close?all;
fd=1e3;
fs=fd*10;
delay=5;
?????????????????????????????????????????????????%分別求四種情況下的升余弦函數
?num1=rcosine(fdfs‘fir/normal‘0delay);
?num2=rcosine(fdfs‘fir/normal‘0.5delay);
?num3=rcosine(fdfs‘fir/normal‘0.75delay);
?num4=rcosine(fdfs‘fir/normal‘1delay);
????????????????????????????????????????????????%分別求四種情況下的升余弦函數的幅度響應
????Hwa=abs(fft(num11000));
????Hwb=abs(fft(num21000));
????Hwc=abs(fft(num31000));
????Hwd=abs(fft(num41000));
????angH1=angle(fft(num11000));
????angH2=angle(fft(num21000));
????angH3=angle(fft(num31000));
????angH4=angle(fft(num41000));
????
????
????t=0:1/fs:1/fs*(length(num1)-1);??????????????%分別畫出四種情況下的升余弦對應的時域波形
????figure(1);?
????plot(tnum1‘-‘);hold?on;text(0.00363-0.21‘\rightarrowα=0‘?‘FontSize‘11)
????plot(tnum2‘:‘);text(0.0034-0.10‘\rightarrowα=0.5‘?‘FontSize‘11);hold?on;
????plot(tnum3‘-.‘);text(0.00358-0.05‘\rightarrowα=0.75‘?‘FontSize‘11);hold?on;
????plot(tnum4‘--‘);text(0.003520.006‘\rightarrowα=1‘?‘FontSize‘11);hold?on;
????ylabel(‘四種升余弦對應的時域波形‘)?;?xlabel(‘時間/s‘)?;title(‘-實線表α=0:點線表α=0.5—.點劃線表α=0.75—?—虛線表α=1‘);
????grid?on;
??
????f=(1:fs/1000:fs)-1;????????????????????????????%分別畫出四種情況下的升余弦對應的幅頻響應
????figure(2);
????plot(fHwa‘-‘);text(43610.2‘\rightarrowα=0‘?‘FontSize‘11);hold?on;
????plot(fHwb‘:‘);text(3519.0‘\rightarrowα=0.5‘?‘FontSize‘11);hold?on;
????plot(fHwc‘-.‘);text(3208.4‘\rightarrowα=0.75‘?‘FontSize‘11);hold?on;
????plot(fHwd‘--‘);text(1919.2‘\rightarrowα=1‘?‘FontSize‘11);hold?on?;
????axis([0?1500?0?12]);ylabel(‘四種升余弦對應的幅度譜‘);xlabel(‘頻率/Hz‘)?;
title(‘-實線表α=0:點線表α=0.5—.點劃線表α=0.75—?—虛線表α=1‘);grid?on;
figure(3);
subplot(221)
????plot(fangH1);axis([0?1500?0?4]);?title(‘α=0‘);ylabel(‘四種升余弦對應的幅度譜‘);xlabel(‘頻率/Hz‘)?;grid?on;
????subplot(222)
????plot(fangH2);axis([0?1500?0?4]);title(‘α=0.5‘);xlabel(‘頻率/Hz‘)?;grid?on;
????subplot(223)
????plot(fangH3);axis([0?1500?0?4]);title(‘α=0.75‘);xlabel(‘頻率/Hz‘)?;grid?on;
????subplot(224)
????plot(fangH4);axis([0?1500?0?4]);title(‘α=1‘);xlabel(‘頻率/Hz‘)?;grid?on;
????
?
N=2000;?????????????????????????????%數字序列長度
????N_sample=10;????????????????????%每碼元抽樣次
- 上一篇:多重網格法求解微分方程-matlab
- 下一篇:音樂專輯數據庫用access做的
評論
共有 條評論