Hi,
I need to find the average wave from several in a trace.
I would like to produce a figure where the action potentials are stacked, with their average incorporated into the figure.
I only managed to find peaks in the waveform and plot them but I don’t know how to proceed.
%What I have so far:

clear
close all

D = load(‘C:UsersrsanaDesktopAPAnalysisS2C1_17082020_WM4peaks.mat’);

F = fieldnames(D);
C = struct2cell(D);
vm = C{1}.values; % Signal Vector
L1 = numel(vm); % Length

% Set Time Vector for ms Interval
t = linspace(0, L1, L1).’*C{1}.interval;

% Determine peaks, time interval, and cutoff for APs
[peaks, int, W] = findpeaks(vm, ‘WidthReference’, ‘halfprom’);

figure
plot(t, vm)
hold on
plot(t(int), peaks, ‘+r’)
hold off
xlabel(‘Time (s)’)
ylabel(‘Voltage (mV)’)

%%%
The above generates labelled peaks. So I am wondering how to extract, stack, and average them.

Thank you,
RoksanaHi,
I need to find the average wave from several in a trace.
I would like to produce a figure where the action potentials are stacked, with their average incorporated into the figure.
I only managed to find peaks in the waveform and plot them but I don’t know how to proceed.
%What I have so far:

clear
close all

D = load(‘C:UsersrsanaDesktopAPAnalysisS2C1_17082020_WM4peaks.mat’);

F = fieldnames(D);
C = struct2cell(D);
vm = C{1}.values; % Signal Vector
L1 = numel(vm); % Length

% Set Time Vector for ms Interval
t = linspace(0, L1, L1).’*C{1}.interval;

% Determine peaks, time interval, and cutoff for APs
[peaks, int, W] = findpeaks(vm, ‘WidthReference’, ‘halfprom’);

figure
plot(t, vm)
hold on
plot(t(int), peaks, ‘+r’)
hold off
xlabel(‘Time (s)’)
ylabel(‘Voltage (mV)’)

%%%
The above generates labelled peaks. So I am wondering how to extract, stack, and average them.

Thank you,
Roksana Hi,
I need to find the average wave from several in a trace.
I would like to produce a figure where the action potentials are stacked, with their average incorporated into the figure.
I only managed to find peaks in the waveform and plot them but I don’t know how to proceed.
%What I have so far:

clear
close all

D = load(‘C:UsersrsanaDesktopAPAnalysisS2C1_17082020_WM4peaks.mat’);

F = fieldnames(D);
C = struct2cell(D);
vm = C{1}.values; % Signal Vector
L1 = numel(vm); % Length

% Set Time Vector for ms Interval
t = linspace(0, L1, L1).’*C{1}.interval;

% Determine peaks, time interval, and cutoff for APs
[peaks, int, W] = findpeaks(vm, ‘WidthReference’, ‘halfprom’);

figure
plot(t, vm)
hold on
plot(t(int), peaks, ‘+r’)
hold off
xlabel(‘Time (s)’)
ylabel(‘Voltage (mV)’)

%%%
The above generates labelled peaks. So I am wondering how to extract, stack, and average them.

Thank you,
Roksana waveform analysis, average action potential MATLAB Answers — New Questions

​

I have plugged in math.sin and sin into my equation and I continue getting errors when trying to run the script. I am lost and I am not sure how to correct this.
(8+6)*abs(4-15) + math.sin(math.log10(18))-(4+6) * math.tan(8/18+1)
Unable to resolve the name ‘math.sin’.
(8+6)*abs(4-15) + sin*(math.log10(18))-(4+6) * math.tan(8/18+1)
Error using sin. Not enough input arguments.I have plugged in math.sin and sin into my equation and I continue getting errors when trying to run the script. I am lost and I am not sure how to correct this.
(8+6)*abs(4-15) + math.sin(math.log10(18))-(4+6) * math.tan(8/18+1)
Unable to resolve the name ‘math.sin’.
(8+6)*abs(4-15) + sin*(math.log10(18))-(4+6) * math.tan(8/18+1)
Error using sin. Not enough input arguments. I have plugged in math.sin and sin into my equation and I continue getting errors when trying to run the script. I am lost and I am not sure how to correct this.
(8+6)*abs(4-15) + math.sin(math.log10(18))-(4+6) * math.tan(8/18+1)
Unable to resolve the name ‘math.sin’.
(8+6)*abs(4-15) + sin*(math.log10(18))-(4+6) * math.tan(8/18+1)
Error using sin. Not enough input arguments. matlab MATLAB Answers — New Questions

