n = 500;
t = linspace(0, 10*pi, n);
fprintf("Time domain:\nex=%.1f, ey=%.1f\n", ex
, ey
) ex
=sum
(abs(fft
(x
)).
^2)/n
;ey
=sum
(abs(fft
(y
)).
^2)/n
; fprintf("Frequency domain (rectwin):\nex=%.1f, ey=%.1f\n", ex
, ey
) w = hamming(n)';
ex=sum(abs(fft(x.*w)).^2)/sum(w.^2);
ey=sum(abs(fft(y.*w)).^2)/sum(w.^2);
fprintf("Frequency domain (hamming):\nex=%.1f, ey=%.1f\n", ex, ey)
biA9IDUwMDsKdCA9IGxpbnNwYWNlKDAsIDEwKnBpLCBuKTsKeCA9IHNpbigyKnBpKnQpOwp5ID0gc2luKDIqcGkqdCArIHBpLzIpOwoKZXg9c3VtKGFicyh4KS5eMik7ICAgIApleT1zdW0oYWJzKHkpLl4yKTsKZnByaW50ZigiVGltZSBkb21haW46XG5leD0lLjFmLCBleT0lLjFmXG4iLCBleCwgZXkpCmV4PXN1bShhYnMoZmZ0KHgpKS5eMikvbjsKZXk9c3VtKGFicyhmZnQoeSkpLl4yKS9uOyAgCmZwcmludGYoIkZyZXF1ZW5jeSBkb21haW4gKHJlY3R3aW4pOlxuZXg9JS4xZiwgZXk9JS4xZlxuIiwgZXgsIGV5KQp3ID0gaGFtbWluZyhuKSc7CmV4PXN1bShhYnMoZmZ0KHguKncpKS5eMikvc3VtKHcuXjIpOwpleT1zdW0oYWJzKGZmdCh5Lip3KSkuXjIpL3N1bSh3Ll4yKTsKZnByaW50ZigiRnJlcXVlbmN5IGRvbWFpbiAoaGFtbWluZyk6XG5leD0lLjFmLCBleT0lLjFmXG4iLCBleCwgZXkp