I succesfully plotted my FFT with MATLAB discussion help. Now I could not remove the DC component at 0Hz. Which shows me a very high amplitude. Can any one suggest me an idea?
data1 = xlsread(‘Reading 1.xlsx’) ; %Loading Sensor data from Excel file
t = data1 (1:512,2); %Selecting Time vector
s = data1 (1:512,3); %Selecting Z axis vibrations
L = numel(t); %Signal length
Ts = mean(diff(t)); %Sampling interval
Fs = 1/Ts; %Sampling frequency
Fn = Fs/2; %Nyquist frequency
FTs = fft(s)/L; %Fast fourier transform (s- data)
Fv = linspace(0,1, fix(L/2)+1)*Fn; %Frequency vector
Iv = 1:numel(Fv); %Index vector

subplot(2, 1, 1); %plotting top pane
plot(t,s); %Acceleration vs time
set(gca,’xlim’,[1 50]); %Scale to fit
grid; %Grids on
title (‘Acceleration vs time’);
xlabel(‘time(s)’);
ylabel(‘Acceleration’);

subplot(2, 1, 2); %Plotting bottom pane
plot(Fv, abs(FTs(Iv))*2,’red’); %FFT – Amplitude vs Frequency
grid
title (‘Fast fourier transform’);
xlabel(‘Frequency (Hz)’);
ylabel (‘Amplitude (m)’);I succesfully plotted my FFT with MATLAB discussion help. Now I could not remove the DC component at 0Hz. Which shows me a very high amplitude. Can any one suggest me an idea?
data1 = xlsread(‘Reading 1.xlsx’) ; %Loading Sensor data from Excel file
t = data1 (1:512,2); %Selecting Time vector
s = data1 (1:512,3); %Selecting Z axis vibrations
L = numel(t); %Signal length
Ts = mean(diff(t)); %Sampling interval
Fs = 1/Ts; %Sampling frequency
Fn = Fs/2; %Nyquist frequency
FTs = fft(s)/L; %Fast fourier transform (s- data)
Fv = linspace(0,1, fix(L/2)+1)*Fn; %Frequency vector
Iv = 1:numel(Fv); %Index vector

subplot(2, 1, 1); %plotting top pane
plot(t,s); %Acceleration vs time
set(gca,’xlim’,[1 50]); %Scale to fit
grid; %Grids on
title (‘Acceleration vs time’);
xlabel(‘time(s)’);
ylabel(‘Acceleration’);

subplot(2, 1, 2); %Plotting bottom pane
plot(Fv, abs(FTs(Iv))*2,’red’); %FFT – Amplitude vs Frequency
grid
title (‘Fast fourier transform’);
xlabel(‘Frequency (Hz)’);
ylabel (‘Amplitude (m)’); I succesfully plotted my FFT with MATLAB discussion help. Now I could not remove the DC component at 0Hz. Which shows me a very high amplitude. Can any one suggest me an idea?
data1 = xlsread(‘Reading 1.xlsx’) ; %Loading Sensor data from Excel file
t = data1 (1:512,2); %Selecting Time vector
s = data1 (1:512,3); %Selecting Z axis vibrations
L = numel(t); %Signal length
Ts = mean(diff(t)); %Sampling interval
Fs = 1/Ts; %Sampling frequency
Fn = Fs/2; %Nyquist frequency
FTs = fft(s)/L; %Fast fourier transform (s- data)
Fv = linspace(0,1, fix(L/2)+1)*Fn; %Frequency vector
Iv = 1:numel(Fv); %Index vector

subplot(2, 1, 1); %plotting top pane
plot(t,s); %Acceleration vs time
set(gca,’xlim’,[1 50]); %Scale to fit
grid; %Grids on
title (‘Acceleration vs time’);
xlabel(‘time(s)’);
ylabel(‘Acceleration’);

