Hello Everyone!

I am having a strange problem with the app designer right now:

Every time I start the app designer and load the .mlapp-File of the Project I’m currently working on, there are multiple instances of the same EditField created. It began with only one copy of an EditField I manually added on the GUI. I didn’t think much about it and simply deleted the second instance. But it seems like this problem has grown exponentially with every start of app designer. Now, when I open up my Project 500 (!!) EditFields with the same name are created! ‘Simply’ deleting all of them is obviously not anymore an acceptable workaround – especially since ‘ctrl+a’ doesn’t work in the Component Browser and you have to selected them all manually…

I have no idea why the app designer creates these multiple instances on the startup – after I deleted all multiple instances, they also didn’t appear in the non-editable Function ‘createComponents(app)’ anymore.

Does anyone of you have an idea on that? Where in my code could be a bug causing this phenomenon? Or could this be an appdesigner-internal bug? The project itself is already quite big with ~260 components even without those multiple and unwanted EditFields…

I’m desperately looking for your help! Thanks a lot in advance.Hello Everyone!

I am having a strange problem with the app designer right now:

Every time I start the app designer and load the .mlapp-File of the Project I’m currently working on, there are multiple instances of the same EditField created. It began with only one copy of an EditField I manually added on the GUI. I didn’t think much about it and simply deleted the second instance. But it seems like this problem has grown exponentially with every start of app designer. Now, when I open up my Project 500 (!!) EditFields with the same name are created! ‘Simply’ deleting all of them is obviously not anymore an acceptable workaround – especially since ‘ctrl+a’ doesn’t work in the Component Browser and you have to selected them all manually…

I have no idea why the app designer creates these multiple instances on the startup – after I deleted all multiple instances, they also didn’t appear in the non-editable Function ‘createComponents(app)’ anymore.

Does anyone of you have an idea on that? Where in my code could be a bug causing this phenomenon? Or could this be an appdesigner-internal bug? The project itself is already quite big with ~260 components even without those multiple and unwanted EditFields…

I’m desperately looking for your help! Thanks a lot in advance. Hello Everyone!

I am having a strange problem with the app designer right now:

Every time I start the app designer and load the .mlapp-File of the Project I’m currently working on, there are multiple instances of the same EditField created. It began with only one copy of an EditField I manually added on the GUI. I didn’t think much about it and simply deleted the second instance. But it seems like this problem has grown exponentially with every start of app designer. Now, when I open up my Project 500 (!!) EditFields with the same name are created! ‘Simply’ deleting all of them is obviously not anymore an acceptable workaround – especially since ‘ctrl+a’ doesn’t work in the Component Browser and you have to selected them all manually…

I have no idea why the app designer creates these multiple instances on the startup – after I deleted all multiple instances, they also didn’t appear in the non-editable Function ‘createComponents(app)’ anymore.

Does anyone of you have an idea on that? Where in my code could be a bug causing this phenomenon? Or could this be an appdesigner-internal bug? The project itself is already quite big with ~260 components even without those multiple and unwanted EditFields…

I’m desperately looking for your help! Thanks a lot in advance. appdesigner, gui MATLAB Answers — New Questions

​

Good morning to everyone,
I have just bought a new laptop, downloaded 2025a, opened a figure I created with 2024b and BOOM I realized (with great regret!) that the "Figure Export Setup" have been removed…
In my academic and professional experience it was an essential feature to produce figures and plots in line with the specific request (journal article, technical report, presentation etc..) and coherent throughout the work for whant concerns font size, figure dimensions and features.
I know that now there’s the "uiexportdlg" https://www.mathworks.com/help/releases/R2025a/matlab/ref/uiexportdlg.html but, in my humble opinion, there is no comparison at all. Is there a way to be able to save profiles for the export of figures as was done in the old “File” > “Export Setup” mode? And, above all, does anyone know the reason why this feature was removed? I simply can’t get the point about it.

