I want to use a hardware interrupt, in particular, the TXE flag on Simulink, however I have found that there does not appear to be a hardware interrupt block in the simulink support package for stm32 nucleo boards.

The hardware interrupt blocks for STM32F4 Discovery and the ARM Cortex-Mx don’t give any errors but I cannot imagine that using them would be correct.

Is there a way to get the TXE flag for the STM32 Nucleo F767ZI on simulink or should I just restart my project on the STM32CubeIDE instead.I want to use a hardware interrupt, in particular, the TXE flag on Simulink, however I have found that there does not appear to be a hardware interrupt block in the simulink support package for stm32 nucleo boards.

The hardware interrupt blocks for STM32F4 Discovery and the ARM Cortex-Mx don’t give any errors but I cannot imagine that using them would be correct.

Is there a way to get the TXE flag for the STM32 Nucleo F767ZI on simulink or should I just restart my project on the STM32CubeIDE instead. I want to use a hardware interrupt, in particular, the TXE flag on Simulink, however I have found that there does not appear to be a hardware interrupt block in the simulink support package for stm32 nucleo boards.

The hardware interrupt blocks for STM32F4 Discovery and the ARM Cortex-Mx don’t give any errors but I cannot imagine that using them would be correct.

Is there a way to get the TXE flag for the STM32 Nucleo F767ZI on simulink or should I just restart my project on the STM32CubeIDE instead. simulink, embedded coder MATLAB Answers — New Questions

​

Hi there! What is the best way to plot a 3D surface whose matrix C = A/B contains values of +/- infinity (whenever entries in B are zero)? Currently I am using surf( ), which results in weird vertical cross sections. I looked into the fillMissing( ) and scatteredInterpolant( ) functions. Is there a best / recommended way to proceed, to plot a smooth surface? Thanks in advance,Hi there! What is the best way to plot a 3D surface whose matrix C = A/B contains values of +/- infinity (whenever entries in B are zero)? Currently I am using surf( ), which results in weird vertical cross sections. I looked into the fillMissing( ) and scatteredInterpolant( ) functions. Is there a best / recommended way to proceed, to plot a smooth surface? Thanks in advance, Hi there! What is the best way to plot a 3D surface whose matrix C = A/B contains values of +/- infinity (whenever entries in B are zero)? Currently I am using surf( ), which results in weird vertical cross sections. I looked into the fillMissing( ) and scatteredInterpolant( ) functions. Is there a best / recommended way to proceed, to plot a smooth surface? Thanks in advance, matlab surface plot singularities MATLAB Answers — New Questions

​

The way my CI is currently setup is that it creates a new project per branch. What I want is to compare the result of two branches using the polyspace-access export command. I tried downloading the result from one project and uploading it to the other but it failed.The way my CI is currently setup is that it creates a new project per branch. What I want is to compare the result of two branches using the polyspace-access export command. I tried downloading the result from one project and uploading it to the other but it failed. The way my CI is currently setup is that it creates a new project per branch. What I want is to compare the result of two branches using the polyspace-access export command. I tried downloading the result from one project and uploading it to the other but it failed. bugfinder, polyspace_access, polyspace_results MATLAB Answers — New Questions

​

I wish to understand the DTFT and DFT using MATLAB.
So I want to write a code which does the following:
plot X1[n]=COS(2*pi*201/1024*n) where n=0 to1023
plot X2[n]=rect(N) where N=1024 i.e. X2[n]=1 for n=0 to 1023; else X2[n]=0
plot multiplication of both i.e. X1(n)*X2(n)
plot magnitude plot of Y1(w)=DTFT of X1(n) for -pi<w<pi
plot magnitude plot of Y2(w)=DTFT of X2(n) for -pi<w<pi
plot convolution of Y1 and Y2 i.e. Z1(w)=Y1(w) ⊗ Y2(w)
plot frequency samples Z1[k] that are values of Z1(w) at w[k]=2*pi*k/N where 0<k<(N-1)
Plz help me with the codeI wish to understand the DTFT and DFT using MATLAB.
So I want to write a code which does the following:
plot X1[n]=COS(2*pi*201/1024*n) where n=0 to1023
plot X2[n]=rect(N) where N=1024 i.e. X2[n]=1 for n=0 to 1023; else X2[n]=0
plot multiplication of both i.e. X1(n)*X2(n)
plot magnitude plot of Y1(w)=DTFT of X1(n) for -pi<w<pi
plot magnitude plot of Y2(w)=DTFT of X2(n) for -pi<w<pi
plot convolution of Y1 and Y2 i.e. Z1(w)=Y1(w) ⊗ Y2(w)
plot frequency samples Z1[k] that are values of Z1(w) at w[k]=2*pi*k/N where 0<k<(N-1)
Plz help me with the code I wish to understand the DTFT and DFT using MATLAB.
So I want to write a code which does the following:
plot X1[n]=COS(2*pi*201/1024*n) where n=0 to1023
plot X2[n]=rect(N) where N=1024 i.e. X2[n]=1 for n=0 to 1023; else X2[n]=0
plot multiplication of both i.e. X1(n)*X2(n)
plot magnitude plot of Y1(w)=DTFT of X1(n) for -pi<w<pi
plot magnitude plot of Y2(w)=DTFT of X2(n) for -pi<w<pi
plot convolution of Y1 and Y2 i.e. Z1(w)=Y1(w) ⊗ Y2(w)
plot frequency samples Z1[k] that are values of Z1(w) at w[k]=2*pi*k/N where 0<k<(N-1)
Plz help me with the code dft, dtft, freqz, fft, matlab code MATLAB Answers — New Questions

