visualise fft plot in 3D

Question

Athira Surendran 2017-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333539-visualise-fft-plot-in-3d

评论： Daniel Esteban Caballeros Tejada 2021-2-10

采纳的回答： Star Strider

在 MATLAB Online 中打开

Hi,

I want to plot all the fft plots in a single figure. something like the attached image

Here's my code

for i=1:m
    [f(i,:) mag(i,:)]=fftplot(2000000,corr_amp(i,:));
    subplot(3,4,i)
    plot(f(i,:),mag(i,:))
    axis([25000,85000,0,1.1]);
    strtime = ['Frequency Spectrum at \delta=', num2str(load(1,i)),'\mum'];
    title(strtime,'fontsize',10)
    xlabel('Frequency (Hz)')
    ylabel('Magnitude')
    grid on
    grid minor
end;

this code gives me all fft plots as separate plots in a single figure, but i want to arrange all the fft plots in 3D (third axes is 'load' variable in the .mat file attached) as shown in image to see the variation accurately.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2017-4-4

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333539-visualise-fft-plot-in-3d#answer_261713

在 MATLAB Online 中打开

Use the ribbon (link) plot.

The Code —

d = load('Athira Surendran sjs1.mat');
amp = d.amp';                                           % Amplitude (Transposed)
load = d.load;                                          % Load
t = d.t;                                                % Time
Ts = mean(diff(t));                                     % Sampling Interval
Fs = 1/Ts;                                              % Sampling Frequency
Fn = Fs/2;                                              % Nyquist Frequency
L = length(t);
FTamp = fft(amp)/L;                                     % Fourier Transform
Fv = linspace(0, 1, fix(L/2)+1)*Fn;                     % Frequency Vector
Iv = 1:length(Fv);                                      % Index Vector
figure(1)
hr = ribbon(Fv, abs(FTamp(Iv,:))*2, 0.1);
grid on
for k1 = 1:length(hr)
    set(hr(k1), 'FaceColor','b', 'EdgeColor','b')
    xpos = get(hr(k1), 'XData');
    xtix(k1) = mean(xpos(1,:));
end
set(gca, 'XTick',xtix, 'XTickLabel',load)
xlabel('Load')
ylabel('Frequency (Hz)')
zlabel('Amplitude')
view([135,  30])

I believe that I chose the correct variables. If not, changing my code to use the correct ones is straightforward.

The Plot —

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

agung pratama 2020-8-15

编辑：agung pratama 2020-8-15

Can i transform a 3D surface that x,y,z axis in pixel into mm?

请先登录，再进行评论。

Answer 2

KSSV 2017-4-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333539-visualise-fft-plot-in-3d#answer_261685

在 MATLAB Online 中打开

You may follow something like this:

x = linspace(0,2*pi) ;
N = 10 ;
pos = 1:N ;
figure
hold on
for i = 1:N
    z = i*sin(x) ;
    y = repmat(pos(i),size(x)) ;
    plot3(x,y,z,'r');
end
view(3)