[MATLAB] complex permittivity for earth ionosphere syntax

[s_hy]

Homework Statement

I have these equations:

ga = exp(-$σ^{R}$Δt/$ε_{0}$*$ε^{R}_{r}$)

Homework Equations

$σ^{R}$=($ε_{0}$*$ω^{2}_{p}$*v)/($ω^{2}$+$v^{2}$)

The Attempt at a Solution

%ionospheric profile for earth-ionosphere waveguide as coefficients
e = 1.602e-19;          % electron charge
ms = 9.109e-31;         % electron mass

for h = 10:150^3;
    v = 1.816*10^11*exp(-0.15*h);
end

load N191204.mat N
A = N(1,:);                   % electron density, N
omegap = A*e.^2/ms*eps;
epsr   = 1-(omegap.^2/omega.^2+v.^2);
sigmar = (eps*omegap.^2*v)/(omega.^2+v.^2);

%property coefficient
for i=1:k                                  
    for j=1:k;
        ga(i,j)=exp(-(sigmar*dt)/(eps*epsr));  
        gb(i,j)=(1-ga(i,j))/sigmar;              
    end;
end;

the problem is, i have problem with matrix dimension.

Warning: Rank deficient, rank = 0, tol = NaN.
> In propagation2d at 60
Error using /
Matrix dimensions must agree.

Error in propagation2d (line 61)
gb(i,j)=(1-ga(i,j))/sigmar;

i declare ga(i,j) in 2d because ga value is depend on height and distance as well but sigmar is in 1d (the value of sigma is depends on v where v depends on h-height only)..

anyone have idea what i should do in order to use the value of sigmar in gb(i,j)..i use MATLAB for this coding.

thank you.

 

[KataKoniK]

it is important to carefully analyze and understand the equations and variables involved in the problem. In this case, it seems like there may be a mismatch in dimensions between the variables used in the equations and those used in the code.

The first step would be to check the dimensions of each variable and ensure that they align with the equations. It may be helpful to write out the equations in a more simplified form to better understand the relationships between the variables.

Once the dimensions have been confirmed, the next step would be to check the logic and syntax of the code to ensure that it accurately reflects the equations. It may also be helpful to break down the code into smaller sections and test them individually to identify where the issue may be occurring.

If the problem persists, it may be beneficial to seek assistance from a colleague or consult other resources, such as textbooks or online forums, for further guidance. It is important to fully understand the problem and any potential errors in order to find an appropriate solution.