​

Hi
i find the scattering field in cylinder and i use the MOM method with the triangle basis and pulse. The |Es
| are the same when then incoming field are Einc= E0*exp(-j*k0*r*cos(φ)) in MOM pulse and Einc= E0*exp(j*k0*r*cos(φ)) in MOM triangle.
This the pulse current
function [Is] = currentMoM()
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
% Constants
gamma_const = 1.781;
dfpm = 2*pi/N;
Phi0 = (0:N-1) * dfpm;
% Precompute pairwise phase differences
deltaPhi = zeros(N, N); % Initialize deltaPhi matrix
for i = 1:N
for j = 1:N
deltaPhi(i, j) = Phi0(i) – Phi0(j); % Calculate pairwise phase differences
end
end
% Compute distance matrix R (N x N matrix)
R = zeros(N, N); % Initialize R matrix
for i = 1:N
for j = 1:N
R(i, j) = ra * sqrt(2 – 2 * cos(deltaPhi(i, j))); % Compute distance matrix R
end
end
% Coefficients
coeif1 = (Z0 * k0 / 4) * ra * dfpm;
coeif2 = (Z0 / 2) * sin(k0 * ra * dfpm/ 2);
coeif3 = (gamma_const * k0 * ra * dfpm) / 4;
% Compute lmn matrix
lmn = zeros(N, N); % Initialize lmn matrix
for i = 1:N
for j = 1:N
if i==j
lmn(i, j) = coeif1 * (1 – 1i * (2 / pi) * (log(coeif3) – 1)); % Diagonal elements
else
lmn(i, j) = coeif2 * besselh(0, 2, k0 * R(i, j)); % Compute each element of lmn
end
end
zmn = lmn;
end
% Handle diagonal elements separately

% Set zmn equal to lmn

% Compute gm vector using a for loop
gm = zeros(N, 1); % Initialize gm vector
for i = 1:N
gm(i) = E_in_z(Phi0(i)); % Calculate each element of gm using E_in_z function
end
% Solve for currents using linsolve
%Is = linsolve(zmn, gm);
Is = zmn gm;
%Is=linsolve(zmn, gm)
end
and the triangke
function [Is] = currentMoM()
% Parameter initialization
[~, N, ra, k0, Z0, ~] = parameter();
gamma_const = 1.781;
Eps = 1e-9;
dftm =2.*pi./N;

% Precompute constants
B = (4./(Z0 * k0));
A = (gamma_const*k0*ra)./2;

% Initialize matrices and vectors
lmn = zeros(N, N); % (N x N) matrix for lmn
gm = zeros(1, N); % (1 x N) vector for gm
zmn = zeros(N, N); % (N x N) matrix for zmn

