| 查看: 438 | 回复: 0 | |||
hailiang金虫 (正式写手)
|
[求助]
请教，对信号做傅里叶分析，当原始信号不是整数或加0.5时结果不对
|
|
请教，对信号做傅里叶分析，当原始信号不是整数或加0.5时结果不对 我用下面的程序对余弦信号做频谱分析，当信号频率是整数或者是正负0.5的时候都正确，就是求出的振幅和频率刚好是原始信号一样，但是如果原始信号频率不是整数或正负0.5时就不正确，请问是精度设置的不够高，还是采样频率和采样点数不对！ 下面的程序是Mathematica程序，Matlab程序试了也一样 A1 = 3; (*% 频率F1信号的幅度*) (*F1=1/(2 Pi); (*% 信号1频率 (Hz)*)*) F1 = 0.11; Fs = 256; (*% 采样频率 (Hz)，要大于信号频率的两倍*) P1 = 100; (*% 信号1相位 (度)*) NN = 2048; (*% 采样点数，通常取2的整数次方*) S = A1*Cos[2*Pi*F1*t + Pi*P1/180]; (*输入原始信号*) Plot[S, {t, 0, 10}, PlotPoints -> 200, PlotLabel -> "原始信号"] Plot[S, {t, 0, NN/Fs}, PlotPoints -> 200, PlotLabel -> "原始信号"] 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模值结果*) 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; %先关闭所有图片 Adc=2; %直流分量幅度 A1=3; %频率F1信号的幅度 F1=0.1; %信号1频率(Hz) Fs=256; %采样频率(Hz) P1=-30; %信号1相位(度) N=1024; %采样点数 t=[0:1/Fs:N/Fs]; %采样时刻 %信号 S=A1*cos(2*pi*F1*t+pi*P1/180); %显示原始信号 plot(S); title('原始信号'); figure; Y = fft(S,N); %做FFT变换 Ayy = (abs(Y)); %取模 plot(Ayy(1:N)); %显示原始的FFT模值结果 title('FFT 模值'); figure; Ayy=Ayy/(N/2); %换算成实际的幅度 Ayy(1)=Ayy(1)/2; F=([1:N]-1)*Fs/N; %换算成实际的频率值 plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果 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)); %计算相位 Pyy(i)=Pyy(i)*180/pi; %换算为角度 end; plot(F(1:N/2),Pyy(1:N/2)); %显示相位图 title('相位-频率曲线图'); |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有79人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有13人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
请教一下JPSP的申请
已经有4人回复
- 新手提RNA，请教几个问题
已经有8人回复
- 陶瓷信回复请教
已经有10人回复
- 有在澳洲悉尼的吗，请教个问题哦
已经有3人回复
- 请教一道物化题
已经有7人回复
- 请教 ĸ ŋ配位模式的区别
已经有3人回复
- 求南京理工2010和2011年《信号系统与数字电路》考研真题及答案
已经有4人回复
- 没有信号肽会影响蛋白的表达过程
已经有9人回复
- 信号与线性系统分析课件 第四版
已经有44人回复
- 光学信号cos(sin(2pi*f*t))的傅里叶展开？
已经有6人回复
- 信号处理方面世界上比较权威的杂志有哪些，麻烦大家略微介绍一下这些杂志
已经有5人回复
- 材料拉曼信号的数据库
已经有3人回复
- 粉末xrd数据，jade分析,txt比raw的信号强是为什么
已经有11人回复
- 【求助】自制的Faraday cage(法拉第笼)测试时, 屏蔽不了手机信号.求好心人帮忙解惑...
已经有6人回复
- 【求助/交流】测序信号弱
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
25