​

Hello, in the photo below the start and end points are stable,the rate of change between them changes.
In the code below i tried to use exponential curves but different rate gives me different end point.
is there some wat the curves will start and end at the same points as shown below?
Thanks.
% Define the range of x values
x = linspace(0, 1, 10); % 0 to 1 with 100 points

% Define different rates for the exponential rise
rate1 = 10; % Faster rise
rate2 = 9.9; % Moderate rise
rate3 = 2; % Slower rise

% Compute the exponential values
y1 = exp(rate1 * x) – 1; % Exponential curve 1
y2 = exp(rate2 * x) – 1; % Exponential curve 2
y3 = exp(rate3 * x) – 1; % Exponential curve 3

% Plot the exponential curves
figure;
plot(x, y1, ‘b’, ‘LineWidth’, 1.5); hold on;
plot(x, y2, ‘r’, ‘LineWidth’, 1.5);
plot(x, y3, ‘g’, ‘LineWidth’, 1.5);

% Add labels and legend
xlabel(‘x-axis’);
ylabel(‘y-axis’);
legend({‘Rate = 10’, ‘Rate = 5’, ‘Rate = 2’}, ‘Location’, ‘NorthEast’);
title(‘Exponential Rise at Different Rates’);
grid on;Hello, in the photo below the start and end points are stable,the rate of change between them changes.
In the code below i tried to use exponential curves but different rate gives me different end point.
is there some wat the curves will start and end at the same points as shown below?
Thanks.
% Define the range of x values
x = linspace(0, 1, 10); % 0 to 1 with 100 points

% Define different rates for the exponential rise
rate1 = 10; % Faster rise
rate2 = 9.9; % Moderate rise
rate3 = 2; % Slower rise

% Compute the exponential values
y1 = exp(rate1 * x) – 1; % Exponential curve 1
y2 = exp(rate2 * x) – 1; % Exponential curve 2
y3 = exp(rate3 * x) – 1; % Exponential curve 3

% Plot the exponential curves
figure;
plot(x, y1, ‘b’, ‘LineWidth’, 1.5); hold on;
plot(x, y2, ‘r’, ‘LineWidth’, 1.5);
plot(x, y3, ‘g’, ‘LineWidth’, 1.5);

% Add labels and legend
xlabel(‘x-axis’);
ylabel(‘y-axis’);
legend({‘Rate = 10’, ‘Rate = 5’, ‘Rate = 2’}, ‘Location’, ‘NorthEast’);
title(‘Exponential Rise at Different Rates’);
grid on; Hello, in the photo below the start and end points are stable,the rate of change between them changes.
In the code below i tried to use exponential curves but different rate gives me different end point.
is there some wat the curves will start and end at the same points as shown below?
Thanks.
% Define the range of x values
x = linspace(0, 1, 10); % 0 to 1 with 100 points

% Define different rates for the exponential rise
rate1 = 10; % Faster rise
rate2 = 9.9; % Moderate rise
rate3 = 2; % Slower rise

% Compute the exponential values
y1 = exp(rate1 * x) – 1; % Exponential curve 1
y2 = exp(rate2 * x) – 1; % Exponential curve 2
y3 = exp(rate3 * x) – 1; % Exponential curve 3

% Plot the exponential curves
figure;
plot(x, y1, ‘b’, ‘LineWidth’, 1.5); hold on;
plot(x, y2, ‘r’, ‘LineWidth’, 1.5);
plot(x, y3, ‘g’, ‘LineWidth’, 1.5);

% Add labels and legend
xlabel(‘x-axis’);
ylabel(‘y-axis’);
legend({‘Rate = 10’, ‘Rate = 5’, ‘Rate = 2’}, ‘Location’, ‘NorthEast’);
title(‘Exponential Rise at Different Rates’);
grid on; rate of rise steady ends MATLAB Answers — New Questions

