skubiai padėti reikalingos šio Matlab inetgration!

[Ana-Maria]

Hi all,
Man reikia apskaičiuoti šio integralas su MATLAB<img src="http://images.elektroda.net/43_1244304476_thumb.jpg" border="0" alt=""/> ir kodas aš nepadeda.Prašome gali kas nors padėti man ir pasakyti, kur yra klaida?

NES = 200;
w = 2 * pi * 10 ^ 9;
M = 1: NES
olätz (m) =- Pi
2 * (m-1) * Pi / (NES-1);
pabaigos;
grondre = 1 / (2 * pi) * trapz (teta, stačiakampius (teta, m) .* cos (k * olätz));
grondim = 1 / (2 * pi) * trapz (teta, stačiakampius (teta, m) .* sin (k * olätz));
grondx = grondre i * grondim;

Ačiū!

 

[AndyECE]

Jūs trūkstamą informaciją:

kas yra "stačiakampės" funkcija?Jeigu jūs darysite langas, jūs negalėsite atlikti jį "w" ilgas (m yra 1e9 ilgio!)

Kaip vadinama funkcija esate integruoti?Gal galėtumėte aiškiai parašyti, kas G (teta; f (w)) yra?

ačiū

 

[Ana-Maria]

Labas,
G (teta, m) = stačiakampio formos, yra pirmoji dalis, integravimo (antroji dalis exp (-j * KB * eta))
ir kodas G (teta, m)

funkcija [G] = stačiakampius (teta, m);
c = 3 * 10 ^ 8;
G0 = 1;
f = 2 * 10 ^ 9;
h = [c / (4 * f) c / (2 * f) 3 * c / f, c / f];
Ordin = 2;
x = m * h (Ordin) / (2 * c) * sin (teta);
ir k = 1:
length (x)
jei x (k) == 0
G (k) = G0;
kitas
G (k) = G0 * sin (x (k)) / x (k);
pabaiga
pabaigašiuo metu f (w) = m ir m = 2 * 10 ^ 9