效果图
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(___,