​

hi,
i have an excel file (say file1.xlsx)numeric data with 2000 rows and 25 colums. and i have 20 such excel sheets in this file with the same number of rows and colums.
i want to create a matlab file which will create graph of sheet1 of this excel file and save it with all the settings(labels, legends, font, etc), so that next time i only select the sheet2 and it will auomatically create graph and so on for the next 20 sheets of the same excel file.
is there any tool in matlab for such plotting or i have to write piece of code.
if coding is the solution then please help me through this.
thank you in advancehi,
i have an excel file (say file1.xlsx)numeric data with 2000 rows and 25 colums. and i have 20 such excel sheets in this file with the same number of rows and colums.
i want to create a matlab file which will create graph of sheet1 of this excel file and save it with all the settings(labels, legends, font, etc), so that next time i only select the sheet2 and it will auomatically create graph and so on for the next 20 sheets of the same excel file.
is there any tool in matlab for such plotting or i have to write piece of code.
if coding is the solution then please help me through this.
thank you in advance hi,
i have an excel file (say file1.xlsx)numeric data with 2000 rows and 25 colums. and i have 20 such excel sheets in this file with the same number of rows and colums.
i want to create a matlab file which will create graph of sheet1 of this excel file and save it with all the settings(labels, legends, font, etc), so that next time i only select the sheet2 and it will auomatically create graph and so on for the next 20 sheets of the same excel file.
is there any tool in matlab for such plotting or i have to write piece of code.
if coding is the solution then please help me through this.
thank you in advance xlsread, plotting MATLAB Answers — New Questions

​

Hi, I took some mesurements of ramp response for furnace to choose parameters for PID, but I do not know how to define model of it. The ambient temperature was about 15°C, ramp was 120°C/h, final temp was 60°C, but it was tested in closed loop with P=1. The data I received are temperature measurements depending on time. Do you have any ideas how to do it ?Hi, I took some mesurements of ramp response for furnace to choose parameters for PID, but I do not know how to define model of it. The ambient temperature was about 15°C, ramp was 120°C/h, final temp was 60°C, but it was tested in closed loop with P=1. The data I received are temperature measurements depending on time. Do you have any ideas how to do it ? Hi, I took some mesurements of ramp response for furnace to choose parameters for PID, but I do not know how to define model of it. The ambient temperature was about 15°C, ramp was 120°C/h, final temp was 60°C, but it was tested in closed loop with P=1. The data I received are temperature measurements depending on time. Do you have any ideas how to do it ? transfer function, model, simulink, matlab, state-space MATLAB Answers — New Questions

​

Write a function that receives an input argument of a number of kilometers (km) and will output the conversions into both miles (mi) and US nautical miles (n.m.). Use the conversions: 1 km = 0.621 mi and 1 km = 0.540 n.m.Write a function that receives an input argument of a number of kilometers (km) and will output the conversions into both miles (mi) and US nautical miles (n.m.). Use the conversions: 1 km = 0.621 mi and 1 km = 0.540 n.m. Write a function that receives an input argument of a number of kilometers (km) and will output the conversions into both miles (mi) and US nautical miles (n.m.). Use the conversions: 1 km = 0.621 mi and 1 km = 0.540 n.m. functio MATLAB Answers — New Questions

​

Hi
I have system (NOT transfer function) that gives vector numbers. Also, I have specific desire output. I would like to add PID controller as a control to reduce the error between the desire and output based on specific variable (d). Can you please help me with this situation
For example,
A=[0 3 12.5 12.8 11.8 12.3 11.9 12]; % practical data
d=1; % variable
OUT=A*d; % output
Vdesire=12; % desire outputHi
I have system (NOT transfer function) that gives vector numbers. Also, I have specific desire output. I would like to add PID controller as a control to reduce the error between the desire and output based on specific variable (d). Can you please help me with this situation
For example,
A=[0 3 12.5 12.8 11.8 12.3 11.9 12]; % practical data
d=1; % variable
OUT=A*d; % output
Vdesire=12; % desire output Hi
I have system (NOT transfer function) that gives vector numbers. Also, I have specific desire output. I would like to add PID controller as a control to reduce the error between the desire and output based on specific variable (d). Can you please help me with this situation
For example,
A=[0 3 12.5 12.8 11.8 12.3 11.9 12]; % practical data
d=1; % variable
OUT=A*d; % output
Vdesire=12; % desire output pid, controller MATLAB Answers — New Questions