subplot(2, 1, 2); %Plotting bottom pane
plot(Fv, abs(FTs(Iv))*2,’red’); %FFT – Amplitude vs Frequency
grid
title (‘Fast fourier transform’);
xlabel(‘Frequency (Hz)’);
ylabel (‘Amplitude (m)’); dc, fft, amplitude MATLAB Answers — New Questions

​

I find, in R2024b (I haven’t checked other versions), then when running sub-sections of code, for example with , the command line prompt stays visible while the computation is in progress. This is not the case when running normally, and I am wondering if it is a deliberate design choice. If so, I wonder if there is some preference setting that will deactivate this, as I find it makes it a bit more difficult to see when Matlab is busy (in spite of the ‘Busy’ notification in the lower left of the Matlab desktop).I find, in R2024b (I haven’t checked other versions), then when running sub-sections of code, for example with , the command line prompt stays visible while the computation is in progress. This is not the case when running normally, and I am wondering if it is a deliberate design choice. If so, I wonder if there is some preference setting that will deactivate this, as I find it makes it a bit more difficult to see when Matlab is busy (in spite of the ‘Busy’ notification in the lower left of the Matlab desktop). I find, in R2024b (I haven’t checked other versions), then when running sub-sections of code, for example with , the command line prompt stays visible while the computation is in progress. This is not the case when running normally, and I am wondering if it is a deliberate design choice. If so, I wonder if there is some preference setting that will deactivate this, as I find it makes it a bit more difficult to see when Matlab is busy (in spite of the ‘Busy’ notification in the lower left of the Matlab desktop). run section, prompt MATLAB Answers — New Questions

​

I have a large file and try to load it into MATLAB and plot it with timestamp.I have a large file and try to load it into MATLAB and plot it with timestamp. I have a large file and try to load it into MATLAB and plot it with timestamp. timestamp, plot MATLAB Answers — New Questions

​

Hello, i have some data and tried to plot them all in one figure without axes. So far so good but I want all subplots to have one common x-axis. Any ideas?Hello, i have some data and tried to plot them all in one figure without axes. So far so good but I want all subplots to have one common x-axis. Any ideas? Hello, i have some data and tried to plot them all in one figure without axes. So far so good but I want all subplots to have one common x-axis. Any ideas? common xy axes subplot MATLAB Answers — New Questions

​

I have a large file and need to import it into 13 collumns.
Any help?I have a large file and need to import it into 13 collumns.
Any help? I have a large file and need to import it into 13 collumns.
Any help? .csv MATLAB Answers — New Questions

​

The installer for a compiled (using R2015b) app, that has been distributed to many customers for several years, is suddenly (noticed in last couple of days) throwing the above exception. This has happened from multiple sites, domestic and international. The installer is set up to download the MCR from MathWorks.

Is anyone aware of any MathWorks, Windows (10, 11) updates that might be causing this?

Thanks! It’s a bit urgent!!The installer for a compiled (using R2015b) app, that has been distributed to many customers for several years, is suddenly (noticed in last couple of days) throwing the above exception. This has happened from multiple sites, domestic and international. The installer is set up to download the MCR from MathWorks.

Is anyone aware of any MathWorks, Windows (10, 11) updates that might be causing this?

Thanks! It’s a bit urgent!! The installer for a compiled (using R2015b) app, that has been distributed to many customers for several years, is suddenly (noticed in last couple of days) throwing the above exception. This has happened from multiple sites, domestic and international. The installer is set up to download the MCR from MathWorks.

Is anyone aware of any MathWorks, Windows (10, 11) updates that might be causing this?

Thanks! It’s a bit urgent!! http, code 403, download server MATLAB Answers — New Questions

​

Hello I have a big timetable I want to calculate the maximum and minimum values
Any ideas how to proceed?Hello I have a big timetable I want to calculate the maximum and minimum values
Any ideas how to proceed? Hello I have a big timetable I want to calculate the maximum and minimum values
Any ideas how to proceed? max, timestamp MATLAB Answers — New Questions

​

