最新信號與系統(tǒng)課程設(shè)計(jì)應(yīng)用matlab實(shí)現(xiàn)連續(xù)信號的采樣與重構(gòu)仿真(編輯修改稿)
2025-07-23 13:56 本頁面
【文章內(nèi)容簡介】 頻譜，并觀察輸出圖形，并做出分析。具體程序以及步驟如下： 周期信號頻譜% 周期與頻譜的關(guān)系實(shí)現(xiàn)程序function [A_sym,B_sym]=CTFSingsym% 采用符號計(jì)算求[0,T]內(nèi)時間函數(shù)的三角級數(shù)展開系數(shù)% 函數(shù)的輸入輸出都是數(shù)值量% Nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù)% A_sym 第1元素是直流項(xiàng)，其后元素依次是1,2,3,...次斜波cos項(xiàng)展開系數(shù)% B_sym 第2,3,4,...元素依次是1,2,3,...次斜波sin項(xiàng)展開系數(shù)% T T=m*tao, 周期信號% Nf 諧波的階數(shù)% Nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù)% m(m=T/tao)周期與脈沖寬度之比，如m=4,8,16,100等% tao 脈寬:tao=T/Msyms t n y 沈 陽 大 學(xué)課程設(shè)計(jì)說明書 if nargin3。Nf=input(39。pleas Input 所需展開的最高諧波次數(shù)：Nf=39。)。endT=input(39。pleas Input 信號的周期T 39。)。if nargin5。Nn=32。endy=time_fun_s(t)。A0=2*int(y,t,0,T)/T。As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T)。Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T)。A_sym(1)=double(vpa(A0,Nn))。for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn))。B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn))。endif nargout==0 S1=fliplr(A_sym) S1(1,k+1)=A_sym(1) S2=fliplr(1/2*S1) S3=fliplr(1/2*B_sym) S3(1,k+1)=0 S4=fliplr(S3) S5=S2i*S4。 N=Nf*2*pi/T。 k2=0:2*pi/T:N。 x=time_fun_e(t) subplot 212 stem(k2,abs(S5))。 title(39。連續(xù)時間函數(shù)周期矩形脈沖的單邊幅度譜39。) axis([0,80,0,]) line([0,80],[0,0])end %function y=time_fun_s(t)% 。它是由符號變量和表達(dá)式寫成syms a a1 沈 陽 大 學(xué)課程設(shè)計(jì)說明書 T=input(39。
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1