I’m not the reactionary kind of software users, but I’m seriously thinking about switching back to 2024b.Good morning to everyone,
I have just bought a new laptop, downloaded 2025a, opened a figure I created with 2024b and BOOM I realized (with great regret!) that the "Figure Export Setup" have been removed…
In my academic and professional experience it was an essential feature to produce figures and plots in line with the specific request (journal article, technical report, presentation etc..) and coherent throughout the work for whant concerns font size, figure dimensions and features.
I know that now there’s the "uiexportdlg" https://www.mathworks.com/help/releases/R2025a/matlab/ref/uiexportdlg.html but, in my humble opinion, there is no comparison at all. Is there a way to be able to save profiles for the export of figures as was done in the old “File” > “Export Setup” mode? And, above all, does anyone know the reason why this feature was removed? I simply can’t get the point about it.

I’m not the reactionary kind of software users, but I’m seriously thinking about switching back to 2024b. Good morning to everyone,
I have just bought a new laptop, downloaded 2025a, opened a figure I created with 2024b and BOOM I realized (with great regret!) that the "Figure Export Setup" have been removed…
In my academic and professional experience it was an essential feature to produce figures and plots in line with the specific request (journal article, technical report, presentation etc..) and coherent throughout the work for whant concerns font size, figure dimensions and features.
I know that now there’s the "uiexportdlg" https://www.mathworks.com/help/releases/R2025a/matlab/ref/uiexportdlg.html but, in my humble opinion, there is no comparison at all. Is there a way to be able to save profiles for the export of figures as was done in the old “File” > “Export Setup” mode? And, above all, does anyone know the reason why this feature was removed? I simply can’t get the point about it.

I’m not the reactionary kind of software users, but I’m seriously thinking about switching back to 2024b. figure export setup MATLAB Answers — New Questions

​

The examples above in MATLAB all produce BESS sets in discrete mode. I want to build a BESS set in phasor mode but don’t know how to do it because the battery block only allows running in discrete mode. I am very confused. Please help me.The examples above in MATLAB all produce BESS sets in discrete mode. I want to build a BESS set in phasor mode but don’t know how to do it because the battery block only allows running in discrete mode. I am very confused. Please help me. The examples above in MATLAB all produce BESS sets in discrete mode. I want to build a BESS set in phasor mode but don’t know how to do it because the battery block only allows running in discrete mode. I am very confused. Please help me. battery_system_management, bess phasor, battery phasor MATLAB Answers — New Questions

​

Hello, in the process of creating a problem for my students in a numerical methods class I found out that the backslash operation produces different results if you solve A^{T}x=b as "x=(A’)b" when compared to first defining "AT=A’" and then solving as "x=ATb". An example is given below. Order of operations suggests that in the first case, A should be transposed and then the backslash operation performed, which is forced in the second case. One would expect the results to then be the same. Why is there a difference at all?
n = 15;
m = 2*n+1;
t = fliplr(cos(pi/(n)*(0:n))).’;

A = [t.^(0:m);
t.^([0,0:m-1])*diag(max(0,0:m))];
b = (1./((0:m)’+1).*((1).^((0:m)’+1)-(-1).^(((0:m)’+1))));

x = (A’)b;

AT = A’;
x2 = ATb;

difference = norm(x-x2)Hello, in the process of creating a problem for my students in a numerical methods class I found out that the backslash operation produces different results if you solve A^{T}x=b as "x=(A’)b" when compared to first defining "AT=A’" and then solving as "x=ATb". An example is given below. Order of operations suggests that in the first case, A should be transposed and then the backslash operation performed, which is forced in the second case. One would expect the results to then be the same. Why is there a difference at all?
n = 15;
m = 2*n+1;
t = fliplr(cos(pi/(n)*(0:n))).’;

A = [t.^(0:m);
t.^([0,0:m-1])*diag(max(0,0:m))];
b = (1./((0:m)’+1).*((1).^((0:m)’+1)-(-1).^(((0:m)’+1))));

x = (A’)b;

AT = A’;
x2 = ATb;

difference = norm(x-x2) Hello, in the process of creating a problem for my students in a numerical methods class I found out that the backslash operation produces different results if you solve A^{T}x=b as "x=(A’)b" when compared to first defining "AT=A’" and then solving as "x=ATb". An example is given below. Order of operations suggests that in the first case, A should be transposed and then the backslash operation performed, which is forced in the second case. One would expect the results to then be the same. Why is there a difference at all?
n = 15;
m = 2*n+1;
t = fliplr(cos(pi/(n)*(0:n))).’;