P0 = 101325; % Atmospheric pressure (Pa)
C1 = 10;
C2 = 62.29;
C3 = 62.29;
C4 = 0.557;
C5 = 2391.3;
C6 = 0.609;

omega_i = 0.008; % Inlet humidity ratio (kg/kg dry air)
Tg_i = 303.15; % Inlet air temperature (K)

% Initial conditions
Ts_init = Tg_i;
omega_s_init = 0.008; % solid humidity ratio
omega_o_init = omega_i;
Tg_o_init = Tg_i;

% Sim parameters
N = 180; % Angular resolution
theta = linspace(0, 180, N); % Angular positions in degrees
dt = 1 / (15.5 * 60) / N; % Time step for rotation
tspan = [0, dt * N]; % Simulation time

% ODE system to access constants
function dYdt = desiccant_wheel_ode(t, Y)
% Unpack variables
omega_o = Y(1);
Ts = Y(2);
Tg_o = Y(3);
omega_s = Y(4);

% Saturation pressure and relative humidity
Ps = exp(5294 / Ts – 6); % Saturation pressure (Pa)
phi = (omega_s / 0.622) * P0 / ((0.622 + omega_s) * Ps);

% Empirical functions
s1 = (0.24 / 1.5) * (phi^0.333) * (0.622 * P0) / ((0.622 + omega_s)^2 * Ps);
s2 = (0.24 / 1.5) * (phi^0.333) * (5294 * phi) / (Ts^2);

% Differential equations
d_omega_dt = C1 * (omega_i – omega_o) + C2 * (omega_s – omega_o);
d_Ts_dt = C4 * C5 * (omega_o – omega_s) + C6 * (Tg_o – Ts);
d_Tg_dt = C1 * (Tg_i – Tg_o) + C3 * (Ts – Tg_o);
d_omega_s_dt = (s2 / s1) * d_Ts_dt + (C4 / s1) * (omega_o – omega_s);

% Output derivatives
dYdt = [d_omega_dt; d_Ts_dt; d_Tg_dt; d_omega_s_dt];
end

% ODE ode45
Y0 = [omega_o_init, Ts_init, Tg_o_init, omega_s_init]; % Initial conditions
[t, Y] = ode45(@desiccant_wheel_ode, tspan, Y0);

% Extract results
omega_o = Y(:, 1) * 1000; % Convert to g/kg dry air
Ts = Y(:, 2) – 273.15; % temp to °C

% Plot data
figure;
plot(theta, omega_o, ‘k’, ‘LineWidth’, 2);
xlabel(‘Degree’);
ylabel(‘Humidity [g/kg dry air]’);
title(‘Adsorption-side Outlet Humidity Ratio’);
grid on;P0 = 101325; % Atmospheric pressure (Pa)
C1 = 10;
C2 = 62.29;
C3 = 62.29;
C4 = 0.557;
C5 = 2391.3;
C6 = 0.609;

omega_i = 0.008; % Inlet humidity ratio (kg/kg dry air)
Tg_i = 303.15; % Inlet air temperature (K)

% Initial conditions
Ts_init = Tg_i;
omega_s_init = 0.008; % solid humidity ratio
omega_o_init = omega_i;
Tg_o_init = Tg_i;

% Sim parameters
N = 180; % Angular resolution
theta = linspace(0, 180, N); % Angular positions in degrees
dt = 1 / (15.5 * 60) / N; % Time step for rotation
tspan = [0, dt * N]; % Simulation time

% ODE system to access constants
function dYdt = desiccant_wheel_ode(t, Y)
% Unpack variables
omega_o = Y(1);
Ts = Y(2);
Tg_o = Y(3);
omega_s = Y(4);

% Saturation pressure and relative humidity
Ps = exp(5294 / Ts – 6); % Saturation pressure (Pa)
phi = (omega_s / 0.622) * P0 / ((0.622 + omega_s) * Ps);

