版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

yjz001

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 397.8
散金: 30
帖子: 29
在線: 9.1小時
蟲號: 2857323
注冊: 2013-12-08
專業(yè): 通信理論與系統(tǒng)

[求助] 對不同取值范圍的信號序列做分?jǐn)?shù)階傅立葉變換（frft）是否要相應(yīng)的改變frft程序？？？已有1人參與

以下是frft的matlab程序代碼，當(dāng)線性調(diào)頻（LFM）信號表達(dá)式是for n=-255:255； x(n+256)=exp(j*fi+j*2*pi*f0*n/fs+j*pi*k*((n/fs)*(n/fs))); end 時信號序列的取值范圍是-255到255，可用下面的程序進(jìn)行frft估計調(diào)頻率和頻率，而若將n的取值范圍改成0到510時，需不需要改動下面的程序？若需要改動，該怎么修改？（我之前覺得只需要改shft的表達(dá)式，即改成shft=（0：N）+1，但結(jié)果貌似不正確），請懂的朋友幫忙解惑一下~~

CODE:

function Faf = frft(f, a)

% The fast Fractional Fourier Transform快速分?jǐn)?shù)傅里葉變換

% input: f = samples of the signal

%        a = fractional power分?jǐn)?shù)冪

% output: Faf = fast Fractional Fourier transform

error(nargchk(2, 2, nargin));

f = f(:);

N = length(f);

shft =  rem((0:N-1)+fix(N/2),N)+1;

sN = sqrt(N);

a = mod(a,4);

% do special cases

if (a==0), Faf = f; return; end;

if (a==2), Faf = flipud(f); return; end;

if (a==1), Faf(shft,1) = fft(f(shft))/sN; return; end 

if (a==3), Faf(shft,1) = ifft(f(shft))*sN; return; end

% reduce to interval 0.5 < a < 1.5

if (a>2.0), a = a-2; f = flipud(f); end

if (a>1.5), a = a-1; f(shft,1) = fft(f(shft))/sN; end

if (a<0.5), a = a+1; f(shft,1) = ifft(f(shft))*sN; end

% the general case for 0.5 < a < 1.5

alpha = a*pi/2;

tana2 = tan(alpha/2);

sina = sin(alpha);

f = [zeros(N-1,1) ; interp(f) ; zeros(N-1,1)];

% chirp premultiplication

chrp = exp(-i*pi/N*tana2/4*(-2*N+2:2*N-2)'.^2);

f = chrp.*f;

% chirp convolution

c = pi/N/sina/4;

Faf = fconv(exp(i*c*(-(4*N-4):4*N-4)'.^2),f);

Faf = Faf(4*N-3:8*N-7)*sqrt(c/pi);

% chirp post multiplication

Faf = chrp.*Faf;

% normalizing constant

Faf = exp(-i*(1-a)*pi/4)*Faf(N:2:end-N+1);

function xint=interp(x)

% sinc interpolation

N = length(x);

y = zeros(2*N-1,1);

y(1:2:2*N-1) = x;

xint = fconv(y(1:2*N-1), sinc([-(2*N-3):(2*N-3)]'/2));

xint = xint(2*N-2:end-2*N+3);

function z = fconv(x,y)

% convolution by fft

N = length([x(:);y(:)])-1;

P = 2^nextpow2(N);

z = ifft( fft(x,P) .* fft(y,P));

z = z(1:N);

回復(fù)此樓

» 猜你喜歡

298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有5人回復(fù)
297求調(diào)劑已經(jīng)有5人回復(fù)
302求調(diào)劑已經(jīng)有10人回復(fù)
266求調(diào)劑已經(jīng)有8人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有17人回復(fù)
344求調(diào)劑已經(jīng)有5人回復(fù)
312求調(diào)劑已經(jīng)有9人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有4人回復(fù)
304求調(diào)劑已經(jīng)有13人回復(fù)

1樓 2016-01-09 22:44:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yjz001

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 397.8
散金: 30
帖子: 29
在線: 9.1小時
蟲號: 2857323
注冊: 2013-12-08
專業(yè): 通信理論與系統(tǒng)

沒有人懂分?jǐn)?shù)階傅立葉變化的嗎？？有的話也好，可以討論討論也好啊

贊一下

回復(fù)此樓

2樓2016-01-10 22:30:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

木下雨

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4
帖子: 3
在線: 47分鐘
蟲號: 2888860
注冊: 2013-12-23

【答案】應(yīng)助回帖

求討論，這個算法在最佳聚焦峰值上的值不等于公式里算出來的值，有沒有誰知道為什么？

贊一下

回復(fù)此樓

3樓2016-03-11 19:35:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

忠乎

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 33.5
散金: 5
帖子: 18
在線: 3小時
蟲號: 5845665
注冊: 2017-03-05
專業(yè): 集成電路制造與封裝

樓主有二維的嘛，我想用來處理圖像啊，求助

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

4樓2017-04-12 14:28:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

半夏星辰

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 159.2
帖子: 41
在線: 5.6小時
蟲號: 6493526
注冊: 2017-05-10
專業(yè): 信號理論與信號處理

最后你懂了嗎

發(fā)自小木蟲Android客戶端

回復(fù)此樓

5樓2017-11-21 00:27:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 yjz001 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 311求調(diào)劑 +10	冬十三 2026-03-15	11/550	2026-03-18 09:04 by zhukairuo
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5	瑪卡巴卡啊哈 2026-03-11	5/250	2026-03-13 21:40 by JourneyLucky
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

yjz001

[求助] 對不同取值范圍的信號序列做分?jǐn)?shù)階傅立葉變換（frft）是否要相應(yīng)的改變frft程序？？？ 已有1人參與

» 猜你喜歡

yjz001

木下雨

【答案】應(yīng)助回帖

忠乎

半夏星辰

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 對不同取值范圍的信號序列做分?jǐn)?shù)階傅立葉變換（frft）是否要相應(yīng)的改變frft程序？？？已有1人參與