How to find wavelength or wavenumber from signal? (MATLAB)

40 次查看（过去 30 天）

显示更早的评论

Michal G 2015-9-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/243825-how-to-find-wavelength-or-wavenumber-from-signal-matlab

回答： mm99 2019-12-27

signal.txt

在 MATLAB Online 中打开

Hello

I am working on my master thesis and it is connected with signal processing.

I have a numerical model of plate in which i have got an extortion as a 5 periodic sin wave.

I have collected displacements in Y direction from all points on surface of my plate (400 points) in one time increment. I put it into Matlab as a vector:

data = load('jednachw100.txt');  
data=reshape(data,[],1);
y=fft2(data);
plot(abs(fft(data)));

and i have recived something like this:

In x axis i have my all sampling points which is 400(all nodes from surface of plate). My task is to get a wavenumber in x-axis but i can't do this.

Thaks for help in advance.

2 个评论
显示无隐藏无

Image Analyst 2015-9-18

What's your definition of wavenumber?

Michal G 2015-9-20

definition of wavenumber is k=2pi/wavelength

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Walter Roberson 2015-9-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/243825-how-to-find-wavelength-or-wavenumber-from-signal-matlab#answer_192914

在 MATLAB Online 中打开

Untested

numpoint = length(data);
curtick = get(gca, 'XTick');
lefthalf = curtick(curtick <= numpoint/2);
leftwaveL = numpoint ./ (lefthalf - 1);
righthalf = curtick(curtick > numpoint/2);
rightwaveL = numpoint ./ (numpoint - righthalf + 1);
waveL = [leftwaveL(:); rightwaveL(:)];
ticklab = cellstr(str2num(waveL));
set(gca, 'xticklabel', ticklab);

20 个评论
显示 18更早的评论隐藏 18更早的评论

Walter Roberson 2015-9-20

I do not see the txt file attached. After you select a file, you have to click to Attach File.

Michal G 2015-9-20

I can't add any more files today:( i sent yopu a private message with my vector.

请先登录，再进行评论。

mm99 2019-12-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/243825-how-to-find-wavelength-or-wavenumber-from-signal-matlab#answer_407910

Hi!

I am doing similar thing, but have some trouble with sine extortion. I get an error: "Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is 1-by-101." This 1-by-1 is the point in which I want the sine to start and 1-by-101 is this sine with hann window. Do you maybe know what is wrong ?

Thanks in advance!