% Empirical functions
s1 = (0.24 / 1.5) * (phi^0.333) * (0.622 * P0) / ((0.622 + omega_s)^2 * Ps);
s2 = (0.24 / 1.5) * (phi^0.333) * (5294 * phi) / (Ts^2);

% Differential equations
d_omega_dt = C1 * (omega_i – omega_o) + C2 * (omega_s – omega_o);
d_Ts_dt = C4 * C5 * (omega_o – omega_s) + C6 * (Tg_o – Ts);
d_Tg_dt = C1 * (Tg_i – Tg_o) + C3 * (Ts – Tg_o);
d_omega_s_dt = (s2 / s1) * d_Ts_dt + (C4 / s1) * (omega_o – omega_s);

% Output derivatives
dYdt = [d_omega_dt; d_Ts_dt; d_Tg_dt; d_omega_s_dt];
end

% ODE ode45
Y0 = [omega_o_init, Ts_init, Tg_o_init, omega_s_init]; % Initial conditions
[t, Y] = ode45(@desiccant_wheel_ode, tspan, Y0);

% Extract results
omega_o = Y(:, 1) * 1000; % Convert to g/kg dry air
Ts = Y(:, 2) – 273.15; % temp to °C

% Plot data
figure;
plot(theta, omega_o, ‘k’, ‘LineWidth’, 2);
xlabel(‘Degree’);
ylabel(‘Humidity [g/kg dry air]’);
title(‘Adsorption-side Outlet Humidity Ratio’);
grid on; P0 = 101325; % Atmospheric pressure (Pa)
C1 = 10;
C2 = 62.29;
C3 = 62.29;
C4 = 0.557;
C5 = 2391.3;
C6 = 0.609;

omega_i = 0.008; % Inlet humidity ratio (kg/kg dry air)
Tg_i = 303.15; % Inlet air temperature (K)

% Initial conditions
Ts_init = Tg_i;
omega_s_init = 0.008; % solid humidity ratio
omega_o_init = omega_i;
Tg_o_init = Tg_i;

% Sim parameters
N = 180; % Angular resolution
theta = linspace(0, 180, N); % Angular positions in degrees
dt = 1 / (15.5 * 60) / N; % Time step for rotation
tspan = [0, dt * N]; % Simulation time

% ODE system to access constants
function dYdt = desiccant_wheel_ode(t, Y)
% Unpack variables
omega_o = Y(1);
Ts = Y(2);
Tg_o = Y(3);
omega_s = Y(4);

% Saturation pressure and relative humidity
Ps = exp(5294 / Ts – 6); % Saturation pressure (Pa)
phi = (omega_s / 0.622) * P0 / ((0.622 + omega_s) * Ps);

% Empirical functions
s1 = (0.24 / 1.5) * (phi^0.333) * (0.622 * P0) / ((0.622 + omega_s)^2 * Ps);
s2 = (0.24 / 1.5) * (phi^0.333) * (5294 * phi) / (Ts^2);

% Differential equations
d_omega_dt = C1 * (omega_i – omega_o) + C2 * (omega_s – omega_o);
d_Ts_dt = C4 * C5 * (omega_o – omega_s) + C6 * (Tg_o – Ts);
d_Tg_dt = C1 * (Tg_i – Tg_o) + C3 * (Ts – Tg_o);
d_omega_s_dt = (s2 / s1) * d_Ts_dt + (C4 / s1) * (omega_o – omega_s);

% Output derivatives
dYdt = [d_omega_dt; d_Ts_dt; d_Tg_dt; d_omega_s_dt];
end

% ODE ode45
Y0 = [omega_o_init, Ts_init, Tg_o_init, omega_s_init]; % Initial conditions
[t, Y] = ode45(@desiccant_wheel_ode, tspan, Y0);

% Extract results
omega_o = Y(:, 1) * 1000; % Convert to g/kg dry air
Ts = Y(:, 2) – 273.15; % temp to °C

