Gaussian-modulated function for Temperature definition

Question

Elisa Revello 2023-3-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1930655-gaussian-modulated-function-for-temperature-definition

回答： Ayush Aniket 2023-8-23

Twall_gaussian function.png

I am trying to implement a code that describes the temperature distribution at the wall (T_wall in input) with a function of the type shown in the attachment (first section equal to Tamb, central section with a Gaussian and final section again with Tamb).

I don't understand how to modify the syntax reported for the gauspuls command with the correct frequency and bandwidth, as showed in the following code.

Tamb = 293;
Tmax = Tamb+200; % peak amplitude for temperature distribution
% t = linspace(0,300); time interval 5 min
t0 = Tamb*ones(1,19);
t1 = [20:100]; % time interval for gaussian T_wall
t2 = Tamb*ones(1,200);
t = [t0 t1 t2]; % time interval definition 
%tc = gauspuls('cutoff',50e3,0.6); 
%t = -tc : 1e-7 : tc; 
%[yi,yq,ye] = gauspuls(t,50e3,0.6); 
%plot(t,ye)
%legend('Envelope');
[yi,yq,ye] = gauspuls(t,50e3,0.6); 
figure(1)
plot(t,ye)
xlabel('t [s]')
ylabel('Wall Temperature [K]')
title('Wall Temperature as a gaussian function')
legend('Wall Temperature')

Thanks in advance, if someone can help me.

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

请先登录，再进行评论。

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

Answer 1

Ayush Aniket 2023-8-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1930655-gaussian-modulated-function-for-temperature-definition#answer_1292022

在 MATLAB Online 中打开

Hi Elisa Revello,

I understand you are trying to plot the graph as attached in your question. For plotting the sinusoidal part of the plot you are trying to use ‘gauspuls’ function of MATLAB. However, based on the image, I think it will not be suitable for your case as ‘gauspuls’ command plots a unit-amplitude Gaussian-modulated sinusoidal RF pulse which will have decaying amplitude for each of period of the sine wave. In the image the sinusoidal portion has 4 pulses with equal amplitude. Hence, I suggest you to use a squared sinusoidal wave for your plot as shown below:

Tamb = 293; % Ambient temperature
Tmax = Tamb + 200; % Maximum temperature at the peak
t0 = 1:19; % Initial section with ambient temperature
t1 = 20:100; % Time interval for the Gaussian section
t2 = 101:200; % Final section with ambient temperature
t = [t0 t1 t2]
F = 1 / 40 % as each sine wave has a period of 40s(before squaring)
y = Tamb + (Tmax -Tamb) * sin(2 * pi *F*t ).^2
Y = [Tamb*ones(1,numel(t0)), y(t1), Tamb*ones(1,numel(t2))]
plot(t,Y)
xlabel('t[s]')
ylabel('Wall Temperature [K]')
yticks([Tamb,Tmax]);
yticklabels({'Tamb','Tmax'});
xticks([20,100]);
xticklabels({'20','100'});
title('Wall Temperature as a gaussian function')
legend('Wall Temperature')