% Compute lmn and zmn matrices using nested for-loops
for index_i = 1:N
for index_j = 1:N
% Precompute integration bounds
xi1_i = (index_i – 1) * dftm;
xi2_i = xi1_i + dftm;
xj1_j = (index_j – 1) * dftm;
xj2_j = xj1_j + dftm;
%log(A) +log(abs(x – y) +Eps))
if index_i == index_j
% Diagonal elements
funa = @(x, y) triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* (1 – 1i * (2./pi).*(log(A) +log(abs(x – y) +Eps))) ;
lmn(index_i, index_j) = integral2(funa, xi1_i, xi2_i, xj1_j, xj2_j,’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’Method’,’iterated’);
else
% Off-diagonal elements
funb = @(x, y)triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* besselh(0, 2, k0.*ra.*sqrt(abs(2 – 2 * cos(x – y))) );
lmn(index_i, index_j) = integral2(funb, xi1_i, xi2_i, xj1_j, xj2_j, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Assign to zmn
zmn(index_i, index_j) = lmn(index_i, index_j);
end

end

% Compute gm vector using a single for-loop
for index_i = 1:N
xi1 = (index_i – 1) * dftm;
xi2 = xi1 + dftm;
func = @(x)B.* triangle_basisn(x, index_i) .* (E_in_z(x));
gm(index_i) = integral(func, xi1, xi2, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Solve the linear system
epsilon_reg = 1e-10; % Regularization parameter (can be adjusted)

zmn = zmn + epsilon_reg * eye(N); % Add small value to the diagonal
% Solve the linear system
Is = linsolve(zmn, gm’);

end

the scattering field
function z = Escat(r, phi)
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
[Is] = currentMoM(); % Retrieve the current distribution

% Preallocate for final result
FINAL = zeros(size(r)); % Initialize FINAL to the same size as r

dfpulse = 2 * pi / N;
Phi0 = (0:N-1) * dfpulse; % Angle steps (in radians)
coif = (k0 * Z0 / 4); % Coefficient to scale the results
% Create a column vector of Phi0 values
x1 = Phi0(:);
x2 = x1 + dfpulse; % Integration limits (x1, x2)

% Loop to compute the scattered field for each combination of rho and phi
for jj = 1:N
% Ensure phi and r are compatible in size for broadcasting

r_squared = r.^2; % Compute r squared
ra_squared = ra.^2;
% Compute the integrand term

% Define the integrand with the properly dimensioned term
integrand = @(y)besselh(0, 2, k0 * sqrt(r_squared + ra_squared – 2*r.*ra .*cos(phi-y)));

% Perform the numerical integration for each jj
FINAL = FINAL -coif * Is(jj).* ra .* integral(integrand, x1(jj), x2(jj),’ArrayValued’,true);
end

% Return the result
z = FINAL;
end
function z = Escat(r,phi)
[~,N,ra,k0,Z0,~] =parameter();
Is=currentMoM();
%Phi0=zeros(N);
FINAL=0;
A=(k0.*Z0./4);
dftm=2.*pi./N;
%index=1:N
%Phi0=(index-2)*dftm;
Phi0=(0:N-1).*dftm;
x1 = Phi0(:); % Create a column vector of x1 values
x2 = x1 + dftm; % x2 values depend on x1
for jj=1:N
F=@(y)besselh(0,2,k0*sqrt(r.^2 + ra.^2-2.*r.*ra.*cos(phi-y))).* triangle_basisn(y,jj);

% Perform the integral
%pts = optimset(‘TolFun’,1e-6); % Adjust tolerance
%integral_result = integral(integrand, lower_limit, upper_limit, ‘ArrayValued’, true);

integral_result=integral(F,x1(jj),x2(jj),’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’ArrayValued’,false);
% Accumulate the result
FINAL = FINAL-A.*Is(jj).*ra.*integral_result;
end
z=FINAL;
end
function [f,N,ra,k0,Z0,lambda] = parameter()
%UNTITLED Summary of this function goes here
c0=3e8;
Z0=120.*pi;
ra=1;
N=60;
f=300e6;
lambda=c0./f;
k0=2*pi./lambda;
end
and teh triangle basis function
function z = triangle_basisn(phi, k)
% Get number of segments (N) and compute segment width (dftm)
[~, N, ~, ~, ~, ~] = parameter();
df = 2*pi/ N; % Segment width

% Compute the left and right boundaries for the triangle basis
left_bound = (k – 2) * df;
right_bound = (k – 1) * df;
center_bound = k * df;

% Check if phi lies within the support of the triangular basis function
if phi >= left_bound & phi <= right_bound
% Left part of the triangle
z = (phi – left_bound) / df;
elseif phi > right_bound & phi <= center_bound
% Right part of the triangle
z = (center_bound – phi) / df;
else
% Outside the support of the triangle
z = 0;
end
end

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(+j.*k0.*ra.*cos(phi));
end
for the triangle

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(-j.*k0.*ra.*cos(phi));
end
for the pulse

thank you
George

the main script is

clear all
clc
[f,N,ra,k0,Z0,lambda] = parameter();
%ph_i=pi/2;
rho=6*lambda ;
phi=-pi:pi/180:pi;
Es=zeros(size(phi));
phi_d=zeros(size(phi));
%phi=0:pi/200:pi/3
tic
parfor jj=1:length(phi)
Es(jj)=abs(Escat(rho,phi(jj)));
phi_d(jj)=phi(jj).*(180/pi);
end

plot(phi_d,Es,’b–‘,’LineWidth’,2)Hi
i find the scattering field in cylinder and i use the MOM method with the triangle basis and pulse. The |Es
| are the same when then incoming field are Einc= E0*exp(-j*k0*r*cos(φ)) in MOM pulse and Einc= E0*exp(j*k0*r*cos(φ)) in MOM triangle.
This the pulse current
function [Is] = currentMoM()
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
% Constants
gamma_const = 1.781;
dfpm = 2*pi/N;
Phi0 = (0:N-1) * dfpm;
% Precompute pairwise phase differences
deltaPhi = zeros(N, N); % Initialize deltaPhi matrix
for i = 1:N
for j = 1:N
deltaPhi(i, j) = Phi0(i) – Phi0(j); % Calculate pairwise phase differences
end
end
% Compute distance matrix R (N x N matrix)
R = zeros(N, N); % Initialize R matrix
for i = 1:N
for j = 1:N
R(i, j) = ra * sqrt(2 – 2 * cos(deltaPhi(i, j))); % Compute distance matrix R
end
end
% Coefficients
coeif1 = (Z0 * k0 / 4) * ra * dfpm;
coeif2 = (Z0 / 2) * sin(k0 * ra * dfpm/ 2);
coeif3 = (gamma_const * k0 * ra * dfpm) / 4;
% Compute lmn matrix
lmn = zeros(N, N); % Initialize lmn matrix
for i = 1:N
for j = 1:N
if i==j
lmn(i, j) = coeif1 * (1 – 1i * (2 / pi) * (log(coeif3) – 1)); % Diagonal elements
else
lmn(i, j) = coeif2 * besselh(0, 2, k0 * R(i, j)); % Compute each element of lmn
end
end
zmn = lmn;
end
% Handle diagonal elements separately

% Set zmn equal to lmn

% Compute gm vector using a for loop
gm = zeros(N, 1); % Initialize gm vector
for i = 1:N
gm(i) = E_in_z(Phi0(i)); % Calculate each element of gm using E_in_z function
end
% Solve for currents using linsolve
%Is = linsolve(zmn, gm);
Is = zmn gm;
%Is=linsolve(zmn, gm)
end
and the triangke
function [Is] = currentMoM()
% Parameter initialization
[~, N, ra, k0, Z0, ~] = parameter();
gamma_const = 1.781;
Eps = 1e-9;
dftm =2.*pi./N;

% Precompute constants
B = (4./(Z0 * k0));
A = (gamma_const*k0*ra)./2;

% Initialize matrices and vectors
lmn = zeros(N, N); % (N x N) matrix for lmn
gm = zeros(1, N); % (1 x N) vector for gm
zmn = zeros(N, N); % (N x N) matrix for zmn

% Compute lmn and zmn matrices using nested for-loops
for index_i = 1:N
for index_j = 1:N
% Precompute integration bounds
xi1_i = (index_i – 1) * dftm;
xi2_i = xi1_i + dftm;
xj1_j = (index_j – 1) * dftm;
xj2_j = xj1_j + dftm;
%log(A) +log(abs(x – y) +Eps))
if index_i == index_j
% Diagonal elements
funa = @(x, y) triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* (1 – 1i * (2./pi).*(log(A) +log(abs(x – y) +Eps))) ;
lmn(index_i, index_j) = integral2(funa, xi1_i, xi2_i, xj1_j, xj2_j,’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’Method’,’iterated’);
else
% Off-diagonal elements
funb = @(x, y)triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* besselh(0, 2, k0.*ra.*sqrt(abs(2 – 2 * cos(x – y))) );
lmn(index_i, index_j) = integral2(funb, xi1_i, xi2_i, xj1_j, xj2_j, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Assign to zmn
zmn(index_i, index_j) = lmn(index_i, index_j);
end

end

% Compute gm vector using a single for-loop
for index_i = 1:N
xi1 = (index_i – 1) * dftm;
xi2 = xi1 + dftm;
func = @(x)B.* triangle_basisn(x, index_i) .* (E_in_z(x));
gm(index_i) = integral(func, xi1, xi2, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Solve the linear system
epsilon_reg = 1e-10; % Regularization parameter (can be adjusted)

zmn = zmn + epsilon_reg * eye(N); % Add small value to the diagonal
% Solve the linear system
Is = linsolve(zmn, gm’);

end

the scattering field
function z = Escat(r, phi)
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
[Is] = currentMoM(); % Retrieve the current distribution

% Preallocate for final result
FINAL = zeros(size(r)); % Initialize FINAL to the same size as r

dfpulse = 2 * pi / N;
Phi0 = (0:N-1) * dfpulse; % Angle steps (in radians)
coif = (k0 * Z0 / 4); % Coefficient to scale the results
% Create a column vector of Phi0 values
x1 = Phi0(:);
x2 = x1 + dfpulse; % Integration limits (x1, x2)

% Loop to compute the scattered field for each combination of rho and phi
for jj = 1:N
% Ensure phi and r are compatible in size for broadcasting

r_squared = r.^2; % Compute r squared
ra_squared = ra.^2;
% Compute the integrand term

% Define the integrand with the properly dimensioned term
integrand = @(y)besselh(0, 2, k0 * sqrt(r_squared + ra_squared – 2*r.*ra .*cos(phi-y)));

% Perform the numerical integration for each jj
FINAL = FINAL -coif * Is(jj).* ra .* integral(integrand, x1(jj), x2(jj),’ArrayValued’,true);
end

% Return the result
z = FINAL;
end
function z = Escat(r,phi)
[~,N,ra,k0,Z0,~] =parameter();
Is=currentMoM();
%Phi0=zeros(N);
FINAL=0;
A=(k0.*Z0./4);
dftm=2.*pi./N;
%index=1:N
%Phi0=(index-2)*dftm;
Phi0=(0:N-1).*dftm;
x1 = Phi0(:); % Create a column vector of x1 values
x2 = x1 + dftm; % x2 values depend on x1
for jj=1:N
F=@(y)besselh(0,2,k0*sqrt(r.^2 + ra.^2-2.*r.*ra.*cos(phi-y))).* triangle_basisn(y,jj);

% Perform the integral
%pts = optimset(‘TolFun’,1e-6); % Adjust tolerance
%integral_result = integral(integrand, lower_limit, upper_limit, ‘ArrayValued’, true);

integral_result=integral(F,x1(jj),x2(jj),’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’ArrayValued’,false);
% Accumulate the result
FINAL = FINAL-A.*Is(jj).*ra.*integral_result;
end
z=FINAL;
end
function [f,N,ra,k0,Z0,lambda] = parameter()
%UNTITLED Summary of this function goes here
c0=3e8;
Z0=120.*pi;
ra=1;
N=60;
f=300e6;
lambda=c0./f;
k0=2*pi./lambda;
end
and teh triangle basis function
function z = triangle_basisn(phi, k)
% Get number of segments (N) and compute segment width (dftm)
[~, N, ~, ~, ~, ~] = parameter();
df = 2*pi/ N; % Segment width

% Compute the left and right boundaries for the triangle basis
left_bound = (k – 2) * df;
right_bound = (k – 1) * df;
center_bound = k * df;

% Check if phi lies within the support of the triangular basis function
if phi >= left_bound & phi <= right_bound
% Left part of the triangle
z = (phi – left_bound) / df;
elseif phi > right_bound & phi <= center_bound
% Right part of the triangle
z = (center_bound – phi) / df;
else
% Outside the support of the triangle
z = 0;
end
end

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(+j.*k0.*ra.*cos(phi));
end
for the triangle

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(-j.*k0.*ra.*cos(phi));
end
for the pulse

thank you
George

the main script is

clear all
clc
[f,N,ra,k0,Z0,lambda] = parameter();
%ph_i=pi/2;
rho=6*lambda ;
phi=-pi:pi/180:pi;
Es=zeros(size(phi));
phi_d=zeros(size(phi));
%phi=0:pi/200:pi/3
tic
parfor jj=1:length(phi)
Es(jj)=abs(Escat(rho,phi(jj)));
phi_d(jj)=phi(jj).*(180/pi);
end

plot(phi_d,Es,’b–‘,’LineWidth’,2) Hi
i find the scattering field in cylinder and i use the MOM method with the triangle basis and pulse. The |Es
| are the same when then incoming field are Einc= E0*exp(-j*k0*r*cos(φ)) in MOM pulse and Einc= E0*exp(j*k0*r*cos(φ)) in MOM triangle.
This the pulse current
function [Is] = currentMoM()
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
% Constants
gamma_const = 1.781;
dfpm = 2*pi/N;
Phi0 = (0:N-1) * dfpm;
% Precompute pairwise phase differences
deltaPhi = zeros(N, N); % Initialize deltaPhi matrix
for i = 1:N
for j = 1:N
deltaPhi(i, j) = Phi0(i) – Phi0(j); % Calculate pairwise phase differences
end
end
% Compute distance matrix R (N x N matrix)
R = zeros(N, N); % Initialize R matrix
for i = 1:N
for j = 1:N
R(i, j) = ra * sqrt(2 – 2 * cos(deltaPhi(i, j))); % Compute distance matrix R
end
end
% Coefficients
coeif1 = (Z0 * k0 / 4) * ra * dfpm;
coeif2 = (Z0 / 2) * sin(k0 * ra * dfpm/ 2);
coeif3 = (gamma_const * k0 * ra * dfpm) / 4;
% Compute lmn matrix
lmn = zeros(N, N); % Initialize lmn matrix
for i = 1:N
for j = 1:N
if i==j
lmn(i, j) = coeif1 * (1 – 1i * (2 / pi) * (log(coeif3) – 1)); % Diagonal elements
else
lmn(i, j) = coeif2 * besselh(0, 2, k0 * R(i, j)); % Compute each element of lmn
end
end
zmn = lmn;
end
% Handle diagonal elements separately

% Set zmn equal to lmn

% Compute gm vector using a for loop
gm = zeros(N, 1); % Initialize gm vector
for i = 1:N
gm(i) = E_in_z(Phi0(i)); % Calculate each element of gm using E_in_z function
end
% Solve for currents using linsolve
%Is = linsolve(zmn, gm);
Is = zmn gm;
%Is=linsolve(zmn, gm)
end
and the triangke
function [Is] = currentMoM()
% Parameter initialization
[~, N, ra, k0, Z0, ~] = parameter();
gamma_const = 1.781;
Eps = 1e-9;
dftm =2.*pi./N;

% Precompute constants
B = (4./(Z0 * k0));
A = (gamma_const*k0*ra)./2;

% Initialize matrices and vectors
lmn = zeros(N, N); % (N x N) matrix for lmn
gm = zeros(1, N); % (1 x N) vector for gm
zmn = zeros(N, N); % (N x N) matrix for zmn

% Compute lmn and zmn matrices using nested for-loops
for index_i = 1:N
for index_j = 1:N
% Precompute integration bounds
xi1_i = (index_i – 1) * dftm;
xi2_i = xi1_i + dftm;
xj1_j = (index_j – 1) * dftm;
xj2_j = xj1_j + dftm;
%log(A) +log(abs(x – y) +Eps))
if index_i == index_j
% Diagonal elements
funa = @(x, y) triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* (1 – 1i * (2./pi).*(log(A) +log(abs(x – y) +Eps))) ;
lmn(index_i, index_j) = integral2(funa, xi1_i, xi2_i, xj1_j, xj2_j,’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’Method’,’iterated’);
else
% Off-diagonal elements
funb = @(x, y)triangle_basisn(x, index_i) .* triangle_basisn(y, index_j) .*ra .* besselh(0, 2, k0.*ra.*sqrt(abs(2 – 2 * cos(x – y))) );
lmn(index_i, index_j) = integral2(funb, xi1_i, xi2_i, xj1_j, xj2_j, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Assign to zmn
zmn(index_i, index_j) = lmn(index_i, index_j);
end

end

% Compute gm vector using a single for-loop
for index_i = 1:N
xi1 = (index_i – 1) * dftm;
xi2 = xi1 + dftm;
func = @(x)B.* triangle_basisn(x, index_i) .* (E_in_z(x));
gm(index_i) = integral(func, xi1, xi2, ‘RelTol’, 1e-6, ‘AbsTol’, 1e-6);
end

% Solve the linear system
epsilon_reg = 1e-10; % Regularization parameter (can be adjusted)

zmn = zmn + epsilon_reg * eye(N); % Add small value to the diagonal
% Solve the linear system
Is = linsolve(zmn, gm’);

end

the scattering field
function z = Escat(r, phi)
% Load parameters
[f, N, ra, k0, Z0, lambda] = parameter();
[Is] = currentMoM(); % Retrieve the current distribution

% Preallocate for final result
FINAL = zeros(size(r)); % Initialize FINAL to the same size as r

dfpulse = 2 * pi / N;
Phi0 = (0:N-1) * dfpulse; % Angle steps (in radians)
coif = (k0 * Z0 / 4); % Coefficient to scale the results
% Create a column vector of Phi0 values
x1 = Phi0(:);
x2 = x1 + dfpulse; % Integration limits (x1, x2)

% Loop to compute the scattered field for each combination of rho and phi
for jj = 1:N
% Ensure phi and r are compatible in size for broadcasting

r_squared = r.^2; % Compute r squared
ra_squared = ra.^2;
% Compute the integrand term

% Define the integrand with the properly dimensioned term
integrand = @(y)besselh(0, 2, k0 * sqrt(r_squared + ra_squared – 2*r.*ra .*cos(phi-y)));

% Perform the numerical integration for each jj
FINAL = FINAL -coif * Is(jj).* ra .* integral(integrand, x1(jj), x2(jj),’ArrayValued’,true);
end

% Return the result
z = FINAL;
end
function z = Escat(r,phi)
[~,N,ra,k0,Z0,~] =parameter();
Is=currentMoM();
%Phi0=zeros(N);
FINAL=0;
A=(k0.*Z0./4);
dftm=2.*pi./N;
%index=1:N
%Phi0=(index-2)*dftm;
Phi0=(0:N-1).*dftm;
x1 = Phi0(:); % Create a column vector of x1 values
x2 = x1 + dftm; % x2 values depend on x1
for jj=1:N
F=@(y)besselh(0,2,k0*sqrt(r.^2 + ra.^2-2.*r.*ra.*cos(phi-y))).* triangle_basisn(y,jj);

% Perform the integral
%pts = optimset(‘TolFun’,1e-6); % Adjust tolerance
%integral_result = integral(integrand, lower_limit, upper_limit, ‘ArrayValued’, true);

integral_result=integral(F,x1(jj),x2(jj),’RelTol’, 1e-6, ‘AbsTol’, 1e-6,’ArrayValued’,false);
% Accumulate the result
FINAL = FINAL-A.*Is(jj).*ra.*integral_result;
end
z=FINAL;
end
function [f,N,ra,k0,Z0,lambda] = parameter()
%UNTITLED Summary of this function goes here
c0=3e8;
Z0=120.*pi;
ra=1;
N=60;
f=300e6;
lambda=c0./f;
k0=2*pi./lambda;
end
and teh triangle basis function
function z = triangle_basisn(phi, k)
% Get number of segments (N) and compute segment width (dftm)
[~, N, ~, ~, ~, ~] = parameter();
df = 2*pi/ N; % Segment width

% Compute the left and right boundaries for the triangle basis
left_bound = (k – 2) * df;
right_bound = (k – 1) * df;
center_bound = k * df;

% Check if phi lies within the support of the triangular basis function
if phi >= left_bound & phi <= right_bound
% Left part of the triangle
z = (phi – left_bound) / df;
elseif phi > right_bound & phi <= center_bound
% Right part of the triangle
z = (center_bound – phi) / df;
else
% Outside the support of the triangle
z = 0;
end
end

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(+j.*k0.*ra.*cos(phi));
end
for the triangle

function z=E_in_z(phi)
%amplitude
[f,N,ra,k0,Z0] = parameter();
E0=1;
%z=E0.*exp(j.*k0.*ra.*cos(phi)+j.*k0.*ra.*sin(phi));
z=E0.*exp(-j.*k0.*ra.*cos(phi));
end
for the pulse

thank you
George

the main script is

clear all
clc
[f,N,ra,k0,Z0,lambda] = parameter();
%ph_i=pi/2;
rho=6*lambda ;
phi=-pi:pi/180:pi;
Es=zeros(size(phi));
phi_d=zeros(size(phi));
%phi=0:pi/200:pi/3
tic
parfor jj=1:length(phi)
Es(jj)=abs(Escat(rho,phi(jj)));
phi_d(jj)=phi(jj).*(180/pi);
end

plot(phi_d,Es,’b–‘,’LineWidth’,2) mom, matrix, scattering MATLAB Answers — New Questions

​

I do most of my work in Linux but sometime have to use Windows for certain applications. I have MATLAB on both but have only recently wanted to use MATLAB on Windows. What is the simplest method for transferring my settings from Linux to Windows? I have my library of m files stored on Dropbox but I do not actually put every subfolder into the path. I also store my `.matlab/` directory in Dropbox (sync setting between my Linux desktop and laptop) and can see my `EmacsModifiedDefaultSet.xml` file with my keyboard preferences but MATLAB says it cannot import that.

I doubt it would be possible to sync the keyboard shortcuts file but if not, please let me know.I do most of my work in Linux but sometime have to use Windows for certain applications. I have MATLAB on both but have only recently wanted to use MATLAB on Windows. What is the simplest method for transferring my settings from Linux to Windows? I have my library of m files stored on Dropbox but I do not actually put every subfolder into the path. I also store my `.matlab/` directory in Dropbox (sync setting between my Linux desktop and laptop) and can see my `EmacsModifiedDefaultSet.xml` file with my keyboard preferences but MATLAB says it cannot import that.

I doubt it would be possible to sync the keyboard shortcuts file but if not, please let me know. I do most of my work in Linux but sometime have to use Windows for certain applications. I have MATLAB on both but have only recently wanted to use MATLAB on Windows. What is the simplest method for transferring my settings from Linux to Windows? I have my library of m files stored on Dropbox but I do not actually put every subfolder into the path. I also store my `.matlab/` directory in Dropbox (sync setting between my Linux desktop and laptop) and can see my `EmacsModifiedDefaultSet.xml` file with my keyboard preferences but MATLAB says it cannot import that.

I doubt it would be possible to sync the keyboard shortcuts file but if not, please let me know. settings, path, windows, linux, preferences MATLAB Answers — New Questions

​

Hi,
i would like to know how using m-script we can create the S-function for a top level subsystem.Manually i am able to achive this by right clicking on the subsystem->C/C++ Code->Generate S-function.Any answers related to S-function viw m-script are welcome
Thanks for your support in advance.Hi,
i would like to know how using m-script we can create the S-function for a top level subsystem.Manually i am able to achive this by right clicking on the subsystem->C/C++ Code->Generate S-function.Any answers related to S-function viw m-script are welcome
Thanks for your support in advance. Hi,
i would like to know how using m-script we can create the S-function for a top level subsystem.Manually i am able to achive this by right clicking on the subsystem->C/C++ Code->Generate S-function.Any answers related to S-function viw m-script are welcome
Thanks for your support in advance. s-function MATLAB Answers — New Questions

​

I am trying to have a play around with simulink and the UAVToolbox but am having issues in simply connecting to the pixhawk.
During the hardware setup for PX4 in simulink, (building firmware, uploading, testing connection etc), I am unable to read accelerometer data, however the odd thing is that I AM able to build and upload the firmware. I am very confused as to why it can connect when uploading, but not when trying to read data. Additionally I dont have this problem when using things such as QGC.
As far as I am aware, the port for uploading is the same as reading (/dev/ttyACM0), regardless of this I have tried various ports and still no connection.
What do I need to do for this to connect and read from my pixhawk?
Note: I am trying to do this with pixhawk 6x as the target to be deployed in to a physical aircraft, NOT SITL or HITL, alot of the examples and tutorials cover only simulated scenarios and they are not helpful in my current situation, so please do not refer me to any of them unless it specifically refers to a connection issue like this, thank you.
System – Ubuntu 20.04
any help would be great,
thanksI am trying to have a play around with simulink and the UAVToolbox but am having issues in simply connecting to the pixhawk.
During the hardware setup for PX4 in simulink, (building firmware, uploading, testing connection etc), I am unable to read accelerometer data, however the odd thing is that I AM able to build and upload the firmware. I am very confused as to why it can connect when uploading, but not when trying to read data. Additionally I dont have this problem when using things such as QGC.
As far as I am aware, the port for uploading is the same as reading (/dev/ttyACM0), regardless of this I have tried various ports and still no connection.
What do I need to do for this to connect and read from my pixhawk?
Note: I am trying to do this with pixhawk 6x as the target to be deployed in to a physical aircraft, NOT SITL or HITL, alot of the examples and tutorials cover only simulated scenarios and they are not helpful in my current situation, so please do not refer me to any of them unless it specifically refers to a connection issue like this, thank you.
System – Ubuntu 20.04
any help would be great,
thanks I am trying to have a play around with simulink and the UAVToolbox but am having issues in simply connecting to the pixhawk.
During the hardware setup for PX4 in simulink, (building firmware, uploading, testing connection etc), I am unable to read accelerometer data, however the odd thing is that I AM able to build and upload the firmware. I am very confused as to why it can connect when uploading, but not when trying to read data. Additionally I dont have this problem when using things such as QGC.
As far as I am aware, the port for uploading is the same as reading (/dev/ttyACM0), regardless of this I have tried various ports and still no connection.
What do I need to do for this to connect and read from my pixhawk?
Note: I am trying to do this with pixhawk 6x as the target to be deployed in to a physical aircraft, NOT SITL or HITL, alot of the examples and tutorials cover only simulated scenarios and they are not helpful in my current situation, so please do not refer me to any of them unless it specifically refers to a connection issue like this, thank you.
System – Ubuntu 20.04
any help would be great,
thanks pixhawk, px4, linux MATLAB Answers — New Questions

​

I’m using N1 210 rev 4 usrp radio and I want to use it in Matlab environment.

My problem is the following. If I type "findsdru" in Matlab command line, then I get the following result.

———- begin libuhd error message output ———-

Device discovery error: input stream error

———- end libuhd error message output ———-

ans =

Platform: ‘N200/N210/USRP2’
IPAddress: ‘192.168.10.2’
SerialNum: ‘E7R26M7UP’
Status: ‘Success’

It seems the Matlab finds the usrp device well, but I couldn’t figure out why the error message "Device discovery error: input stream error" occurs.

How can I fix this?I’m using N1 210 rev 4 usrp radio and I want to use it in Matlab environment.

My problem is the following. If I type "findsdru" in Matlab command line, then I get the following result.

———- begin libuhd error message output ———-

Device discovery error: input stream error

———- end libuhd error message output ———-

ans =

Platform: ‘N200/N210/USRP2’
IPAddress: ‘192.168.10.2’
SerialNum: ‘E7R26M7UP’
Status: ‘Success’

It seems the Matlab finds the usrp device well, but I couldn’t figure out why the error message "Device discovery error: input stream error" occurs.

How can I fix this? I’m using N1 210 rev 4 usrp radio and I want to use it in Matlab environment.

My problem is the following. If I type "findsdru" in Matlab command line, then I get the following result.

———- begin libuhd error message output ———-

Device discovery error: input stream error

———- end libuhd error message output ———-

ans =

Platform: ‘N200/N210/USRP2’
IPAddress: ‘192.168.10.2’
SerialNum: ‘E7R26M7UP’
Status: ‘Success’

It seems the Matlab finds the usrp device well, but I couldn’t figure out why the error message "Device discovery error: input stream error" occurs.

How can I fix this? usrp, error MATLAB Answers — New Questions

​

I am trying to implement the "Highway lane change" example with my own controller. Initially I designed my controller such that if I can provide the controller a reference path ( in the form of a N*2 matrix) the controller will be able to follow that path. Now I want to extend my project such that the reference path will be generated autonomously with the help of sensor fusion. I can see the similar thing has been done in the "Highway lane change" project. But for my controller I need to extract the reference X and Y values. But I am struggling to find those 2 values. I have tries with a matlab function block with the following code –
function [X_ref, Y_ref] = extractWaypoints(TrajectoryInfo)
% Extract X and Y waypoints from the planned trajectory
optimalTrajIndex = TrajectoryInfo.OptimalTrajectoryIndex;
trajectory = TrajectoryInfo.GlobalTrajectory(optimalTrajIndex).Trajectory;
% Extract X and Y positions
X_ref = trajectory(:, 1);
Y_ref = trajectory(:, 2);
end
But with this I am getting the error, "Error:’Input Port 1′ of ‘MotionPlanner/Ref_waypoints’ expects a bus but receives a nonbus signal from ‘Output Port 1’ of ‘MotionPlanner/MATLAB Function2’."
I am not sure how to work on this example.I am trying to implement the "Highway lane change" example with my own controller. Initially I designed my controller such that if I can provide the controller a reference path ( in the form of a N*2 matrix) the controller will be able to follow that path. Now I want to extend my project such that the reference path will be generated autonomously with the help of sensor fusion. I can see the similar thing has been done in the "Highway lane change" project. But for my controller I need to extract the reference X and Y values. But I am struggling to find those 2 values. I have tries with a matlab function block with the following code –
function [X_ref, Y_ref] = extractWaypoints(TrajectoryInfo)
% Extract X and Y waypoints from the planned trajectory
optimalTrajIndex = TrajectoryInfo.OptimalTrajectoryIndex;
trajectory = TrajectoryInfo.GlobalTrajectory(optimalTrajIndex).Trajectory;
% Extract X and Y positions
X_ref = trajectory(:, 1);
Y_ref = trajectory(:, 2);
end
But with this I am getting the error, "Error:’Input Port 1′ of ‘MotionPlanner/Ref_waypoints’ expects a bus but receives a nonbus signal from ‘Output Port 1’ of ‘MotionPlanner/MATLAB Function2’."
I am not sure how to work on this example. I am trying to implement the "Highway lane change" example with my own controller. Initially I designed my controller such that if I can provide the controller a reference path ( in the form of a N*2 matrix) the controller will be able to follow that path. Now I want to extend my project such that the reference path will be generated autonomously with the help of sensor fusion. I can see the similar thing has been done in the "Highway lane change" project. But for my controller I need to extract the reference X and Y values. But I am struggling to find those 2 values. I have tries with a matlab function block with the following code –
function [X_ref, Y_ref] = extractWaypoints(TrajectoryInfo)
% Extract X and Y waypoints from the planned trajectory
optimalTrajIndex = TrajectoryInfo.OptimalTrajectoryIndex;
trajectory = TrajectoryInfo.GlobalTrajectory(optimalTrajIndex).Trajectory;
% Extract X and Y positions
X_ref = trajectory(:, 1);
Y_ref = trajectory(:, 2);
end
But with this I am getting the error, "Error:’Input Port 1′ of ‘MotionPlanner/Ref_waypoints’ expects a bus but receives a nonbus signal from ‘Output Port 1’ of ‘MotionPlanner/MATLAB Function2’."
I am not sure how to work on this example. highway lane change controller modification MATLAB Answers — New Questions

​