fazės masyvas beaforming naudojant MATLAB

S

sasamdyu

Guest
I need help on my kodavimo nes i cant figure out dalį, kuri yra klaidinga, nes man gali atrodyti, kad gauti išvesties Noriu ...

funkcijos [] = matricas;
aiški visiems;
beveik visi;
CLC;
CLF;

k = 2 * pi;
Mx = 5;
Ny = 5;
DX = 0,25;
dy = 0,25;
thmax2 = 0;
phimax2 = 0;
dtor = pi/180;
teta = linspace (1.180.180);
Phi = linspace (1.180.180);

betax =- k * Dx * SIN (dtor * thmax2) * cos (dtor * phimax2);
betay =- k * dy * sin (dtor * thmax2) * sin (dtor * phimax2);

[Theta, Phi] = meshgrid (teta, phi);
AF3 = afnf (teta, Phi, Mx, NY, dx, dy, betax, betay);
spherical_plot (AF3, teta, PHI)
ß = pavadinimas ('3 D Sferinės sklypas masyvo Factor "," Šrifto dydis ", 15);
pav
rect_plot (AF3)

% rutulio sklypas
funkcija spherical_plot (R, teta, PHI)
dtor = pi/180;
x = abs (R) .* sin (teta * dtor) .* cos (PHI * dtor);
y = abs (R) .* sin (teta * dtor) .* sin (PHI * dtor);
z = abs (R) .* cos (Theta * dtor);
surf (x, y, z);
Žiūrėti (0.180);
C = [.8, .8, ,8]; paletės (C); ašies ne vienodas;

funkcija rect_plot (R)
S = abs (R);
% Akių (-1 * (Length (s) / 2-1): 1: length (s) / 2, -1 * (Length (s) / 2-1): 1: length (s) / 2, S / norma (S))
% Akių (-179:1:179, -45:1:45, S / norma (S))
tinkleliai (-89:1:90, -89:1:90, S / norma (S));

% akių (-45:1:45, -179:1:179, S / norma (S));

Pavadinimas ( "Stačiakampiai Sklypas Planar Array spinduliavimo dėsningumo)
xlabel (kampas teta-xaxis indeksuotų m laipsniais)
ylabel (kampas Phi-yaxis indeksuoti N laipsniais)
zlabel (Normalizuotas Amlitude)

funkcija [F10] = afnf (theta1, phi1, Mx, NY, dx, dy, betax, betay)
k = 2 * pi;
dtor = pi/180;
AFx1 (MX) = 30;
AFy1 (NY) = 30;
už teta = 1:1:180;
už Phi = 1:1:180;
elementx = 1;
Elementai = 1;
AFx1 (elementx) = 0;
AFy1 (elementai) = 0;
AXF = 0;
psix = k * DX * sin (teta * dtor-pi / 2) * cos (phi * dtor-pi / 2) betax;
psiy = k * dy * sin (teta * dtor-pi / 4) * sin (phi * dtor-pi / 2) betay;
už elementx = 1: MX
AFx2 (elementx) = exp (I * elementx * psix);
AFx3 (elementx) = AFx2 (elementx) AFx1 (elementx);
AFx1 (elementx 1) = AFx3 (elementx);
pabaigoje

AXF = AFx1 (MX);
už elementy = 1: Naujas;
AFy2 (elementai) = exp (I * Elementai * psiy) * AXF;
AFy3 (elementai) = AFy2 (elementai) AFy1 (elementai);
AFy1 (elementy 1) = AFy3 (elementai);
pabaigoje
AYF = AFy1 (NY);
AF (PHI, teta) = AYF;
pabaigoje
pabaigoje

F10 = AF;

 
kurios produkcija tu nori?

Aš išbandžiau savo kodą ir atrodo, kad jis veikia gerai ..
Atsiprašome, bet jums reikia prisijungęs, kad galėtumėte peržiūrėti šį priedą

 
angolaX rašė:

kurios produkcija tu nori?Aš išbandžiau savo kodą ir atrodo, kad jis veikia gerai ..
 
Hi vyras, aš neturiu laiko patikrinti savo programą, tačiau čia atteched yra kodas, kad aš naudoju iki sklypo 3D stačiakampio masyvo.
Jūs galite tiesiog ištrinti weightning W gauti Uniforme matrica aš tiesiog HPE, kad jis jums padėjo.

Kodas

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

%%%%%%%%%%% Rayonnement 2D par Dolph-Chebyshev %%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

% perdirbimo ir stačiakampio masyvas Nx * Ny antenos, su masyvo

% atitinkamai dx ir dy tarpais išilgai X ir Y ašis su RDB SLR

% su Dolph chebychev wieghtningfunkcija Cheb_2D (nx, ny, dx, dy, RDB)dir_theta = 20;

dir_phi = 30;dir_ux = sin ((dir_theta/180) * pi) * cos ((dir_phi/180) * pi);

dir_uy = sin ((dir_theta/180) * pi) * SIN ((dir_phi/180) * pi);

%%%%%% Calcul des coefs de Cheb korespondentų au reseau%%%%

w = planar_dolphcheby (dx, dy, nx, ny, RDB);%%% panaudojimas de ces coefs pour le calcul du Facteur reseau %%%%%%%%%%%%%%%%%% CAS paire %%%%%%%%%%%%%%%%%%%%%%%

keistai = rem (Nx, 2);

if (nelyginis)

hh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

P = fix (Nx / 2);

, n = M: M

M = M: M

B (x, y) = B (x, y) W (N M 1, M M 1) * exp (-j * 2 * pi * (n * UX * dx M * UY * dy) ) * exp (-j * 2 * pi * (n * dir_ux * dx M * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)kitashh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

M = 1: Nx

, n = 1: Ny

B (x, y) = B (x, y) W (m, n) * exp (-j * 2 * pi * (m * UX * dx n * UY * dy)) * exp (-j * 2 * pi * (m * dir_ux * dx n * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)pabaigaB = abs (b) / max (max (ABS (B)));

dydis (B);

tinklelis

surf (ux_grid, uy_grid, (20 * log 10 (B ))');

axis ([-1 1 -1 1 -70 0])

 
Albert wrote:

Hi vyras, aš neturiu laiko patikrinti savo programą, tačiau čia atteched yra kodas, kad aš naudoju iki sklypo 3D stačiakampio masyvo.

Jūs galite tiesiog ištrinti weightning W gauti Uniforme matrica aš tiesiog HPE, kad jis jums padėjo.Kodas

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

%%%%%%%%%%% Rayonnement 2D par Dolph-Chebyshev %%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

% perdirbimo ir stačiakampio masyvas Nx * Ny antenos, su masyvo

% atitinkamai dx ir dy tarpais išilgai X ir Y ašis su RDB SLR

% su Dolph chebychev wieghtningfunkcija Cheb_2D (nx, ny, dx, dy, RDB)dir_theta = 20;

dir_phi = 30;dir_ux = sin ((dir_theta/180) * pi) * cos ((dir_phi/180) * pi);

dir_uy = sin ((dir_theta/180) * pi) * SIN ((dir_phi/180) * pi);

%%%%%% Calcul des coefs de Cheb korespondentų au reseau%%%%

w = planar_dolphcheby (dx, dy, nx, ny, RDB);%%% panaudojimas de ces coefs pour le calcul du Facteur reseau %%%%%%%%%%%%%%%%%% CAS paire %%%%%%%%%%%%%%%%%%%%%%%

keistai = rem (Nx, 2);

if (nelyginis)

hh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

M = fix (Nx / 2);

, n = M: M

M = M: M

B (x, y) = B (x, y) W (N M 1, M M 1) * exp (-j * 2 * pi * (n * UX * dx M * UY * dy) ) * exp (-j * 2 * pi * (n * dir_ux * dx M * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)kitashh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

M = 1: Nx

, n = 1: Ny

B (x, y) = B (x, y) W (m, n) * exp (-j * 2 * pi * (m * UX * dx n * UY * dy)) * exp (-j * 2 * pi * (m * dir_ux * dx n * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)pabaigaB = abs (b) / max (max (ABS (B)));

dydis (B);

tinklelis

surf (ux_grid, uy_grid, (20 * log 10 (B ))');

axis ([-1 1 -1 1 -70 0])

 
Beje, tai klaida, kad sasamdyu pasireiškia:

"Function apibrėžimai neleidžiama eilutę arba scenarijus"

ieškoti google atrodo, kad funkcijos negali būti įdedama viduje scenarijus.kas nors žino kaip išspręsti šią klaidą?

 
Albert wrote:

Hi vyras, aš neturiu laiko patikrinti savo programą, tačiau čia atteched yra kodas, kad aš naudoju iki sklypo 3D stačiakampio masyvo.

Jūs galite tiesiog ištrinti weightning W gauti Uniforme matrica aš tiesiog HPE, kad jis jums padėjo.Kodas

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

%%%%%%%%%%% Rayonnement 2D par Dolph-Chebyshev %%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%

% perdirbimo ir stačiakampio masyvas Nx * Ny antenos, su masyvo

% atitinkamai dx ir dy tarpais išilgai X ir Y ašis su RDB SLR

% su Dolph chebychev wieghtningfunkcija Cheb_2D (nx, ny, dx, dy, RDB)dir_theta = 20;

dir_phi = 30;dir_ux = sin ((dir_theta/180) * pi) * cos ((dir_phi/180) * pi);

dir_uy = sin ((dir_theta/180) * pi) * SIN ((dir_phi/180) * pi);

%%%%%% Calcul des coefs de Cheb korespondentų au reseau%%%%

w = planar_dolphcheby (dx, dy, nx, ny, RDB);%%% panaudojimas de ces coefs pour le calcul du Facteur reseau %%%%%%%%%%%%%%%%%% CAS paire %%%%%%%%%%%%%%%%%%%%%%%

keistai = rem (Nx, 2);

if (nelyginis)

hh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

M = fix (Nx / 2);

, n = M: M

M = M: M

B (x, y) = B (x, y) W (N M 1, M M 1) * exp (-j * 2 * pi * (n * UX * dx M * UY * dy) ) * exp (-j * 2 * pi * (n * dir_ux * dx M * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)kitashh = waitbar (0, "Perdirbimas ...');

ux_grid = -1:0.01:1;

uy_grid = -1:0.01:1;

kai x = 1: length (ux_grid)

waitbar (x / ilgis (ux_grid));

UX = ux_grid (x);

už y = 1: length (uy_grid)uy = uy_grid (y);

B (x, y) = 0;

M = 1: Nx

, n = 1: Ny

B (x, y) = B (x, y) W (m, n) * exp (-j * 2 * pi * (m * UX * dx n * UY * dy)) * exp (-j * 2 * pi * (m * dir_ux * dx n * dir_uy * dy));

pabaiga

pabaiga

pabaiga

pabaiga

arti (HH)pabaigaB = abs (b) / max (max (ABS (B)));

dydis (B);

tinklelis

surf (ux_grid, uy_grid, (20 * log 10 (B ))');

axis ([-1 1 -1 1 -70 0])

 
hello ..nors aš vis dar cant paleisti kodavimo ir man ...Help ...Skelbimas po 6 minučių:Hello Albert ...U davė man šį kodavimo bet kažkaip i cant paleisti programą ..Help ...Kodas

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%% Rayonnement 2D par Dolph-Chebyshev %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%% perdirbimo ir stačiakampio masyvas Nx * Ny antenos, su masyvo% atitinkamai dx ir dy tarpais išilgai X ir Y ašis su RDB SLR% su Dolph chebychev wieghtning
funkcija Cheb_2D (nx, ny, dx, dy, RDB)
dir_theta = 20;dir_phi = 30;
dir_ux = sin ((dir_theta/180) * pi) * cos ((dir_phi/180) * pi);dir_uy = sin ((dir_theta/180) * pi) * SIN ((dir_phi/180) * pi);%%%%%% Calcul des coefs de Cheb korespondentų au reseau%%%%w = planar_dolphcheby (dx, dy, nx, ny, RDB);
%%% panaudojimas de ces coefs pour le calcul du Facteur reseau %%%%%%%%%
%%%%%%%%% CAS paire %%%%%%%%%%%%%%%%%%%%%%%keistai = rem (Nx, 2);if (nelyginis)hh = waitbar (0, \ "Processing ... \");ux_grid = -1:0.01:1;uy_grid = -1:0.01:1;kai x = 1: length (ux_grid)waitbar (x / ilgis (ux_grid));UX = ux_grid (x);už y = 1: length (uy_grid)
uy = uy_grid (y);B (x, y) = 0;M = fix (Nx / 2);, n = M: MM = M: MB (x, y) = B (x, y) W (N M 1, M M 1) * exp (-j * 2 * pi * (n * UX * dx M * UY * dy) ) * exp (-j * 2 * pi * (n * dir_ux * dx M * dir_uy * dy));pabaigapabaigapabaigapabaigaarti (HH)kitas
hh = waitbar (0, \ "Processing ... \");ux_grid = -1:0.01:1;uy_grid = -1:0.01:1;kai x = 1: length (ux_grid)waitbar (x / ilgis (ux_grid));UX = ux_grid (x);už y = 1: length (uy_grid)
uy = uy_grid (y);B (x, y) = 0;M = 1: Nx, n = 1: NyB (x, y) = B (x, y) W (m, n) * exp (-j * 2 * pi * (m * UX * dx n * UY * dy)) * exp (-j * 2 * pi * (m * dir_ux * dx n * dir_uy * dy));pabaigapabaigapabaigapabaigaarti (HH)
pabaiga
B = abs (b) / max (max (ABS (B)));dydis (B);tinklelissurf (ux_grid, uy_grid, (20 * log 10 (B)) \ ');axis ([-1 1 -1 1 -70 0])
 

Welcome to EDABoard.com

Sponsor

Back
Top