​

Hello All
I am using the Industrial Robot Models in Simscape by Steve Miller, and need help editing a robot model. Specifically, I would like to modify link dimensions, add joints, or adjust dynamics. Any guidance or tips would be greatly appreciated!
Thank you in advance
ChintanHello All
I am using the Industrial Robot Models in Simscape by Steve Miller, and need help editing a robot model. Specifically, I would like to modify link dimensions, add joints, or adjust dynamics. Any guidance or tips would be greatly appreciated!
Thank you in advance
Chintan Hello All
I am using the Industrial Robot Models in Simscape by Steve Miller, and need help editing a robot model. Specifically, I would like to modify link dimensions, add joints, or adjust dynamics. Any guidance or tips would be greatly appreciated!
Thank you in advance
Chintan industrial robot models, robotic manipulator, 6-dof pick and place manipulator, simscape model MATLAB Answers — New Questions

​

I have a Simulink model that uses variables set in Matlab code. However it keeps giving me the same error.
I have this function
function y = tclabsim(t, x0, u, p)
% Ensure input sizes match
if length(t) ~= length(u)
error(‘Time vector t and input profile u must have the same length.’);
end

% Extract parameters from p
U = p(1); % Heat transfer coefficient
alpha = p(2); % Ambient heat loss parameter
tau = p(3); % Time constant

% Initial conditions
x0_T = x0; % Initial temperature

% Ensure t and u are column vectors
t = t(:);
u = u(:);

% Create input matrices
inputMatrix = [t, u];
U_matrix = [0, U];
alpha_matrix = [0, alpha];
tau_matrix = [0, tau];
x0_T_matrix = [0, x0_T];

% Set simulation stop time
stopTime = t(end);

% Assign all required variables to base workspace
assignin(‘base’, ‘x0_T’, x0_T_matrix);
assignin(‘base’, ‘U’, U_matrix);
assignin(‘base’, ‘alpha’, alpha_matrix);
assignin(‘base’, ‘tau’, tau_matrix);
assignin(‘base’, ‘inputMatrix’, inputMatrix);

% Run the Simulink model
simOut = sim(‘SingleHeaterTCLab_1.slx’, ‘SrcWorkspace’, ‘current’, …
‘StopTime’, num2str(stopTime));

% Extract the output from the simulation
y = simOut.Ts; % Sensor temperature output
end
And I call the function like this
% Time vector (0 to 1000 seconds with 1-second steps)
t = (0:999)’; % Make sure t is a column vector

% Initial temperature (20°C)
x0 = 20;

% Heater input (step response: 50% heater power after 100s)
u = zeros(size(t));
u(t >= 100) = 50;

% Parameters [U, alpha, tau]
p = [0.01, 0.01, 200]; % Example values for U, alpha, tau

% Simulate the system
y = tclabsim(t, x0, u, p);
However, every time I try to run it gives the same error.

I have tried to load them as simple double and I have also tried a matrix with t on the first column and the same value in every line of the second column. They also gave the same error.
The block in Simulink has these properties. And I have also tried to change the sample time to -1.

This is the entire simulink if it helps.

How can I fix it?
Thank youI have a Simulink model that uses variables set in Matlab code. However it keeps giving me the same error.
I have this function
function y = tclabsim(t, x0, u, p)
% Ensure input sizes match
if length(t) ~= length(u)
error(‘Time vector t and input profile u must have the same length.’);
end

% Extract parameters from p
U = p(1); % Heat transfer coefficient
alpha = p(2); % Ambient heat loss parameter
tau = p(3); % Time constant

% Initial conditions
x0_T = x0; % Initial temperature

% Ensure t and u are column vectors
t = t(:);
u = u(:);

% Create input matrices
inputMatrix = [t, u];
U_matrix = [0, U];
alpha_matrix = [0, alpha];
tau_matrix = [0, tau];
x0_T_matrix = [0, x0_T];

