t=[-20:0.1:20];%syms?t;
f=10*exp（i*t.^2）;
Faf1=frft（f1.58）;
subplot（321）;
plot（tabs（f））;
title（‘輸入函數(shù)‘）;
xlabel（‘t‘）;
ylabel（‘f‘）
subplot（323）;
plot（tabs（Faf1））;
title（‘分?jǐn)?shù)階傅立葉變換輸出結(jié)果‘）;
xlabel（‘t‘）;
ylabel（‘a(chǎn)bs（Faf1）‘）
Faf2=frft（Faf1-1.58）;%分?jǐn)?shù)階傅立葉反變換
subplot（325）;
plot（tabs（Faf2））;
title（‘再進(jìn)行分?jǐn)?shù)階傅立葉反變換后輸出結(jié)果‘）;
xlabel（‘t‘）;
ylabel（‘a(chǎn)bs（Faf2）‘）
f=awgn（f2）;%加入高斯白噪聲
Faf1=frft（f1.58）;
subplot（322）;
plot（tabs（f））;
xlabel（‘t‘）;
ylabel（‘f‘）
title（‘加入高白的輸入函數(shù)‘）;
subplot（324）;
plot（tabs（Faf1））;
xlabel（‘t‘）;
ylabel（‘a(chǎn)bs（Faf1）‘）
title（‘加入高白的分?jǐn)?shù)階傅立葉變換輸出結(jié)果‘）;
Faf2=frft（Faf11.58）;%分?jǐn)?shù)階傅立葉反變換
subplot（326）;
plot（tabs（Faf2））;
title（‘再進(jìn)行分?jǐn)?shù)階傅立葉反變換后輸出結(jié)果（加入高白）‘）;
xlabel（‘t‘）;
ylabel（‘a(chǎn)bs（Faf2）‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-04-10?21:02??frft\
?????文件?????????841??2013-04-10?20:13??frft\chirp.asv
?????文件?????????840??2013-04-10?20:13??frft\chirp.m
?????文件?????????160??2013-03-31?20:41??frft\danpin.asv
?????文件?????????245??2013-04-10?21:02??frft\danpin.m
?????文件?????????254??2013-04-01?17:42??frft\duopin.m
?????文件????????1706??2013-03-29?16:40??frft\frft.asv
?????文件????????1706??2013-03-29?16:43??frft\frft.m
?????文件?????????392??2013-04-01?15:35??frft\gaobai.asv
?????文件?????????392??2013-04-01?15:35??frft\gaobai.m