功率譜的案例
Matlab實例：頻譜、功率譜和功率譜密度計算
本實例以高斯信號為例，計算幅度譜、相位譜、雙邊功率譜、雙邊功率譜密度、單邊功率譜、單邊功率譜密度。高斯信號的半波全寬FWHM=50ps，中心點位于2.5ns處。
MATLAB程序代碼：
%================
%作者：yoyoba
%Email：stuyou@126.com
%================
clc;
clear;
FWHM=50e-12;
%高斯信號FWHM寬度，為50ps
time_window=100*FWHM;
%高斯信號的采樣窗口寬度，該值決定了傅里葉變換后的頻率分辨率
Ns=2048;
%采樣點
dt=time_window/(Ns-1);
%采樣時間間隔
t=0:dt:time_window;
%采樣時間
gauss_time=exp(-0.5*(2*sqrt(2*log(2))*(t-2.5e-9)/FWHM).^2);
%高斯脈沖，中心位于2.5ns處。
plot(t*1e+9,gauss_time,'linewidth',2.5);
xlabel('Time/ns');
ylabel('Amplitude/V');
title('Gauss pulse');
%================
%以下計算雙邊譜、雙邊功率譜、雙邊功率譜密度
%================
gauss_spec=fftshift(fft(ifftshift(gauss_time)));
%傅里葉變換，并且進行fftshift移位操作。
展開 3-7 功率信號、能量信號及其功率譜與能量譜
從而信號的總功率為
若信號
的頻譜函數(shù)為
，則由帕塞瓦爾定理
可知
比較式(3-86)和式(3-87)，有
所以對于非周期功率信號可用功率密度函數(shù)
描述其功率的頻率特性，并稱之為功率譜函數(shù)。功率譜
是
的偶函數(shù)，它僅取決于頻譜函數(shù)的模量，而與相位無關，單位為瓦·秒 (
)。
對于周期信號，很容易求得其功率譜:
。
三、能量頻譜
對于能量信號，同樣可用密度的概念表示信號能量在各頻率點的分布情況，即定義單位頻率內(nèi)的信號能量為能量密度函數(shù)，記為
，從而信號
的總能量為
若信號f(t)的頻譜函數(shù)為F(jω)，由帕塞瓦爾定理，有
因此
故對于能量信號可用能量密度函數(shù)
描述其能量的頻率特性，并稱之為能量譜函數(shù)。能量譜
只與信號的
有關，而與相位無關，單位為焦/赫(
)。
由以上討論可知，信號的功率或能量既可在時域內(nèi)求得，也可以在頻域內(nèi)求得。它反映能量守恒定理在信號分析中的體現(xiàn)，也是信號的時域特性與頻域特性的一個重要關系。
例3-21 求信號
的能量。
展開 ansys workbench 隨機振動功率譜密度轉(zhuǎn)換公式
隨機振動功率譜密度轉(zhuǎn)換公式.pdf
Abaqus隨機響應分析中PSD的定義
隨機振動指未來任一給定時刻的瞬時值不能預先確定的機械振動，無法用確定的函數(shù)而須用概率統(tǒng)計方法定量描述其運動規(guī)律的振動，因此在進行隨機響應分析時隨機激勵以PSD（功率譜密度）的形式進行輸入。
1.PSD的的定義
功率譜密度譜是一種概率統(tǒng)計方法，是對隨機變量均方值的量度。一般用于隨機振動分析，連續(xù)瞬態(tài)響應只能通過概率分布函數(shù)進行描述，即出現(xiàn)某水平響應所對應的概率。功率譜密度的定義是單位頻帶內(nèi)的“功率”（均方值）。
功率譜密度是結(jié)構(gòu)在隨機動態(tài)載荷激勵下響應的統(tǒng)計結(jié)果，是一條功率譜密度值—頻率值的關系曲線，其中功率譜密度可以是位移功率譜密度、速度功率譜密度、加速度功率譜密度、力功率譜密度等形式。數(shù)學上，功率譜密度值—頻率值的關系曲線下的面積就是均方值，當均值為零時均方值等于方差，即響應標準偏差的平方值。
2. Abaqus中PSD的定義
Abauqs中通過Random response分析步進行基于模態(tài)的隨機響應分析。
?PSD的類型
Abaqus中支持以下類型的PSD施加：
① 集中載荷
② 分布載荷
③ 基礎運動（BASE MOTION）
BASE MOTION的類型分別為加速度、速度和位移。
關鍵字：*BASE MOTION
?PSD的定義
PSD為模型數(shù)據(jù)必須在分析步之前定義
① 定義隨機激勵PSD曲線的關鍵字為：
*PSD-DEFINITION, NAME=psd-name
② CAE環(huán)境定義PSD,如圖1所示。
其中：
TYPE = BASE 為基礎激勵
G為參考加速度，默認為1，用戶可以根據(jù)實際自行定義。
展開 
試驗洞悉工程 | 開孔建筑內(nèi)壓風洞試驗
圖6 孔口周邊外壓系數(shù)與內(nèi)壓系數(shù)時程曲線對比
3.1.3 內(nèi)外壓響應功率譜相關性
為了進一步分析開孔結(jié)構(gòu)內(nèi)壓動力響應的特性及其與孔口周邊外壓響應之間的關系，圖7 對比了在0° 正風向角和70° 斜風向角下，開孔尺寸為53 mm×53 mm 時模型內(nèi)外壓系數(shù)的功率譜曲線。從圖中可以看出，雖然內(nèi)外壓響應功率譜的同步性沒有內(nèi)外壓響應時程的同步性那么高，但是在0°正風向角下仍然表現(xiàn)出了一定的相關性，并且彼此之間相差不大。而在70°斜風向角下，內(nèi)外壓系數(shù)功率譜函數(shù)則有很大的區(qū)別，尤其是在內(nèi)壓功率譜函數(shù)曲線上，出現(xiàn)了2 個共振峰，內(nèi)壓系數(shù)功率譜在這2個共振峰的頻率附近明顯大于外壓系數(shù)功率譜。根據(jù)相關研究的結(jié)論，這2 個共振峰可能是由于斜風剪切流在孔口處產(chǎn)生的漩渦脫落與內(nèi)壓Helmholtz共振導致的。
圖7 孔口周邊外壓系數(shù)與內(nèi)壓系數(shù)功率譜對比
圖8 內(nèi)壓系數(shù)導納曲線
為了更直觀地對比內(nèi)外風壓系數(shù)功率譜函數(shù)的大小，將內(nèi)壓系數(shù)功率譜與外壓系數(shù)功率譜的比值定義為內(nèi)壓導納，并將0°正風向角和70°斜風向角的開孔模型內(nèi)壓導納曲線繪于圖8中。由圖可見，在0°正風向角下，總體上內(nèi)外壓系數(shù)功率譜相差不大，內(nèi)壓導納的值基本在1左右的范圍內(nèi)波動。但是在70°斜風向角下，內(nèi)壓導納波動很大并且也出現(xiàn)了2 個明顯的峰值，其最大值甚至超過了10，這意味著內(nèi)壓系數(shù)功率譜可能比外壓系數(shù)功率譜大10 倍以上。因此，開孔結(jié)構(gòu)的內(nèi)壓響應在功率譜上與孔口周邊外壓在斜風向角下存在很大的差異。
3.2 湍流強度的影響
建筑物周邊不同類型的地貌會導致來流湍流強度發(fā)生很大的變化，對開孔建筑的內(nèi)壓響應具有不可忽略的影響。圖9 為來流風速為10 m·s-1時，邊長53 mm 的方形中心開孔結(jié)構(gòu)在3 種不同湍流度的均勻湍流風場中平均內(nèi)壓系數(shù)隨風向角的變化曲線。
展開 隨機振動中Grms值的計算
功率譜密度 (PSD)：功率譜密度指，隨機信號的各個頻率分量所包含的功率（或稱能量）在頻域上是怎樣分布的，通常用PSD表示，單位為g2/Hz。它在頻域上分布的曲線圖稱譜圖（簡稱譜）。橫坐標為頻率，縱坐標為功率譜密度g2/Hz（稱功率譜）。
加速度譜密度 (ASD)：表示隨機信號的各個頻率分量所包含的加速度方均值在頻域上是怎樣分布的。通常用ASD表示，單位：m2/s3或是(m/s2)2/Hz。
加速度譜密度與功率譜密度換算關系：
三
計算Grms
多邊形譜
上圖中含有平直譜 (A1)，梯形譜 (A1+A3)，上升譜 (A2)，下降譜 (A3)。
1. 平值譜的計算，PSD值乘以頻率之差得到A1 的面積。
2. 上升譜面積A2計算公式
式中，功率譜密度：PSDj =0.1；初始頻率：fi =10；末端頻率：fj =100；斜率：m =3。
3. 下降譜面積A3計算公式
式中，功率譜密度：PSDi =0.1；初始頻率：fi =1000；末端頻率：fj =2000；斜率：m=6。
總面積：
上述值經(jīng)過開方可以得到：
這幾個公式可以在EXCEL表中輸入，在需要用的時候只要將相關數(shù)據(jù)輸入即可直接得到。
參考資料：
鄭志國、王宇峰 2009《隨機振動中的參數(shù)介紹及計算方法》
展開 聯(lián)合ABAQUS與Fe-safe的隨機振動疲勞分析（隨機疲勞理論及有限元軟件操作講解） ￥25
2.1 隨機過程基礎理論與功率譜密度PSD
隨機過程按統(tǒng)計特性隨時間變化情況，可以分為兩類: 平穩(wěn)隨機振動和非平穩(wěn)隨機振動。 為了描述平穩(wěn)過程隨時間變化的特性，一般采用相關函數(shù)。而功率譜密度可以在頻域范圍內(nèi)描述平穩(wěn)過程，能夠表達隨機過程的能量分布情況，是隨機過程在頻域內(nèi)的重要參數(shù)。
自相關函數(shù)Rx(τ)和功率譜密度Sx(ω)符合維納一辛欽定理:
為了方便工程中的運算，根據(jù) ω = 2πf，把圓頻率換算為頻率f。
2.2 應力功率譜密度
結(jié)構(gòu)激勵信號的功率譜密度為W( f )，則應力功率譜密度為：
其中，H ( f )為應力頻響函數(shù)。
2.3 應力功率譜密度慣性矩與不規(guī)則因子
在時域中，常用一些統(tǒng)計參數(shù)來描述一個隨機應力應變時間歷程中1s的樣本，如圖2所示，圖中E[0]為樣本中自下而上穿越均值的次數(shù)，E[P]為樣本中出現(xiàn)峰值的次數(shù)。
E[0]和E[P]可以通過功率譜密度函數(shù)的i階慣性矩mi換算得到。慣性矩即為功率譜密度函數(shù)曲線下包括的面積。如圖3所示。
應力功率譜密度的第 i 階慣性矩 mi 為
不規(guī)則因子計算公式為
不規(guī)則因子γ，來分辨平穩(wěn)過程X(t) 是符合窄帶過程還是寬帶過程。當不規(guī)則因子γ接近0時，平穩(wěn)過程X(t) 符合寬帶隨機過程。當不規(guī)則因子γ趨向于1時，平穩(wěn)過程X(t)符合窄帶隨機過程。
展開 HyperWorks仿真——隨機振動理論簡介
對隨機信號來說，當時，自相關函數(shù)即為時域上的平均功率；當非常小時，和是相差很小的，當慢慢增大時，和逐漸不再相關，自相關函數(shù)快速衰減為零。
當隨機信號中具有某種周期性信號時，其自相關函數(shù)會將其識別出來，并且忽略其相位信息。
根據(jù)維納-辛欽定理，信號的自相關函數(shù)與功率譜密度是一對傅里葉變換對：
所以隨機信號的功率譜密度可以通過對其自相關函數(shù)的傅里葉變換求得。
注：這里有一個問題，前面說隨機信號不一定滿足傅里葉變換的條件，后面說對信號進行自相關處理后就可以進行傅里葉變換進而求得功率譜密度，但是對功率譜密度的定義是對隨機信號進行傅里葉變換得到的。也就是功率譜密度是對一個不一定滿足變換條件的函數(shù)進行傅里葉變換得到的。這個地方我還不是很理解。
總結(jié)一下，隨機振動的信號輸入是隨機的，不能作為振動分析的輸入，但是其分布是滿足正態(tài)分布的，可以作為隨機振動分析的輸入。為了確定隨機振動信號的分布，需要知道信號的功率譜密度，這時可以通過求隨機信號的自相關函數(shù)，對其進行傅里葉變換來得到，進而得到了隨機振動的輸入。
定義好隨機振動的輸入后，我們需要看下隨機振動的分析結(jié)果。如前面所說，隨機振動分析的輸出是按照概率分布來的。得到的應力符合正態(tài)分布。如下圖所示，求出的結(jié)果為1σ應力，表明應力值小于這個值的概率為68.27%，通常情況下會選取3σ應力作為結(jié)果，其值大于可能出現(xiàn)的99.73%的應力。
二、參考文獻
如何理解隨機振動的功率譜密度. J Pan. 知乎
《OptiStruct結(jié)構(gòu)分析與工程應用》，劉勇，陳斌，羅峰. 機械工業(yè)出版社
文章來源：結(jié)構(gòu)仿真學習
展開 基于頻率響應的電池包隨機振動仿真
頻率響應的設置
隨機振動的設置
（1）RANDPS卡片：用于定義功率譜度的放大因子。
本例的動力電池包通過剛性工裝連接在振動試驗臺上，默認四個支撐托腳受到振動臺的激勵都相同，所以此電池包的隨機振動是一個自功率譜密度加載過程。RANDPS卡片如下圖所示。
其中：1處J（1）設置激勵的工況；2處的K（1）設置激勵的工況；由于該動力電池包為自功率譜密度加載，所以J（1）= K（1）。3和4處X（1）、Y（1）設置功率放大系數(shù)，功率放大因子由X和Y控制，如下式所示：
Sjk（F）=（X+iY）G（F）
本例采用MM-MPa-N單位制，設置時已經(jīng)進行了手動轉(zhuǎn)化，故而功率譜密度放大因子X+iY=1,所以X=1，Y=0。
圖中5處的TID（1）設置功率譜密度對應的卡片。
（2）TABRAND1卡片：用于定義加速度功率譜密度。
輸入動力電池包在垂直路面方向上隨機振動加速度功率譜密度。加速度功率譜密度的輸入有三種方式：
1）定義曲線的轉(zhuǎn)折點。
2)利用EXCEL將載荷譜編制成兩列，保存為CSV格式，然后通過Utility選項卡下的TABLE Create載入。此方法適合載荷譜數(shù)據(jù)量較大的情況。
3)直接在通過上圖中的3處創(chuàng)建載荷譜。
（3）TABDAMP1卡片：定義模態(tài)阻尼。
阻尼對頻率響應的峰值影響很大，鑒于條件有限，本例將各階模態(tài)阻尼取為0.1；讀者也可以通過PARAM下的G_V1定義固定的結(jié)構(gòu)阻尼。
展開 基于nCodeDesignlife的電池箱疲勞壽命研究
摘要：針對電池箱對振動疲勞耐久性能的要求，結(jié)合國標中隨機振動的加速度功率譜密度函數(shù)和材料的S-N曲線，采用Miner線性累積損傷理論和Dirlik疲勞壽命計算方法，對電池箱進行隨機振動條件下的疲勞壽命分析。
關鍵詞：電池箱；疲勞壽命;隨機振動；nCodeDesignlife
路面不平度產(chǎn)生的隨機振動是造成電動汽車零部件發(fā)生疲勞破壞的主要因素[1]。電池箱對電池組起防護和保障安全的作用，其疲勞耐久性能對于保障車輛及乘員安全至關重要。
電池箱的疲勞壽命分析主要有時域和頻域方法。時域法采用經(jīng)典的雨流循環(huán)計數(shù)統(tǒng)計載荷信息，容易丟失載荷數(shù)據(jù)，計算量大，工程應用不是很廣泛;頻域法從概率統(tǒng)計的角度統(tǒng)計載荷信息，采用功率譜密度（PSD）描述隨機振動載荷在各個頻率成分上的統(tǒng)計特性[2]，廣泛應用于航天、海工、汽車等領域。目前電池箱的疲勞壽命研究[3-5]大多致力于時域振動和定頻振動，隨機振動的疲勞壽命研究較少。本文基于頻率響應分析研究電池箱隨機振動的疲勞壽命，為電池箱的疲勞壽命分析提供一種高效的方法。
1 隨機振動疲勞壽命分析方法
本文參照電池箱振動測試的國標選取加速度功率譜密度，避免建立整車模型和提取加速度載荷譜的復雜過程。
1.1 加速度功率譜密度
隨機振動無法用確定的函數(shù)關系式表示，只能通過概率統(tǒng)計的方法表示。在頻域內(nèi)，采用功率譜密度函數(shù)表示隨機振動在各個頻率的統(tǒng)計特性[6]。功率譜密度函數(shù)Sx(ω)為自相關函數(shù)R(τ)的傅里葉變換公式為：
GB/T 31467.3-2015[7]中規(guī)定電池箱隨機振動的加速度功率譜密度如圖1所示。
圖1 加速度功率譜密度
1.2 疲勞壽命計算方法
Dirlik計算方法是疲勞仿真軟件nCodeDesignlife所采用的方法。
展開 matlab如何做頻譜分析！
感覺matlab做的都是幅值的平方，用FFT做，應該算功率譜把，有的取模，有的取模方。
不知道有沒有能做出信號頻譜的程序，如取信號 x=sin(2*pi*60*t)+1/2*sin(2*pi*120*t)
應該包含兩個頻率成分一個是60Hz幅值1，一個是120Hz幅值1/2,我用matlab編的頻率成分對，就是幅值不對，幅值取的應該是功率譜把！ 對于功率譜，功率譜密度，還有它的單位，老是搞不清楚，也看到不少解釋的文章，希望大家不厭其煩在給我解釋一下，也希望能把我說的這個例子用matlab實現(xiàn)一下，主要幅值我編的，老是不是頻譜，而用我就一般采用origin做頻譜了，matlab的信號工具箱，也一樣頻率很準，幅值也不知道是什么單位！ 我總結(jié)一下matlab畫功率譜的函數(shù)不下4個，希望有關專家給總結(jié)一下在解釋一下！
展開 
abaqus隨機響應分析
加速度的單位是：米/秒平方
加速度均方值的單位為：米平方/秒四次方，若用G，則為G平方
均方根加速度的單位：米/秒平方，若用G
加速度功率譜密度的單位：米平方/秒三次方，若用G，則為（G平方/Hz）
若對加速度功率譜密度曲線在頻域內(nèi)進行積分，積分的結(jié)果就是PSD曲線下方的面積，單位為米平方/秒四次方，就是加速度均方值。其表征了振動的強度。其開方的單位為（米/秒平方）：為均方根加速度值。(對于均值，均方值，均方根值的概念，以及傅里葉變換可以參考機械工程測試基礎。)
隨機響應分析的輸入是PSD，輸出為帶前綴的“R”，RMS為均方根應力，RA為均方根加速度。還可以在歷史輸出中輸出某些SET集內(nèi)節(jié)點的加速度PSD曲線和位移psd曲線。
三、 隨機響應分析實例：
1、 模型描述
2、 網(wǎng)格劃分及模態(tài)：
（1）、厚度方向放一層單元，使用一階單元C3D8R，1階模態(tài)的頻率為8。8691HZ。
（2）、厚度方向放一層單元，使用二階單元C3D20R，1階模態(tài)的頻率為85。512 HZ。
（3）、厚度方向放二層單元，使用一階單元C3D8R，1階模態(tài)的頻率為74。149HZ。
在第三種情況下，通過增加厚度方向的網(wǎng)格層數(shù)，第一階模態(tài)的頻率會接近于第二種情況的結(jié)果。在通過上面的對比可以看出，在對薄板型的結(jié)構(gòu)進行模態(tài)分析時，厚度方向要多放幾層單元或者使用二階單元才能獲得比較可信的結(jié)果。
3、 分析結(jié)果：
從圖8和圖1的對比可以看出，激勵的加速度功率譜密度曲線與某節(jié)點的響應的加速度功率譜密度曲線具有相似的變化趨勢。
ABAQUS隨機響應分析.pdf
展開 隨機振動分析-abaqus（附一個電池包計算案例） ￥20
目錄
一、隨機振動的定義、特點及常見場景
二、隨機振動的數(shù)學特征--正態(tài)分布
三、 隨機振動信號為什么要用功率譜密度（PSD）表達？
四、如何將時域隨機振動曲線轉(zhuǎn)換得到功率譜密度曲線
五、 隨機振動分析理論
附.常見功率譜密度曲線給出形式
附.以dB/oct形式給出的功率譜密度曲線如何計算
附.國標中定義的PSD譜總均方根加速度值是如何計算的？
六. 隨機振動分析案例-abaqus
第一步：計算結(jié)構(gòu)模態(tài)，輸出位移和應力。
第二步：隨機振動分析
2.1 定義輸出頻率上下限和模態(tài)阻尼
2.2 定義PSD載荷及加載
2.3 定義輸出
2.4 隨機振動計算頭文件設置
2.5 隨機振動分析結(jié)果
2.6 隨機振動σ應力結(jié)果評價
展開 模態(tài)應力、頻響應力和PSD應力
隨機激勵載荷是頻域信號表征的，通常采用功率譜密度函數(shù)（PSD）表示。隨機激勵下的響應是統(tǒng)計意義下描述的響應，在任何瞬時響應具體大小未知，但其大小的概率可以給出。
隨機響應分析通常是作為頻率響應分析的后處理。首先施加單位正弦激勵，得到頻響函數(shù)，然后根據(jù)用戶給定的載荷條件（形式為各載荷的自功率譜密度和互功率譜密度）。輸出為響應的功率譜密度、響應的自相關函數(shù)或響應的均方值，其中應力響應的功率譜密度我們簡稱為PSD應力。
與頻響函數(shù)一樣，響應的PSD也是以頻率為橫坐標的函數(shù)，但PSD是正的實函數(shù)，并沒有相位的概念。
一個線性系統(tǒng)，承受平穩(wěn)隨機激勵載荷x(t)，其PSD為Sx(ω)，響應信號為y(t)，其PSD為Sy(ω)，則可按下式計算y(t)的功率譜密度，
其中H(ω)為載荷x(t)所對應的頻響函數(shù)。
如果有多個隨機載荷信號x1(t),x2(t),…,xN(t)共同作用，則
其中Sxr,xs(ω)當r≠s時，為載荷xr(t)和xs(t)的互功率譜密度函數(shù)（CPSD）；當r=s時，為載荷xr(t)的功率譜密度函數(shù)。
對于互不相關的多個載荷信號，其互功率譜密度全部為零，則
所以，只要有各載荷的頻響函數(shù)Hr(ω)，我們就可以根據(jù)激勵載荷的PSD矩陣來求得結(jié)構(gòu)上各點的PSD應力。
展開 干貨｜漫話噪聲系數(shù)
當外部電阻與被測器件的輸入阻抗匹配時，器件或系統(tǒng)得到的輸入噪聲功率為
通常取室溫 T=T0=290K，此時這個電阻的噪聲功率譜密度為
在此條件下，噪聲系數(shù)為
通常將上述式子中的 PNO(total)/B 作為輸出端的噪聲功率譜密度，記為 DNO，則有
改用分貝的形式，就是
2、利用噪聲指數(shù)計算系統(tǒng)的輸出噪聲
根據(jù)噪聲指數(shù)的定義，系統(tǒng)輸出端的噪聲功率譜密度為
然而值得注意的是：上述公式僅在信號源的噪聲僅為匹配的信號源內(nèi)阻的熱噪聲條件下成立，或者說信 號源的噪聲功率譜密度被規(guī)定為
這個規(guī)定對于生產(chǎn)廠標定噪聲系數(shù) 是必須的，但在計算一個實際的射頻系統(tǒng)的輸出噪聲時，由于輸入系統(tǒng)的噪聲功率譜密度未必就是 kT0，所 以很多時候不能直接套用 14 式計算射頻系統(tǒng)的輸出噪聲功率，也不能直接套用 4 式計算射頻系統(tǒng)的輸出信 噪比。
一個噪聲指數(shù)為 NF(dB)的系統(tǒng)，其增益為 GP(dB)，信號源輸入的噪聲功率譜密度為 DNS(mW/Hz)，則 此系統(tǒng)內(nèi)部噪聲在輸出端的功率譜密度為
信號源輸入的噪聲在輸出端的功率譜密度為
或者寫成分貝的形式：
3、等效噪聲溫度
在低噪聲系統(tǒng)中，通常用等效噪聲溫度 Te 替代噪聲系數(shù)。
展開 