A = [t.^(0:m);
t.^([0,0:m-1])*diag(max(0,0:m))];
b = (1./((0:m)’+1).*((1).^((0:m)’+1)-(-1).^(((0:m)’+1))));

x = (A’)b;

AT = A’;
x2 = ATb;

difference = norm(x-x2) backslash, transpose MATLAB Answers — New Questions

​

After set all the data into int16, by using HDL coder advisor i got this error message:

"Error:Error: HDL code generation from Stateflow failed: Entry event for state ‘Target_Search’ implied in ‘aero_guidance_model/Guidance/subsystem/Guidance Processor’. This is not allowed for an HDL target.

HDL code generation failed for ‘Guidance/subsystem/Guidance Processor’ (#59). Consult the Diagnostic Viewer for additional information."

Is it possible to convert the stateflow inside guidance subsystem into HDL for running hardware in the loop? Needing advice to call a stop or continue if its possible.

Example can be obtained https://www.mathworks.com/help/simulink/slref/designing-a-guidance-system-in-matlab-and-simulink.htmlAfter set all the data into int16, by using HDL coder advisor i got this error message:

"Error:Error: HDL code generation from Stateflow failed: Entry event for state ‘Target_Search’ implied in ‘aero_guidance_model/Guidance/subsystem/Guidance Processor’. This is not allowed for an HDL target.

HDL code generation failed for ‘Guidance/subsystem/Guidance Processor’ (#59). Consult the Diagnostic Viewer for additional information."

Is it possible to convert the stateflow inside guidance subsystem into HDL for running hardware in the loop? Needing advice to call a stop or continue if its possible.

Example can be obtained https://www.mathworks.com/help/simulink/slref/designing-a-guidance-system-in-matlab-and-simulink.html After set all the data into int16, by using HDL coder advisor i got this error message:

"Error:Error: HDL code generation from Stateflow failed: Entry event for state ‘Target_Search’ implied in ‘aero_guidance_model/Guidance/subsystem/Guidance Processor’. This is not allowed for an HDL target.

HDL code generation failed for ‘Guidance/subsystem/Guidance Processor’ (#59). Consult the Diagnostic Viewer for additional information."

Is it possible to convert the stateflow inside guidance subsystem into HDL for running hardware in the loop? Needing advice to call a stop or continue if its possible.

Example can be obtained https://www.mathworks.com/help/simulink/slref/designing-a-guidance-system-in-matlab-and-simulink.html hdl coder, hardware-in-the-loop MATLAB Answers — New Questions

​

Hi there
I was customising the standard MATLAB Bode plot to make the gridlines, major and minor bold,
change the x and y label size and bold them and also customise the axis.
Then I tried to plot the straight line approximation, but I can only seem to access the phase plot, and not the magnitude plot.
For the two axes, ax(1) and ax(2), I have called them mag_ax and phase_ax.
I can access the other properties fine, but when I tries to add the lines, my code doesn’t recognise the mag_ax.
Any advice greatly appreciated

close all

s = tf(‘s’);

G_1 = tf(1/(0.1*s + 1));

figure(‘WindowState’,’maximized’);

gg = bodeplot(G_1);

title(‘Bode plot of $frac{1}{0.1s+1}$’,…
‘Interpreter’,’latex’,’FontSize’,18,’FontWeight’,’bold’);

gg.Responses.Color = "blue";
gg.Responses.LineWidth = 2;
gg.XLabel.FontWeight = "bold";

% Get the handles of the axes
ax = findobj(gcf,’type’,’axes’);
mag_ax = ax(1);
phase_ax = ax(2);

mag_ax.XLim = [0.1 1000];
mag_ax.YLim= [-60 +60];

phase_ax.XLim = [0.1 1000];
phase_ax.YLim= [-90 +90];

phase_ax.YTick = (-90:45:90);

grid on
grid minor
phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.MinorGridLineStyle = ‘-‘;
mag_ax.MinorGridLineStyle = ‘-‘;

phase_ax.MinorGridLineWidth = 0.75;
mag_ax.MinorGridLineWidth = 0.75;

phase_ax.MinorGridAlpha = 0.5;
mag_ax.MinorGridAlpha = 0.5;

phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.GridAlpha = 1;
mag_ax.GridAlpha = 1;

mag_ax.YLabel.FontWeight = "bold";
mag_ax.YLabel.FontSize = 12;

phase_ax.YLabel.FontWeight = "bold";
phase_ax.YLabel.FontSize = 12;

for k = 1:length(phase_ax)
xticks(phase_ax(k), [0.1 1 10 100 1000]); % Set desired tick positions
end
for k = 1:length(phase_ax)
phase_ax(k).XTickLabel = arrayfun(@(x) sprintf(‘%.1f’, x), phase_ax(k).XTick, ‘UniformOutput’, false);
end

set(phase_ax,’fontweight’,’bold’)
set(mag_ax,’fontweight’,’bold’)

gg.Responses.Visible = ‘on’;
gg.Responses.Color = ‘blue’;

mag_ax.NextPlot = ‘add’;
line([0.1 10.0],[0 0],’linewidth’,2,’color’,’red’);
line([10.0 1000],[0 -40],’linewidth’,2,’color’,’red’);

phase_ax.NextPlot = ‘add’;
line([0.1 1.0],[0 0],’linewidth’,2,’color’,’red’);
line([1.0 100],[0 -90],’linewidth’,2,’color’,’red’);
line([100 1000],[-90 -90],’linewidth’,2,’color’,’red’);Hi there
I was customising the standard MATLAB Bode plot to make the gridlines, major and minor bold,
change the x and y label size and bold them and also customise the axis.
Then I tried to plot the straight line approximation, but I can only seem to access the phase plot, and not the magnitude plot.
For the two axes, ax(1) and ax(2), I have called them mag_ax and phase_ax.
I can access the other properties fine, but when I tries to add the lines, my code doesn’t recognise the mag_ax.
Any advice greatly appreciated

close all

s = tf(‘s’);

G_1 = tf(1/(0.1*s + 1));

figure(‘WindowState’,’maximized’);

gg = bodeplot(G_1);

title(‘Bode plot of $frac{1}{0.1s+1}$’,…
‘Interpreter’,’latex’,’FontSize’,18,’FontWeight’,’bold’);

gg.Responses.Color = "blue";
gg.Responses.LineWidth = 2;
gg.XLabel.FontWeight = "bold";

% Get the handles of the axes
ax = findobj(gcf,’type’,’axes’);
mag_ax = ax(1);
phase_ax = ax(2);

mag_ax.XLim = [0.1 1000];
mag_ax.YLim= [-60 +60];

phase_ax.XLim = [0.1 1000];
phase_ax.YLim= [-90 +90];

phase_ax.YTick = (-90:45:90);

grid on
grid minor
phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.MinorGridLineStyle = ‘-‘;
mag_ax.MinorGridLineStyle = ‘-‘;

phase_ax.MinorGridLineWidth = 0.75;
mag_ax.MinorGridLineWidth = 0.75;

phase_ax.MinorGridAlpha = 0.5;
mag_ax.MinorGridAlpha = 0.5;

phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.GridAlpha = 1;
mag_ax.GridAlpha = 1;

mag_ax.YLabel.FontWeight = "bold";
mag_ax.YLabel.FontSize = 12;

phase_ax.YLabel.FontWeight = "bold";
phase_ax.YLabel.FontSize = 12;

for k = 1:length(phase_ax)
xticks(phase_ax(k), [0.1 1 10 100 1000]); % Set desired tick positions
end
for k = 1:length(phase_ax)
phase_ax(k).XTickLabel = arrayfun(@(x) sprintf(‘%.1f’, x), phase_ax(k).XTick, ‘UniformOutput’, false);
end

set(phase_ax,’fontweight’,’bold’)
set(mag_ax,’fontweight’,’bold’)

gg.Responses.Visible = ‘on’;
gg.Responses.Color = ‘blue’;

mag_ax.NextPlot = ‘add’;
line([0.1 10.0],[0 0],’linewidth’,2,’color’,’red’);
line([10.0 1000],[0 -40],’linewidth’,2,’color’,’red’);

phase_ax.NextPlot = ‘add’;
line([0.1 1.0],[0 0],’linewidth’,2,’color’,’red’);
line([1.0 100],[0 -90],’linewidth’,2,’color’,’red’);
line([100 1000],[-90 -90],’linewidth’,2,’color’,’red’); Hi there
I was customising the standard MATLAB Bode plot to make the gridlines, major and minor bold,
change the x and y label size and bold them and also customise the axis.
Then I tried to plot the straight line approximation, but I can only seem to access the phase plot, and not the magnitude plot.
For the two axes, ax(1) and ax(2), I have called them mag_ax and phase_ax.
I can access the other properties fine, but when I tries to add the lines, my code doesn’t recognise the mag_ax.
Any advice greatly appreciated

close all

s = tf(‘s’);

G_1 = tf(1/(0.1*s + 1));

figure(‘WindowState’,’maximized’);

gg = bodeplot(G_1);

title(‘Bode plot of $frac{1}{0.1s+1}$’,…
‘Interpreter’,’latex’,’FontSize’,18,’FontWeight’,’bold’);

gg.Responses.Color = "blue";
gg.Responses.LineWidth = 2;
gg.XLabel.FontWeight = "bold";

% Get the handles of the axes
ax = findobj(gcf,’type’,’axes’);
mag_ax = ax(1);
phase_ax = ax(2);

mag_ax.XLim = [0.1 1000];
mag_ax.YLim= [-60 +60];

phase_ax.XLim = [0.1 1000];
phase_ax.YLim= [-90 +90];

phase_ax.YTick = (-90:45:90);

grid on
grid minor
phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.MinorGridLineStyle = ‘-‘;
mag_ax.MinorGridLineStyle = ‘-‘;

phase_ax.MinorGridLineWidth = 0.75;
mag_ax.MinorGridLineWidth = 0.75;

phase_ax.MinorGridAlpha = 0.5;
mag_ax.MinorGridAlpha = 0.5;

phase_ax.GridLineWidth = 1;
mag_ax.GridLineWidth = 1;

phase_ax.GridAlpha = 1;
mag_ax.GridAlpha = 1;

mag_ax.YLabel.FontWeight = "bold";
mag_ax.YLabel.FontSize = 12;

phase_ax.YLabel.FontWeight = "bold";
phase_ax.YLabel.FontSize = 12;

for k = 1:length(phase_ax)
xticks(phase_ax(k), [0.1 1 10 100 1000]); % Set desired tick positions
end
for k = 1:length(phase_ax)
phase_ax(k).XTickLabel = arrayfun(@(x) sprintf(‘%.1f’, x), phase_ax(k).XTick, ‘UniformOutput’, false);
end

set(phase_ax,’fontweight’,’bold’)
set(mag_ax,’fontweight’,’bold’)

gg.Responses.Visible = ‘on’;
gg.Responses.Color = ‘blue’;

mag_ax.NextPlot = ‘add’;
line([0.1 10.0],[0 0],’linewidth’,2,’color’,’red’);
line([10.0 1000],[0 -40],’linewidth’,2,’color’,’red’);

phase_ax.NextPlot = ‘add’;
line([0.1 1.0],[0 0],’linewidth’,2,’color’,’red’);
line([1.0 100],[0 -90],’linewidth’,2,’color’,’red’);
line([100 1000],[-90 -90],’linewidth’,2,’color’,’red’); bodeplot, matlab MATLAB Answers — New Questions

​

Dear all
I need help to add an algoritm or modified the exsit one to get the figures annexed in file:the following is my code:
clc; clear; close all;

%% —————- 1. Load Data ————————————
filename = ‘ThreeFaultModel_Modified.xlsx’;
data = xlsread(filename);

x = data(:,1);
g_Bouguer = data(:,2);
x = x(:); g_Bouguer = g_Bouguer(:);

%% —————- 2. User Inputs ———————————
N_layers = input(‘Enter number of layers: ‘);

rho_profiles = zeros(1, N_layers);
for i = 1:N_layers
rho_profiles(i) = input([‘Density of layer ‘ num2str(i) ‘: ‘]);
end

rho_basement = input(‘Basement density: ‘);
rho_contrasts = rho_profiles – rho_basement;

%% —————- 3. Inversion ————————
G = 6.67e-3;

z_inv = bsxfun(@rdivide, g_Bouguer, (2*pi*G*rho_contrasts));

%% —————- 4. Detect ONE Main Fault ————————
VG = gradient(z_inv(:,1), x);

% اختيار أقوى تغير (fault)
[~, idx_fault] = max(abs(VG));

MainFault.x = x(idx_fault);
MainFault.dip = atan2d(VG(idx_fault),1); % تقدير تقريبي
MainFault.throw = max(z_inv(:,1)) * 0.3; % تقدير بسيط
MainFault.type = ‘Normal’;

if MainFault.dip > 0
MainFault.type = ‘Reverse’;
end
if abs(MainFault.dip) > 85
MainFault.type = ‘Vertical’;
end

%% —————- 5. Apply ONE Fault to All Layers —————-
z_layers = z_inv;

for L = 1:N_layers
if strcmp(MainFault.type,’Reverse’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) + MainFault.throw;

elseif strcmp(MainFault.type,’Normal’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) – MainFault.throw;

else
% Vertical → no vertical shift
end
end

%% —————- 6. Plot Depth ————————
figure; hold on;
for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end
set(gca,’YDir’,’reverse’);
grid on;
title(‘Depth with Single Fault’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
legend(arrayfun(@(k) sprintf(‘Layer %d’,k),1:N_layers,’UniformOutput’,false));

%% —————- 7. Plot Single Fault Plane ————————
z_max = max(z_layers(:));

if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’r–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’r–‘,’LineWidth’,2);
end

%% —————- 8. Dip Line from Surface ————————
figure; hold on;

for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
grid on;
title(‘Single Fault Dip Line’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);

z1 = max(z_layers(:,1));
dip = MainFault.dip;
x0 = MainFault.x;

if abs(dip) > 89
x1 = x0;
else
x1 = x0 + z1 / tand(dip);
end

plot([x0 x1],[0 z1],’r–‘,’LineWidth’,2);
text((x0+x1)/2,0.3,’Main Fault’,’Color’,’r’);

%% —————- 9. Geological Section ————————
figure; hold on;

colors = lines(N_layers);

for i = N_layers:-1:1
if i == 1
top = zeros(size(x));
else
top = z_layers(:,i-1);
end

bot = z_layers(:,i);

fill([x;flipud(x)], [top;flipud(bot)], colors(i,:), …
‘EdgeColor’,’k’,’FaceAlpha’,0.6);
end

% رسم الصدع
if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’k–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’k–‘,’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
title(‘Geological Cross Section with Single Fault’);

legend(arrayfun(@(i) sprintf(‘Layer %d’,i),1:N_layers,’UniformOutput’,false));
we use the following paramters:
Enter number of layers: 3
Density of layer 1: 2.20
Density of layer 2: 2.25
Density of layer 3: 2.30
Basement density: 2.67Dear all
I need help to add an algoritm or modified the exsit one to get the figures annexed in file:the following is my code:
clc; clear; close all;

%% —————- 1. Load Data ————————————
filename = ‘ThreeFaultModel_Modified.xlsx’;
data = xlsread(filename);

x = data(:,1);
g_Bouguer = data(:,2);
x = x(:); g_Bouguer = g_Bouguer(:);

%% —————- 2. User Inputs ———————————
N_layers = input(‘Enter number of layers: ‘);

rho_profiles = zeros(1, N_layers);
for i = 1:N_layers
rho_profiles(i) = input([‘Density of layer ‘ num2str(i) ‘: ‘]);
end

rho_basement = input(‘Basement density: ‘);
rho_contrasts = rho_profiles – rho_basement;

%% —————- 3. Inversion ————————
G = 6.67e-3;

z_inv = bsxfun(@rdivide, g_Bouguer, (2*pi*G*rho_contrasts));

%% —————- 4. Detect ONE Main Fault ————————
VG = gradient(z_inv(:,1), x);

% اختيار أقوى تغير (fault)
[~, idx_fault] = max(abs(VG));

MainFault.x = x(idx_fault);
MainFault.dip = atan2d(VG(idx_fault),1); % تقدير تقريبي
MainFault.throw = max(z_inv(:,1)) * 0.3; % تقدير بسيط
MainFault.type = ‘Normal’;

if MainFault.dip > 0
MainFault.type = ‘Reverse’;
end
if abs(MainFault.dip) > 85
MainFault.type = ‘Vertical’;
end

%% —————- 5. Apply ONE Fault to All Layers —————-
z_layers = z_inv;

for L = 1:N_layers
if strcmp(MainFault.type,’Reverse’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) + MainFault.throw;

elseif strcmp(MainFault.type,’Normal’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) – MainFault.throw;

else
% Vertical → no vertical shift
end
end

%% —————- 6. Plot Depth ————————
figure; hold on;
for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end
set(gca,’YDir’,’reverse’);
grid on;
title(‘Depth with Single Fault’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
legend(arrayfun(@(k) sprintf(‘Layer %d’,k),1:N_layers,’UniformOutput’,false));

%% —————- 7. Plot Single Fault Plane ————————
z_max = max(z_layers(:));

if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’r–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’r–‘,’LineWidth’,2);
end

%% —————- 8. Dip Line from Surface ————————
figure; hold on;

for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
grid on;
title(‘Single Fault Dip Line’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);

z1 = max(z_layers(:,1));
dip = MainFault.dip;
x0 = MainFault.x;

if abs(dip) > 89
x1 = x0;
else
x1 = x0 + z1 / tand(dip);
end

plot([x0 x1],[0 z1],’r–‘,’LineWidth’,2);
text((x0+x1)/2,0.3,’Main Fault’,’Color’,’r’);

%% —————- 9. Geological Section ————————
figure; hold on;

colors = lines(N_layers);

for i = N_layers:-1:1
if i == 1
top = zeros(size(x));
else
top = z_layers(:,i-1);
end

bot = z_layers(:,i);

fill([x;flipud(x)], [top;flipud(bot)], colors(i,:), …
‘EdgeColor’,’k’,’FaceAlpha’,0.6);
end

% رسم الصدع
if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’k–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’k–‘,’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
title(‘Geological Cross Section with Single Fault’);

legend(arrayfun(@(i) sprintf(‘Layer %d’,i),1:N_layers,’UniformOutput’,false));
we use the following paramters:
Enter number of layers: 3
Density of layer 1: 2.20
Density of layer 2: 2.25
Density of layer 3: 2.30
Basement density: 2.67 Dear all
I need help to add an algoritm or modified the exsit one to get the figures annexed in file:the following is my code:
clc; clear; close all;

%% —————- 1. Load Data ————————————
filename = ‘ThreeFaultModel_Modified.xlsx’;
data = xlsread(filename);

x = data(:,1);
g_Bouguer = data(:,2);
x = x(:); g_Bouguer = g_Bouguer(:);

%% —————- 2. User Inputs ———————————
N_layers = input(‘Enter number of layers: ‘);

rho_profiles = zeros(1, N_layers);
for i = 1:N_layers
rho_profiles(i) = input([‘Density of layer ‘ num2str(i) ‘: ‘]);
end

rho_basement = input(‘Basement density: ‘);
rho_contrasts = rho_profiles – rho_basement;

%% —————- 3. Inversion ————————
G = 6.67e-3;

z_inv = bsxfun(@rdivide, g_Bouguer, (2*pi*G*rho_contrasts));

%% —————- 4. Detect ONE Main Fault ————————
VG = gradient(z_inv(:,1), x);

% اختيار أقوى تغير (fault)
[~, idx_fault] = max(abs(VG));

MainFault.x = x(idx_fault);
MainFault.dip = atan2d(VG(idx_fault),1); % تقدير تقريبي
MainFault.throw = max(z_inv(:,1)) * 0.3; % تقدير بسيط
MainFault.type = ‘Normal’;

if MainFault.dip > 0
MainFault.type = ‘Reverse’;
end
if abs(MainFault.dip) > 85
MainFault.type = ‘Vertical’;
end

%% —————- 5. Apply ONE Fault to All Layers —————-
z_layers = z_inv;

for L = 1:N_layers
if strcmp(MainFault.type,’Reverse’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) + MainFault.throw;

elseif strcmp(MainFault.type,’Normal’)
z_layers(idx_fault:end, L:end) = z_layers(idx_fault:end, L:end) – MainFault.throw;

else
% Vertical → no vertical shift
end
end

%% —————- 6. Plot Depth ————————
figure; hold on;
for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end
set(gca,’YDir’,’reverse’);
grid on;
title(‘Depth with Single Fault’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
legend(arrayfun(@(k) sprintf(‘Layer %d’,k),1:N_layers,’UniformOutput’,false));

%% —————- 7. Plot Single Fault Plane ————————
z_max = max(z_layers(:));

if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’r–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’r–‘,’LineWidth’,2);
end

%% —————- 8. Dip Line from Surface ————————
figure; hold on;

for i = 1:N_layers
plot(x, z_layers(:,i),’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
grid on;
title(‘Single Fault Dip Line’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);

z1 = max(z_layers(:,1));
dip = MainFault.dip;
x0 = MainFault.x;

if abs(dip) > 89
x1 = x0;
else
x1 = x0 + z1 / tand(dip);
end

plot([x0 x1],[0 z1],’r–‘,’LineWidth’,2);
text((x0+x1)/2,0.3,’Main Fault’,’Color’,’r’);

%% —————- 9. Geological Section ————————
figure; hold on;

colors = lines(N_layers);

for i = N_layers:-1:1
if i == 1
top = zeros(size(x));
else
top = z_layers(:,i-1);
end

bot = z_layers(:,i);

fill([x;flipud(x)], [top;flipud(bot)], colors(i,:), …
‘EdgeColor’,’k’,’FaceAlpha’,0.6);
end

% رسم الصدع
if strcmp(MainFault.type,’Vertical’)
plot([MainFault.x MainFault.x],[0 z_max],’k–‘,’LineWidth’,2);
else
dx = z_max / tand(MainFault.dip);

if MainFault.dip > 0
x2 = MainFault.x – dx;
else
x2 = MainFault.x + dx;
end

plot([MainFault.x x2],[0 z_max],’k–‘,’LineWidth’,2);
end

set(gca,’YDir’,’reverse’);
xlabel(‘Distance (km)’);
ylabel(‘Depth (km)’);
title(‘Geological Cross Section with Single Fault’);

legend(arrayfun(@(i) sprintf(‘Layer %d’,i),1:N_layers,’UniformOutput’,false));
we use the following paramters:
Enter number of layers: 3
Density of layer 1: 2.20
Density of layer 2: 2.25
Density of layer 3: 2.30
Basement density: 2.67 matlab, drawing, plane MATLAB Answers — New Questions

​

I want to draw relation between T ( y-axis) and x ( x-axis) with change parameter B1.

Constant:
B2=0.01;B4=0.1;
vector

Parameter
B1=0.1:0.1:0.5I want to draw relation between T ( y-axis) and x ( x-axis) with change parameter B1.

Constant:
B2=0.01;B4=0.1;
vector

Parameter
B1=0.1:0.1:0.5 I want to draw relation between T ( y-axis) and x ( x-axis) with change parameter B1.

Constant:
B2=0.01;B4=0.1;
vector

Parameter
B1=0.1:0.1:0.5 plot MATLAB Answers — New Questions

​

I am working on project on UAV. MATLAB code is working fine, but generated output does not show any UAVs. I can share some code here if someone is interested.I am working on project on UAV. MATLAB code is working fine, but generated output does not show any UAVs. I can share some code here if someone is interested. I am working on project on UAV. MATLAB code is working fine, but generated output does not show any UAVs. I can share some code here if someone is interested. aerial vehicle MATLAB Answers — New Questions

​

Suppose there are two square symmetric matrices [K] and [M]. Since my work is related with structural analysis so they are typically stiffness and mass matrices. I am doing wave propagation study for a narrow band input load. In frequency domain, load amplitude is significant only between frequency [omega1 omega2]. So I am interested to find all eigen value and corresponding eigen vector only in this range.

The eigen value problem statement is
[K]=lambda[M]Suppose there are two square symmetric matrices [K] and [M]. Since my work is related with structural analysis so they are typically stiffness and mass matrices. I am doing wave propagation study for a narrow band input load. In frequency domain, load amplitude is significant only between frequency [omega1 omega2]. So I am interested to find all eigen value and corresponding eigen vector only in this range.

The eigen value problem statement is
[K]=lambda[M] Suppose there are two square symmetric matrices [K] and [M]. Since my work is related with structural analysis so they are typically stiffness and mass matrices. I am doing wave propagation study for a narrow band input load. In frequency domain, load amplitude is significant only between frequency [omega1 omega2]. So I am interested to find all eigen value and corresponding eigen vector only in this range.

The eigen value problem statement is
[K]=lambda[M] all eigen value between given range MATLAB Answers — New Questions

​