这是我要运行的代码。
tmax = 10;
dt = 0.01;
N = tmax/dt;
t = 0:dt:(tmax-dt);
u1 = sin(2*pi*15*t) + sin(2*pi*40*t);
U = (1/N)*fft(u1);
amp = abs(U);
pha = unwrap(angle(U)) * 180/pi ;
df = 1/tmax;
fmax = 1/dt-df;
f = 0:df:fmax;
u2 = N*ifft(U);
resp = [t; u1; real(u2); f; real(U); imag(U); amp; pha]';
save resp.dat resp -ascii
这是我遇到的错误。
使用保存错误
无法写入文件resp.dat:权限被拒绝。
ch01_fft错误(第27行)
保存resp.dat resp -ascii