% Set simulation stop time
stopTime = t(end);

% Assign all required variables to base workspace
assignin(‘base’, ‘x0_T’, x0_T_matrix);
assignin(‘base’, ‘U’, U_matrix);
assignin(‘base’, ‘alpha’, alpha_matrix);
assignin(‘base’, ‘tau’, tau_matrix);
assignin(‘base’, ‘inputMatrix’, inputMatrix);

% Run the Simulink model
simOut = sim(‘SingleHeaterTCLab_1.slx’, ‘SrcWorkspace’, ‘current’, …
‘StopTime’, num2str(stopTime));

% Extract the output from the simulation
y = simOut.Ts; % Sensor temperature output
end
And I call the function like this
% Time vector (0 to 1000 seconds with 1-second steps)
t = (0:999)’; % Make sure t is a column vector

% Initial temperature (20°C)
x0 = 20;

% Heater input (step response: 50% heater power after 100s)
u = zeros(size(t));
u(t >= 100) = 50;

% Parameters [U, alpha, tau]
p = [0.01, 0.01, 200]; % Example values for U, alpha, tau

% Simulate the system
y = tclabsim(t, x0, u, p);
However, every time I try to run it gives the same error.

I have tried to load them as simple double and I have also tried a matrix with t on the first column and the same value in every line of the second column. They also gave the same error.
The block in Simulink has these properties. And I have also tried to change the sample time to -1.

This is the entire simulink if it helps.

How can I fix it?
Thank you I have a Simulink model that uses variables set in Matlab code. However it keeps giving me the same error.
I have this function
function y = tclabsim(t, x0, u, p)
% Ensure input sizes match
if length(t) ~= length(u)
error(‘Time vector t and input profile u must have the same length.’);
end

% Extract parameters from p
U = p(1); % Heat transfer coefficient
alpha = p(2); % Ambient heat loss parameter
tau = p(3); % Time constant

% Initial conditions
x0_T = x0; % Initial temperature

% Ensure t and u are column vectors
t = t(:);
u = u(:);

% Create input matrices
inputMatrix = [t, u];
U_matrix = [0, U];
alpha_matrix = [0, alpha];
tau_matrix = [0, tau];
x0_T_matrix = [0, x0_T];

% Set simulation stop time
stopTime = t(end);

% Assign all required variables to base workspace
assignin(‘base’, ‘x0_T’, x0_T_matrix);
assignin(‘base’, ‘U’, U_matrix);
assignin(‘base’, ‘alpha’, alpha_matrix);
assignin(‘base’, ‘tau’, tau_matrix);
assignin(‘base’, ‘inputMatrix’, inputMatrix);

% Run the Simulink model
simOut = sim(‘SingleHeaterTCLab_1.slx’, ‘SrcWorkspace’, ‘current’, …
‘StopTime’, num2str(stopTime));

% Extract the output from the simulation
y = simOut.Ts; % Sensor temperature output
end
And I call the function like this
% Time vector (0 to 1000 seconds with 1-second steps)
t = (0:999)’; % Make sure t is a column vector

% Initial temperature (20°C)
x0 = 20;

% Heater input (step response: 50% heater power after 100s)
u = zeros(size(t));
u(t >= 100) = 50;

% Parameters [U, alpha, tau]
p = [0.01, 0.01, 200]; % Example values for U, alpha, tau

% Simulate the system
y = tclabsim(t, x0, u, p);
However, every time I try to run it gives the same error.

I have tried to load them as simple double and I have also tried a matrix with t on the first column and the same value in every line of the second column. They also gave the same error.
The block in Simulink has these properties. And I have also tried to change the sample time to -1.

This is the entire simulink if it helps.

How can I fix it?
Thank you simulink, matlab, workspace MATLAB Answers — New Questions

​

I’m working on a MATLAB project where I need to analyze and plot the relative difference in coverage area for different QAM orders and semi-angles (ϕ_1/2) for a wireless communication simulation. I’ve encountered two persistent errors:
Error using xticks (line 33)
Value must be a numeric vector whose values increase.
Error in untitled11 (line 84)
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
and One more thing is the legend has multiples 15, 30, 45 as shown in the below image

