 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 380 | 回復(fù): 0 | ||
hailiang金蟲(chóng) (正式寫(xiě)手)
|
[求助]
請(qǐng)教，對(duì)信號(hào)做傅里葉分析，當(dāng)原始信號(hào)不是整數(shù)或加0.5時(shí)結(jié)果不對(duì)
|
|
請(qǐng)教，對(duì)信號(hào)做傅里葉分析，當(dāng)原始信號(hào)不是整數(shù)或加0.5時(shí)結(jié)果不對(duì) 我用下面的程序?qū)τ嘞倚盘?hào)做頻譜分析，當(dāng)信號(hào)頻率是整數(shù)或者是正負(fù)0.5的時(shí)候都正確，就是求出的振幅和頻率剛好是原始信號(hào)一樣，但是如果原始信號(hào)頻率不是整數(shù)或正負(fù)0.5時(shí)就不正確，請(qǐng)問(wèn)是精度設(shè)置的不夠高，還是采樣頻率和采樣點(diǎn)數(shù)不對(duì)！ 下面的程序是Mathematica程序，Matlab程序試了也一樣 A1 = 3; (*% 頻率F1信號(hào)的幅度*) (*F1=1/(2 Pi); (*% 信號(hào)1頻率 (Hz)*)*) F1 = 0.11; Fs = 256; (*% 采樣頻率 (Hz)，要大于信號(hào)頻率的兩倍*) P1 = 100; (*% 信號(hào)1相位 (度)*) NN = 2048; (*% 采樣點(diǎn)數(shù)，通常取2的整數(shù)次方*) S = A1*Cos[2*Pi*F1*t + Pi*P1/180]; (*輸入原始信號(hào)*) Plot[S, {t, 0, 10}, PlotPoints -> 200, PlotLabel -> "原始信號(hào)"] Plot[S, {t, 0, NN/Fs}, PlotPoints -> 200, PlotLabel -> "原始信號(hào)"] b = Range[0, NN/Fs, 1/Fs]; g = S /. t -> b; S1 = Take[g, Length[g] - 1]; Y = Fourier[S1, FourierParameters -> {1, -1}]; Ayy = Abs[Y]; ListLinePlot[Ayy, PlotRange -> All] (*以下顯示換算后的FFT模值結(jié)果*) Ayy = Ayy/(NN/2); Ayy[[1]] = Ayy[[1]]/2; (*Ayy[[1]]=0;*) F = Range[1, NN]; F = (F - 1)*Fs/NN; FF = Array[h3, {NN/2, 2}]; Do[h3[i, 1] = F[]; h3[i, 2] = Ayy[], {i, 1, NN/2}]; ListLinePlot[FF, PlotRange -> All, Frame -> True, FrameLabel -> {"頻率/Hz", "振幅"}] 下面是matlab程序 close all; %先關(guān)閉所有圖片 Adc=2; %直流分量幅度 A1=3; %頻率F1信號(hào)的幅度 F1=0.1; %信號(hào)1頻率(Hz) Fs=256; %采樣頻率(Hz) P1=-30; %信號(hào)1相位(度) N=1024; %采樣點(diǎn)數(shù) t=[0:1/Fs:N/Fs]; %采樣時(shí)刻 %信號(hào) S=A1*cos(2*pi*F1*t+pi*P1/180); %顯示原始信號(hào) plot(S); title('原始信號(hào)'); figure; Y = fft(S,N); %做FFT變換 Ayy = (abs(Y)); %取模 plot(Ayy(1:N)); %顯示原始的FFT模值結(jié)果 title('FFT 模值'); figure; Ayy=Ayy/(N/2); %換算成實(shí)際的幅度 Ayy(1)=Ayy(1)/2; F=([1:N]-1)*Fs/N; %換算成實(shí)際的頻率值 plot(F(1:N/2),Ayy(1:N/2)); %顯示換算后的FFT模值結(jié)果 title('幅度-頻率曲線圖'); Y1=Ayy(1:N/2); [da1,dda1]=max(Y1); Y1(dda1)=0; [da2,dda2]=max(Y1); Y1(dda2)=0; [da3,dda3]=max(Y1); da1 da2 da3 figure; Pyy=[1:N/2]; for i=1:N/2 Pyy(i)=phase(Y(i)); %計(jì)算相位 Pyy(i)=Pyy(i)*180/pi; %換算為角度 end; plot(F(1:N/2),Pyy(1:N/2)); %顯示相位圖 title('相位-頻率曲線圖'); |
回復(fù)此樓» 猜你喜歡
328求調(diào)劑
已經(jīng)有7人回復(fù)
-
340求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有15人回復(fù)
-
291求調(diào)劑
已經(jīng)有19人回復(fù)
-
331環(huán)境科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有8人回復(fù)
-
315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有6人回復(fù)
-
291求調(diào)劑
已經(jīng)有7人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
請(qǐng)教一下JPSP的申請(qǐng)
已經(jīng)有4人回復(fù)
- 新手提RNA，請(qǐng)教幾個(gè)問(wèn)題
已經(jīng)有8人回復(fù)
- 陶瓷信回復(fù)請(qǐng)教
已經(jīng)有10人回復(fù)
- 有在澳洲悉尼的嗎，請(qǐng)教個(gè)問(wèn)題哦
已經(jīng)有3人回復(fù)
- 請(qǐng)教一道物化題
已經(jīng)有7人回復(fù)
- 請(qǐng)教 ĸ ŋ配位模式的區(qū)別
已經(jīng)有3人回復(fù)
- 求南京理工2010和2011年《信號(hào)系統(tǒng)與數(shù)字電路》考研真題及答案
已經(jīng)有4人回復(fù)
- 沒(méi)有信號(hào)肽會(huì)影響蛋白的表達(dá)過(guò)程
已經(jīng)有9人回復(fù)
- 信號(hào)與線性系統(tǒng)分析課件 第四版
已經(jīng)有44人回復(fù)
- 光學(xué)信號(hào)cos(sin(2pi*f*t))的傅里葉展開(kāi)？
已經(jīng)有6人回復(fù)
- 信號(hào)處理方面世界上比較權(quán)威的雜志有哪些，麻煩大家略微介紹一下這些雜志
已經(jīng)有5人回復(fù)
- 材料拉曼信號(hào)的數(shù)據(jù)庫(kù)
已經(jīng)有3人回復(fù)
- 粉末xrd數(shù)據(jù)，jade分析,txt比raw的信號(hào)強(qiáng)是為什么
已經(jīng)有11人回復(fù)
- 【求助】自制的Faraday cage(法拉第籠)測(cè)試時(shí), 屏蔽不了手機(jī)信號(hào).求好心人幫忙解惑...
已經(jīng)有6人回復(fù)
- 【求助/交流】測(cè)序信號(hào)弱
已經(jīng)有6人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hailiang 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +4 | hanamiko 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 289求調(diào)劑 +8 | 新時(shí)代材料 2026-03-27 | 8/400 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +4 | 麻辣魷魚(yú) 2026-03-27 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | 白QF 2026-03-21 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +3 | 26要上岸 2026-03-27 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 311求調(diào)劑 +7 | lin0039 2026-03-26 | 7/350 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛(ài)學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|