• Remix
  • Share
  • New Entry

on 21 Oct 2021
  • 2
  • 28
  • 1
  • 0
  • 278
clear;
close all;
X=3e2;
x=(0:X-1)/X;
u=rand(size(x));
b=[0:X/2-1,0,-X/2+1:-1];
D=2i*pi*b/X;
l=-D.^2-D.^4;
I=ones(1,X);
A=I+l/20;
C=-D.*fft(u.*u)/2;
t=(0:999)*.1;
U(1,:)=u;
u=fft(u);
for n=1:1e3
P=C;
C=-D.*fft(real(ifft(u)).^2)/2;
u=(A.*u+.15*C-P/20).*(I-l/20).^(-1);
U(n,:)=fft(u);
end
pcolor(real(U));
shading interp
Remix Tree