% Plot data
figure;
plot(theta, omega_o, ‘k’, ‘LineWidth’, 2);
xlabel(‘Degree’);
ylabel(‘Humidity [g/kg dry air]’);
title(‘Adsorption-side Outlet Humidity Ratio’);
grid on; desiccant wheel MATLAB Answers — New Questions

​

Dear Matlab users,

I am currently having a problem to define the nonlinear constraints.
The point of problem is with nonlinear constraints, it doesnt’ keep the number of initial population to fitnessfunction and left even only one individual.
But without nonlinear constraints option it keeps 45 individuals.
Please take a look this codes and may very greatful, when you write me any small tips.

size(InitPop) = 45
Vectorized= ‘on’

1. Starting InitialPopulation [45,9]
options = optimoptions(@ga, ‘PopulationSize’, 45, ‘InitialPopulation’, InitPop, …
‘PopulationType’, ‘doubleVector’, ”MaxStallGenerations’, 10, …
‘FunctionTolerance’, 1e-03, ‘MaxGenerations’, 100, ‘Vectorized’, ‘on’);

2.

[Value, fval, ~, output] = ga(@fitnessfunction, nvars, [], [], [], [], lb,ub, @nonlincon, [], options);

3. Definition of nonlinear constraints
function [c, ceq] = nonlincon(x)
c = 3 – (x(:, 1) + 3 * tan(x(:, 2)));% 3 – (x1 + 3tan(x2)) <= 0
ceq = [];
end

Thank you so much for time

Best regards.Dear Matlab users,

I am currently having a problem to define the nonlinear constraints.
The point of problem is with nonlinear constraints, it doesnt’ keep the number of initial population to fitnessfunction and left even only one individual.
But without nonlinear constraints option it keeps 45 individuals.
Please take a look this codes and may very greatful, when you write me any small tips.

size(InitPop) = 45
Vectorized= ‘on’

1. Starting InitialPopulation [45,9]
options = optimoptions(@ga, ‘PopulationSize’, 45, ‘InitialPopulation’, InitPop, …
‘PopulationType’, ‘doubleVector’, ”MaxStallGenerations’, 10, …
‘FunctionTolerance’, 1e-03, ‘MaxGenerations’, 100, ‘Vectorized’, ‘on’);

2.

[Value, fval, ~, output] = ga(@fitnessfunction, nvars, [], [], [], [], lb,ub, @nonlincon, [], options);

3. Definition of nonlinear constraints
function [c, ceq] = nonlincon(x)
c = 3 – (x(:, 1) + 3 * tan(x(:, 2)));% 3 – (x1 + 3tan(x2)) <= 0
ceq = [];
end

Thank you so much for time

Best regards. Dear Matlab users,

I am currently having a problem to define the nonlinear constraints.
The point of problem is with nonlinear constraints, it doesnt’ keep the number of initial population to fitnessfunction and left even only one individual.
But without nonlinear constraints option it keeps 45 individuals.
Please take a look this codes and may very greatful, when you write me any small tips.

size(InitPop) = 45
Vectorized= ‘on’

1. Starting InitialPopulation [45,9]
options = optimoptions(@ga, ‘PopulationSize’, 45, ‘InitialPopulation’, InitPop, …
‘PopulationType’, ‘doubleVector’, ”MaxStallGenerations’, 10, …
‘FunctionTolerance’, 1e-03, ‘MaxGenerations’, 100, ‘Vectorized’, ‘on’);

2.

[Value, fval, ~, output] = ga(@fitnessfunction, nvars, [], [], [], [], lb,ub, @nonlincon, [], options);

3. Definition of nonlinear constraints
function [c, ceq] = nonlincon(x)
c = 3 – (x(:, 1) + 3 * tan(x(:, 2)));% 3 – (x1 + 3tan(x2)) <= 0
ceq = [];
end

Thank you so much for time

Best regards. ga, genetic algorithm, optimization MATLAB Answers — New Questions

​

