扫频输入信号 matlab,扫频信号PSD分析

本帖最后由 Robinsonw 于 2016-1-13 21:47 编辑

大家好，我现在有2个扫频信号，一个输入a1Fsw，一个输出a1Fmf，然后希望通过PSD分析来确定它们的传递函数，采样频率1000Hz

%%主程序

a1=importdata('C:\XXXXXXX\CANOE000129.mat');读取实验结果

a1Fmf_Spannung=a1.Kraft_Signal;  %输出力(做实验的时候跟输入位移a1Weg有关)

a1Fsw_Spannung=a1.Fz_HL;         %输入力

a1Weg_Spannung=a1.Weg_Signal;    %输入位移(正弦线性扫频)

a1t=a1.Zeit_absolut; %数据采集时间变化

a1Weg=10./(2^16).*5.*a1Weg_Spannung*2; %把电压信号变成真实值

Punktnull1=mean(a1Weg(2:6000,1)); %计算出扫频信号和零点偏移值

a1WegNeu=a1Weg-Punktnull1;        % 把位移信号移动到零点

Anfangszahl=1113; %读出扫频信号起始点

Endzahl=18432; %读出扫频信号结束点

Input_tf=a1Fsw; %a1Fsw作为输入

Output_tf=a1Fmf;%a1Fmf作为输出

Nfft=8192;

t_start=1;

overlap=0.8;

window='hann';

x1=Input_tf;

x2=Outp