How to create two graph plots in one, plus scale aside?

조회 수: 2 (최근 30일)
cbentter
cbentter 2017년 5월 4일
편집: cbentter 2017년 5월 13일
How to do this figure in MatLab?
Partial code:
M = 256; sigma = 6.0;
h = sin(+2*pi*sigma*log((1-([0:M-1]/M))-0.2));
h(1:37) = 0; h(193:M) = 0;
H= [fft(h) fft(h) fft(h)];
k = ([0:floor(M/2)-1 -floor(M/2):-1]);
S(1,:)=sum(ifft(H(1:M)))/M*ones(1,M);
for n=1:floor(M/2)-1;
W=(1./M)*exp(-n^2*k.^2/(2*M^2))...
.*exp(-2*pi*i*(n*k/M-sign(n) ...
*sigma*log(sigma+abs(n)*k/M)));
W(M/2+1:M-ceil(sigma*M/abs(n))+1)=0;
W = W/sum(W);
S(n+1,:) = ifft(H(n+M+1:n+M+M).* fft(W));
end
figure(10), contour(abs(S),'ShowText','off');
colorbar
Source code and picture, article DOI: 10.1016/j.sigpro.2004.03.015
Title: "Time-local Fourier analysis with a scalable, phase-modulated analyzing function: the S-transform with a complex window"

채택된 답변

Santhana Raj
Santhana Raj 2017년 5월 4일
The trick is to use position property in subplot.
figure, subplot('position',[0.1 0.3 0.8 0.7]),contour(abs(S),'ShowText','off');xlim([0 M]);
subplot('position',[0.1 0.1 0.8 0.1]),plot(1:M,h),xlim([0 M]);
you can modify the parameters of position vector to get the exact to what you want.

추가 답변 (1개)

cbentter
cbentter 2017년 5월 4일
Similar figure and complete soluction basead on 'Santhana Raj' answer
len = 512;
t(1:len) = [0:0.25:127.75];
h(1:128) = 0;
% h = zeros(1,512);
h(129:512) = exp(-4*[0:383]/256).*sin(2*pi*[0:383]*20.4/512);
h(157:512) = h(157:512)+ exp(-5*[0:355]/256).*sin(2*pi*[0:355]*30.7/512);
h(269:512) = h(269:512) + exp(-4*[0:243]/256).*sin(2*pi*[0:243]*25.3/512);
h(397:512) = h(397:512) + exp(-4*[0:115]/256).*sin(2*pi*[0:115]*15.6/512);
figure(10), plot(h), xlim([0 len]);
xlabel('Sample Point Number'), ylabel('Amplitude');
M = 512; H = [fft(h) fft(h)];
g = 1.0; gf = 0.5;
gb = ((pi-4)*gf+sqrt((8*pi+3*pi^2)*gf- ...
(4*pi^2-8*pi)*g))/2/(pi-2);
t = [[0:floor(M/2)-1]/gb [-floor(M/2):-1]/gf];
W = [1 zeros(1,M-1)];
for m=0:floor(M/4)-1;
STR(m+1,:) = ifft(H(m+1:m+M) .* W);
w = abs(m+1)/sqrt(2.*pi)*2/(gb+gf) ...
* exp(-(m+1)^2*t.^2/(2*M^2));
W = fft(w/sum(w));
end
figure(11), contourf(abs(STR),'ShowText','off');
ylabel('Frequency (Hz)');
% shading interp;
colormap(gray)
colormap(flipud(colormap))
shading flat;
colorbar
figure(12)
% subplot(2,1,1),
subplot('position',[0.1 0.45 0.77 0.5]),
contourf(abs(STR),'ShowText','off');
colormap(gray)
colormap(flipud(colormap))
shading flat;
colorbar
colorbar('position',[0.9 0.15 0.032 0.8]);
% subplot(2,1,2),
subplot('position',[0.1 0.15 0.77 0.2]),
plot(h), xlim([0 len]);
xlabel('Sample Point Number'), ylabel('Amplitude');
reference: PII. S1064827500369803
  댓글 수: 2
KSSV
KSSV 2017년 5월 5일
@cbentter, it is better to give credits to the user who have helped you. I prefer accepting Santha Raj's answer rather then own-self.
cbentter
cbentter 2017년 5월 13일
편집: cbentter 2017년 5월 13일
See below Santhana Raj comment.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 2-D and 3-D Plots에 대해 자세히 알아보기

제품

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by