I’m trying to run the attached code for a Cassegrain antenna, changing the frequency from 18 GHz to the GPS frequencies 1575.42 MHz and 1227.6 MHz. I changed the frequencies and the following errors came up: "Error in em.MeshGeometry/updateMesh, Error in em.MeshGeometry/getMesh, and Error in em.EmStructures/analyze". Can you help resolve the errors?
Rp=0.3175;
fp=0.2536;
Rsub=0.033;
fhyp=0.1416;
ant=cassegrain(‘Radius’,[Rp Rsub],’FocalLength’,[fp fhyp]);
show(ant);

Exciter=design(hornConical,17.7e9);
Exciter.FeedWidth=3.4e-3;
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
show(Exciter);

ant=reflectorParabolic(‘Radius’,0.3175);
ant.Exciter=design(hornConical,17.7e9);
ant.Exciter.Tilt=90;
figure;
pattern(ant,18e9);

ant=cassegrain;
s=sparameters(ant,linspace(18e9,18.8e9,25));
figure;rfplot(s);

figure;
impedance(ant,linspace(18e9,18.8e9,25));

current(ant,18e9,’scale’,’log10′);

Exciter=design(horn,16.2e9);
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
ant.Exciter=Exciter;
show(ant);

figure;
pattern(ant,18e9);I’m trying to run the attached code for a Cassegrain antenna, changing the frequency from 18 GHz to the GPS frequencies 1575.42 MHz and 1227.6 MHz. I changed the frequencies and the following errors came up: "Error in em.MeshGeometry/updateMesh, Error in em.MeshGeometry/getMesh, and Error in em.EmStructures/analyze". Can you help resolve the errors?
Rp=0.3175;
fp=0.2536;
Rsub=0.033;
fhyp=0.1416;
ant=cassegrain(‘Radius’,[Rp Rsub],’FocalLength’,[fp fhyp]);
show(ant);

Exciter=design(hornConical,17.7e9);
Exciter.FeedWidth=3.4e-3;
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
show(Exciter);

ant=reflectorParabolic(‘Radius’,0.3175);
ant.Exciter=design(hornConical,17.7e9);
ant.Exciter.Tilt=90;
figure;
pattern(ant,18e9);

ant=cassegrain;
s=sparameters(ant,linspace(18e9,18.8e9,25));
figure;rfplot(s);

figure;
impedance(ant,linspace(18e9,18.8e9,25));

current(ant,18e9,’scale’,’log10′);

Exciter=design(horn,16.2e9);
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
ant.Exciter=Exciter;
show(ant);

figure;
pattern(ant,18e9); I’m trying to run the attached code for a Cassegrain antenna, changing the frequency from 18 GHz to the GPS frequencies 1575.42 MHz and 1227.6 MHz. I changed the frequencies and the following errors came up: "Error in em.MeshGeometry/updateMesh, Error in em.MeshGeometry/getMesh, and Error in em.EmStructures/analyze". Can you help resolve the errors?
Rp=0.3175;
fp=0.2536;
Rsub=0.033;
fhyp=0.1416;
ant=cassegrain(‘Radius’,[Rp Rsub],’FocalLength’,[fp fhyp]);
show(ant);

Exciter=design(hornConical,17.7e9);
Exciter.FeedWidth=3.4e-3;
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
show(Exciter);

ant=reflectorParabolic(‘Radius’,0.3175);
ant.Exciter=design(hornConical,17.7e9);
ant.Exciter.Tilt=90;
figure;
pattern(ant,18e9);

ant=cassegrain;
s=sparameters(ant,linspace(18e9,18.8e9,25));
figure;rfplot(s);

figure;
impedance(ant,linspace(18e9,18.8e9,25));

current(ant,18e9,’scale’,’log10′);

Exciter=design(horn,16.2e9);
Exciter.Tilt=270;
Exciter.TiltAxis=[0 1 0];
ant.Exciter=Exciter;
show(ant);

figure;
pattern(ant,18e9); cassegrain errors MATLAB Answers — New Questions

​