效果图

MATLAB源码

clc;clear;close all;f=1000;k=255;colormap([ones(22,3);[ones(1,k);linspace(.9,.6,k);linspace(.4,0,k)]']);X=linspace(-164,164,f);[t,r]=cart2pol(X,X');h=r;h(r<10)=1;s=rescale(abs(ifft2(h.^-2.*cos(7*rand(f)))),15,270);s(r<11)=r(r<11);t=sin(t*10);s(t>=-1&t<-.95|r>90)=1;s(r>95&r<100)=80;image(rescale(200*t+80*s,1,k))camva(3.8)

知识回顾

rescale

数组元素的缩放范围

语法

B = rescale(A)

B = rescale(A,l,u)

B = rescale(___,