Expected Outcome:
Correct x-axis ticks: The x-axis should display the QAM orders {16,32,64,128,256,512,1024}.
Unique legend entries: Each angle (ϕ1/2 =15 ∘ ,30 ∘ ,45 ∘) should appear exactly once in the legend.
Any assistance, please?
here is below my code:

clear variables;
close all;
clc;

%% Define angles and file paths
angles = [15, 30, 45]; % Define angles
angleLabels = ["15°", "30°", "45°"]; % Labels for legends
numAngles = length(angles);

% Load data from the first angle to initialize QAM_order
load([‘workspace_withBeamsteering_’, num2str(angles(1)), ‘.mat’]); % Replace with actual file path
QAM_order = unique(sort(double(QAM_order(:)))); % Ensure QAM_order is numeric, sorted, and unique
numQAMOrders = length(QAM_order); % Dynamically determine number of QAM orders

% Debug: Display QAM_order
disp(‘QAM_order:’);
disp(QAM_order);

% Initialize storage for results
coverageRelativeDiffAll = zeros(numAngles, numQAMOrders); % Initialize storage for relative differences

%% Process data for each angle
for a = 1:numAngles
angle = angles(a);

% Load data for the perfect beamsteering case
load([‘workspace_withBeamsteering_’, num2str(angle), ‘.mat’]); % Replace with actual file path
berMapPerfectCase = berMap;
coveragePerfectCase = coverage;

% Load data for the beamsteering with misalignment case
load([‘workspace_’, num2str(angle), ‘_withBeamsteering_Misalignment.mat’]); % Replace with actual file path

% Compute relative difference in coverage area for this angle
for q = 1:numQAMOrders
% Ensure correct size for coverageMap
[N_directions, N_rx, N_ry] = size(berMap{q});
coverageMap = zeros(N_rx, N_ry);

% Accumulate coverage for all directions
for i_n = 1:N_directions
currentCoverage = squeeze(berMap{q}(i_n, :, :)) < 3.8e-3;
coverageMap = coverageMap + double(currentCoverage);
end

% Compute the intersection map (points covered in all directions)
intersectionMap = (coverageMap == N_directions);

% Calculate misalignment coverage
coverageMisalignment = sum(intersectionMap(:)) * (L * W) / numel(intersectionMap);

% Compute relative difference
coverageRelativeDiffAll(a, q) = (coverageMisalignment – coveragePerfectCase(q)) / coveragePerfectCase(q);
end
end

%% Plot Relative Difference in Coverage Area for All Angles
figure;
hold on;

% Define markers and colors for each angle
markers = {‘-o’, ‘-s’, ‘-^’}; % Different markers for angles
colors = {‘b’, ‘r’, ‘g’}; % Blue, Red, Green

