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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

yjz001

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 397.8
散金: 30
帖子: 29
在線: 9.1小時(shí)
蟲號: 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 時(shí)信號序列的取值范圍是-255到255，可用下面的程序進(jìn)行frft估計(jì)調(diào)頻率和頻率，而若將n的取值范圍改成0到510時(shí)，需不需要改動下面的程序？若需要改動，該怎么修改？（我之前覺得只需要改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ù)此樓

» 猜你喜歡

314求調(diào)劑已經(jīng)有8人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有5人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有6人回復(fù)
化工學(xué)碩306求調(diào)劑已經(jīng)有10人回復(fù)
286求調(diào)劑已經(jīng)有8人回復(fù)
材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦已經(jīng)有9人回復(fù)
311求調(diào)劑已經(jīng)有12人回復(fù)
271材料工程求調(diào)劑已經(jīng)有3人回復(fù)

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

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

半夏星辰

新蟲 (初入文壇)

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

最后你懂了嗎

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

回復(fù)此樓

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

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

查看全部 5 個(gè)回答

yjz001

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 397.8
散金: 30
帖子: 29
在線: 9.1小時(shí)
蟲號: 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)助回帖

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

贊一下

回復(fù)此樓

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

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

忠乎

新蟲 (初入文壇)

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

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

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

贊一下

回復(fù)此樓

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

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

查看全部 5 個(gè)回答

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

[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 材料工程327求調(diào)劑 +3	xiaohe12w 2026-03-11	3/150	2026-03-14 20:20 by ms629
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant

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

24小時(shí)熱門版塊排行榜

yjz001

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

» 猜你喜歡

半夏星辰

yjz001

木下雨

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

忠乎

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

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