for a = 1:numAngles
plot(QAM_order, coverageRelativeDiffAll(a, 🙂 * 100, …
markers{a}, ‘LineWidth’, 1.5, ‘Color’, colors{a}, …
‘DisplayName’, [‘$phi_{1/2} = $ ‘, num2str(angles(a)), ‘$^circ$’]);
end

hold off;

% Add labels, title, and legend
xlabel(‘QAM Order (M)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
ylabel(‘Relative Difference in Coverage Area (%)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
title(‘Relative Difference in Coverage Area vs QAM Order’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
legend(‘Location’, ‘northeast’, ‘FontSize’, 11, ‘Interpreter’, ‘latex’);

% Set grid and font appearance
grid on;
set(gca, ‘FontSize’, 11, ‘FontWeight’, ‘bold’);

% Set x-axis to show QAM order directly
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
xticklabels(string(QAM_order)); % Display QAM order values directly

% Set x-axis limits to start at the smallest QAM order
xlim([min(QAM_order), max(QAM_order)]); % Start from the minimum and go to the maximum QAM orderI’m working on a MATLAB project where I need to analyze and plot the relative difference in coverage area for different QAM orders and semi-angles (ϕ_1/2) for a wireless communication simulation. I’ve encountered two persistent errors:
Error using xticks (line 33)
Value must be a numeric vector whose values increase.
Error in untitled11 (line 84)
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
and One more thing is the legend has multiples 15, 30, 45 as shown in the below image

Expected Outcome:
Correct x-axis ticks: The x-axis should display the QAM orders {16,32,64,128,256,512,1024}.
Unique legend entries: Each angle (ϕ1/2 =15 ∘ ,30 ∘ ,45 ∘) should appear exactly once in the legend.
Any assistance, please?
here is below my code:

clear variables;
close all;
clc;

%% Define angles and file paths
angles = [15, 30, 45]; % Define angles
angleLabels = ["15°", "30°", "45°"]; % Labels for legends
numAngles = length(angles);

% Load data from the first angle to initialize QAM_order
load([‘workspace_withBeamsteering_’, num2str(angles(1)), ‘.mat’]); % Replace with actual file path
QAM_order = unique(sort(double(QAM_order(:)))); % Ensure QAM_order is numeric, sorted, and unique
numQAMOrders = length(QAM_order); % Dynamically determine number of QAM orders

% Debug: Display QAM_order
disp(‘QAM_order:’);
disp(QAM_order);

% Initialize storage for results
coverageRelativeDiffAll = zeros(numAngles, numQAMOrders); % Initialize storage for relative differences

%% Process data for each angle
for a = 1:numAngles
angle = angles(a);

% Load data for the perfect beamsteering case
load([‘workspace_withBeamsteering_’, num2str(angle), ‘.mat’]); % Replace with actual file path
berMapPerfectCase = berMap;
coveragePerfectCase = coverage;

% Load data for the beamsteering with misalignment case
load([‘workspace_’, num2str(angle), ‘_withBeamsteering_Misalignment.mat’]); % Replace with actual file path

% Compute relative difference in coverage area for this angle
for q = 1:numQAMOrders
% Ensure correct size for coverageMap
[N_directions, N_rx, N_ry] = size(berMap{q});
coverageMap = zeros(N_rx, N_ry);

% Accumulate coverage for all directions
for i_n = 1:N_directions
currentCoverage = squeeze(berMap{q}(i_n, :, :)) < 3.8e-3;
coverageMap = coverageMap + double(currentCoverage);
end

% Compute the intersection map (points covered in all directions)
intersectionMap = (coverageMap == N_directions);

% Calculate misalignment coverage
coverageMisalignment = sum(intersectionMap(:)) * (L * W) / numel(intersectionMap);

% Compute relative difference
coverageRelativeDiffAll(a, q) = (coverageMisalignment – coveragePerfectCase(q)) / coveragePerfectCase(q);
end
end

%% Plot Relative Difference in Coverage Area for All Angles
figure;
hold on;

% Define markers and colors for each angle
markers = {‘-o’, ‘-s’, ‘-^’}; % Different markers for angles
colors = {‘b’, ‘r’, ‘g’}; % Blue, Red, Green

for a = 1:numAngles
plot(QAM_order, coverageRelativeDiffAll(a, 🙂 * 100, …
markers{a}, ‘LineWidth’, 1.5, ‘Color’, colors{a}, …
‘DisplayName’, [‘$phi_{1/2} = $ ‘, num2str(angles(a)), ‘$^circ$’]);
end

hold off;

% Add labels, title, and legend
xlabel(‘QAM Order (M)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
ylabel(‘Relative Difference in Coverage Area (%)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
title(‘Relative Difference in Coverage Area vs QAM Order’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
legend(‘Location’, ‘northeast’, ‘FontSize’, 11, ‘Interpreter’, ‘latex’);

% Set grid and font appearance
grid on;
set(gca, ‘FontSize’, 11, ‘FontWeight’, ‘bold’);

% Set x-axis to show QAM order directly
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
xticklabels(string(QAM_order)); % Display QAM order values directly

% Set x-axis limits to start at the smallest QAM order
xlim([min(QAM_order), max(QAM_order)]); % Start from the minimum and go to the maximum QAM order I’m working on a MATLAB project where I need to analyze and plot the relative difference in coverage area for different QAM orders and semi-angles (ϕ_1/2) for a wireless communication simulation. I’ve encountered two persistent errors:
Error using xticks (line 33)
Value must be a numeric vector whose values increase.
Error in untitled11 (line 84)
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
and One more thing is the legend has multiples 15, 30, 45 as shown in the below image

Expected Outcome:
Correct x-axis ticks: The x-axis should display the QAM orders {16,32,64,128,256,512,1024}.
Unique legend entries: Each angle (ϕ1/2 =15 ∘ ,30 ∘ ,45 ∘) should appear exactly once in the legend.
Any assistance, please?
here is below my code:

clear variables;
close all;
clc;

%% Define angles and file paths
angles = [15, 30, 45]; % Define angles
angleLabels = ["15°", "30°", "45°"]; % Labels for legends
numAngles = length(angles);

% Load data from the first angle to initialize QAM_order
load([‘workspace_withBeamsteering_’, num2str(angles(1)), ‘.mat’]); % Replace with actual file path
QAM_order = unique(sort(double(QAM_order(:)))); % Ensure QAM_order is numeric, sorted, and unique
numQAMOrders = length(QAM_order); % Dynamically determine number of QAM orders

% Debug: Display QAM_order
disp(‘QAM_order:’);
disp(QAM_order);

% Initialize storage for results
coverageRelativeDiffAll = zeros(numAngles, numQAMOrders); % Initialize storage for relative differences

%% Process data for each angle
for a = 1:numAngles
angle = angles(a);

% Load data for the perfect beamsteering case
load([‘workspace_withBeamsteering_’, num2str(angle), ‘.mat’]); % Replace with actual file path
berMapPerfectCase = berMap;
coveragePerfectCase = coverage;

% Load data for the beamsteering with misalignment case
load([‘workspace_’, num2str(angle), ‘_withBeamsteering_Misalignment.mat’]); % Replace with actual file path

% Compute relative difference in coverage area for this angle
for q = 1:numQAMOrders
% Ensure correct size for coverageMap
[N_directions, N_rx, N_ry] = size(berMap{q});
coverageMap = zeros(N_rx, N_ry);

% Accumulate coverage for all directions
for i_n = 1:N_directions
currentCoverage = squeeze(berMap{q}(i_n, :, :)) < 3.8e-3;
coverageMap = coverageMap + double(currentCoverage);
end

% Compute the intersection map (points covered in all directions)
intersectionMap = (coverageMap == N_directions);

% Calculate misalignment coverage
coverageMisalignment = sum(intersectionMap(:)) * (L * W) / numel(intersectionMap);

% Compute relative difference
coverageRelativeDiffAll(a, q) = (coverageMisalignment – coveragePerfectCase(q)) / coveragePerfectCase(q);
end
end

%% Plot Relative Difference in Coverage Area for All Angles
figure;
hold on;

% Define markers and colors for each angle
markers = {‘-o’, ‘-s’, ‘-^’}; % Different markers for angles
colors = {‘b’, ‘r’, ‘g’}; % Blue, Red, Green

for a = 1:numAngles
plot(QAM_order, coverageRelativeDiffAll(a, 🙂 * 100, …
markers{a}, ‘LineWidth’, 1.5, ‘Color’, colors{a}, …
‘DisplayName’, [‘$phi_{1/2} = $ ‘, num2str(angles(a)), ‘$^circ$’]);
end

hold off;

% Add labels, title, and legend
xlabel(‘QAM Order (M)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
ylabel(‘Relative Difference in Coverage Area (%)’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
title(‘Relative Difference in Coverage Area vs QAM Order’, ‘FontSize’, 12, ‘FontWeight’, ‘bold’);
legend(‘Location’, ‘northeast’, ‘FontSize’, 11, ‘Interpreter’, ‘latex’);

% Set grid and font appearance
grid on;
set(gca, ‘FontSize’, 11, ‘FontWeight’, ‘bold’);

% Set x-axis to show QAM order directly
xticks(QAM_order); % Ensure QAM_order is numeric and sorted
xticklabels(string(QAM_order)); % Display QAM order values directly

% Set x-axis limits to start at the smallest QAM order
xlim([min(QAM_order), max(QAM_order)]); % Start from the minimum and go to the maximum QAM order plot, xticks, legend MATLAB Answers — New Questions

​