MATLAB | 迟到的圣诞树绘制

fig=figure('Units','normalized','Position',[.1,.05,.4,.8]); ax=axes(fig); hold(ax,'on'); ax.XTick=[]; ax.YTick=[]; ax.XColor='none'; ax.YColor='none'; ax.XLim=[0,600]; ax.YLim=[0,800]; ax.DataAspectRatio=[1 1 1]; ax.Color=[109,128,205]./255; drawMist() % 绘制雾气 drawMandS() % 绘制星星和月亮 drawMountains() % 绘制山脉和雪地 drawXMasTree() % 绘制圣诞树

% 雾气绘制函数
function drawMist() colMist=[181,137,190]./255; a=500;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 二维柏林式噪声生成
    alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colMist(1); CMesh(:,:,2)=colMist(2); CMesh(:,:,3)=colMist(3); % 越向上、雾气越透明
    dy=(500:-1:1)./500.*.8+.2; image([0,600],[200,800],CMesh,'AlphaData',alp.*(dy')); end

绘制效果：

% 星星月亮绘制函数
function drawMandS() % 绘制双层圆形当月亮
    t=0:pi/500:2*pi; fill(cos(t).*25+450,sin(t).*25+670,[236,235,253]./255,'EdgeColor','none') fill(cos(t).*30+450,sin(t).*30+670,[236,235,253]./255,'EdgeColor','none','FaceAlpha',.3) % 绘制散点当星星 scatter(randi([0,600],[15,1]),randi([200,800],[15,1]),4,'filled',... 'CData',[1 1 1],'MarkerFaceAlpha',.6); end

绘制效果：

% 山峰绘制函数
function drawMountains() % 绘制山峰- - - - - - - - - - - - - - - - - - - - - - - - - - - - - colFurther=[240 34 75]; % 远山的颜色
    colCloser=[219 43 60]; % 近山的颜色

    a=800;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 8层山 for i=1:8 % 每次都生成一次二维柏林式噪声，并取其中一行的数据
        h=abs(ifftn(exp(5i*rand(a))./r.^1.05)).*(c*c').*10; nh=(8-i)*30+h(400,:)+180; hm=ceil(max(nh)); CMesh=zeros([hm,800,3]); % 颜色矩阵构造，
        tcol=colFurther+(colCloser-colFurther)./8.*(i); tcol=hsv2rgb(tcol./[360,100,100]); CMesh(:,:,1)=tcol(1); CMesh(:,:,2)=tcol(2); CMesh(:,:,3)=tcol(3); % 用nan数值框出山的轮廓
        alp=ones(hm,800); alp((1:hm)'>nh)=nan; % 绘制山峰 image([-50,650],[0,hm],CMesh,'AlphaData',alp.*0.98); end % 画个雪地- - - - - - - - - - - - - - - - - - - - - - - - - - - - - nh=h(800,:); fill([linspace(0,600,length(nh)),600,0],[nh./1.5+140,0,0],[218,222,234]./255,... 'EdgeColor','none') % 雪地附近雾气- - - - - - - - - - - - - - - - - - - - - - - - - - - colMist=[155,172,198]./255.*1.2; a=500;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 二维柏林式噪声生成
    alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colMist(1); CMesh(:,:,2)=colMist(2); CMesh(:,:,3)=colMist(3); % 越向上、雾气越透明
    dy=(500:-1:1)./500.*.6+.2; image([0,600],[0,300],CMesh,'AlphaData',alp.*(dy')); end

绘制效果：

% 圣诞树绘制函数
function drawXMasTree() % 绘制树干- - - - - - - - - - - - - - - - - - - - - - - - - - - - - c=zeros(1,4,3); c(1,1:2,1)=163./255;c(1,1:2,2)=165./255;c(1,1:2,3)=183./255; c(1,3:4,1)=85./255 ;c(1,3:4,2)=95./255 ;c(1,3:4,3)=130./255; patch([283,317,317,283],[100,100,180,180],c,'FaceColor','interp','EdgeColor','none') fill([283:317,317,283],[100+rand(1,35).*10,90,90],[207,217,235]./255,'EdgeColor','none') % 绘制树干粗糙花纹- - - - - - - - - - - - - - - - - - - - - - - - - a=400;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); alp=alp(1:100,:); colTrunk=[160,162,177]./255; CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colTrunk(1); CMesh(:,:,2)=colTrunk(2); CMesh(:,:,3)=colTrunk(3); dy=(100:-1:1)./100.*.8+.2; image([283,317],[100,180],CMesh,'AlphaData',alp.*(dy')); % 绘制树叶- - - - - - - - - - - - - - - - - - - - - - - - - - - - - colTree=[64,88,124]./255; a=100;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; CMesh=zeros([a,a,3]); CMesh(:,:,1)=colTree(1); CMesh(:,:,2)=colTree(2); CMesh(:,:,3)=colTree(3); jList=[1,2,3,5,6,8,10,12,14,16,16,18,18,20,20,24,24,26,26,28,... 28,30,30,32,32,32,32,34,32,32,30,26,26,20]; for i=1:34 % 二维柏林式噪声生成
        alp=abs(ifftn(exp(7i*rand(a))./r.^.9)).*(c*c'); alp=alp./max(alp,[],'all'); alp(alp<.5)=0; alp(alp>=.5)=.7; for j=(1:jList(i))-jList(i)/2 image([270,330]+j*6+(rand(1)-.5).*6,... [200,260]+400-i*13+(rand(1)-.5).*3,CMesh,'AlphaData',alp); end
    end % 绘制五角星- - - - - - - - - - - - - - - - - - - - - - - - - - - - t=linspace(pi/2,2*pi+pi/2,11); l=ones(1,11);l(2:2:11)=.47; fill(cos(t).*26.*l+306,sin(t).*26.*l+640,[237,231,183]./255,... 'EdgeColor','none','FaceAlpha',.9) % 绘制其他装饰物- - - - - - - - - - - - - - - - - - - - - - - - - - % 绘制彩带
    siny=-abs(sin(linspace(0,2.3*pi,100)))-linspace(0,2,100); plot(linspace(200,400,100),siny.*20+300,'Color',[198,137,159]./255,'LineWidth',2) siny=-abs(sin([linspace(0,pi,100),linspace(0,pi,40)]))-linspace(0,-5,140); plot(linspace(220,380,140),siny.*20+330,'Color',[198,137,159]./255,'LineWidth',2) siny=-abs(sin(linspace(0,1.2*pi,100)))-linspace(0,1.6,100); plot(linspace(260,350,100),siny.*20+520,'Color',[198,137,159]./255,'LineWidth',2) % 绘制装饰圆球 scatter([-50,20,20,50]+300,[-80,-25,70,-160]+400,50,'filled','CData',[152,173,176]./255); scatter([-45,0,20,-60]+300,[110,30,140,-160]+400,50,'filled','CData',[109,133,212]./255); scatter([-30,70]+300,[90,-180]+400,50,'filled','CData',[203,202,181]./255); end

绘制效果：

function XmasTree fig=figure('Units','normalized','Position',[.1,.05,.4,.8]); ax=axes(fig); hold(ax,'on'); ax.XTick=[]; ax.YTick=[]; ax.XColor='none'; ax.YColor='none'; ax.XLim=[0,600]; ax.YLim=[0,800]; ax.DataAspectRatio=[1 1 1]; ax.Color=[109,128,205]./255; drawMist() % 绘制雾气 drawMandS() % 绘制星星和月亮 drawMountains() % 绘制山脉和雪地 drawXMasTree() % 绘制圣诞树 % ========================================================================= % 雾气绘制函数
    function drawMist() colMist=[181,137,190]./255; a=500;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 二维柏林式噪声生成
        alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colMist(1); CMesh(:,:,2)=colMist(2); CMesh(:,:,3)=colMist(3); % 越向上、雾气越透明
        dy=(500:-1:1)./500.*.8+.2; image([0,600],[200,800],CMesh,'AlphaData',alp.*(dy')); end % ------------------------------------------------------------------------- % 星星月亮绘制函数
    function drawMandS() % 绘制双层圆形当月亮
        t=0:pi/500:2*pi; fill(cos(t).*25+450,sin(t).*25+670,[236,235,253]./255,'EdgeColor','none') fill(cos(t).*30+450,sin(t).*30+670,[236,235,253]./255,'EdgeColor','none','FaceAlpha',.3) % 绘制散点当星星 scatter(randi([0,600],[15,1]),randi([200,800],[15,1]),4,'filled',... 'CData',[1 1 1],'MarkerFaceAlpha',.6); end % ------------------------------------------------------------------------- % 圣诞树绘制函数
    function drawXMasTree() % 绘制树干- - - - - - - - - - - - - - - - - - - - - - - - - - - - - c=zeros(1,4,3); c(1,1:2,1)=163./255;c(1,1:2,2)=165./255;c(1,1:2,3)=183./255; c(1,3:4,1)=85./255 ;c(1,3:4,2)=95./255 ;c(1,3:4,3)=130./255; patch([283,317,317,283],[100,100,180,180],c,'FaceColor','interp','EdgeColor','none') fill([283:317,317,283],[100+rand(1,35).*10,90,90],[207,217,235]./255,'EdgeColor','none') % 绘制树干粗糙花纹- - - - - - - - - - - - - - - - - - - - - - - - - a=400;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); alp=alp(1:100,:); colTrunk=[160,162,177]./255; CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colTrunk(1); CMesh(:,:,2)=colTrunk(2); CMesh(:,:,3)=colTrunk(3); dy=(100:-1:1)./100.*.8+.2; image([283,317],[100,180],CMesh,'AlphaData',alp.*(dy')); % 绘制树叶- - - - - - - - - - - - - - - - - - - - - - - - - - - - - colTree=[64,88,124]./255; a=100;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; CMesh=zeros([a,a,3]); CMesh(:,:,1)=colTree(1); CMesh(:,:,2)=colTree(2); CMesh(:,:,3)=colTree(3); jList=[1,2,3,5,6,8,10,12,14,16,16,18,18,20,20,24,24,26,26,28,... 28,30,30,32,32,32,32,34,32,32,30,26,26,20]; for i=1:34 % 二维柏林式噪声生成
            alp=abs(ifftn(exp(7i*rand(a))./r.^.9)).*(c*c'); alp=alp./max(alp,[],'all'); alp(alp<.5)=0; alp(alp>=.5)=.7; for j=(1:jList(i))-jList(i)/2 image([270,330]+j*6+(rand(1)-.5).*6,... [200,260]+400-i*13+(rand(1)-.5).*3,CMesh,'AlphaData',alp); end
        end % 绘制五角星- - - - - - - - - - - - - - - - - - - - - - - - - - - - t=linspace(pi/2,2*pi+pi/2,11); l=ones(1,11);l(2:2:11)=.47; fill(cos(t).*26.*l+306,sin(t).*26.*l+640,[237,231,183]./255,... 'EdgeColor','none','FaceAlpha',.9) % 绘制其他装饰物- - - - - - - - - - - - - - - - - - - - - - - - - - % 绘制彩带
        siny=-abs(sin(linspace(0,2.3*pi,100)))-linspace(0,2,100); plot(linspace(200,400,100),siny.*20+300,'Color',[198,137,159]./255,'LineWidth',2) siny=-abs(sin([linspace(0,pi,100),linspace(0,pi,40)]))-linspace(0,-5,140); plot(linspace(220,380,140),siny.*20+330,'Color',[198,137,159]./255,'LineWidth',2) siny=-abs(sin(linspace(0,1.2*pi,100)))-linspace(0,1.6,100); plot(linspace(260,350,100),siny.*20+520,'Color',[198,137,159]./255,'LineWidth',2) % 绘制装饰圆球 scatter([-50,20,20,50]+300,[-80,-25,70,-160]+400,50,'filled','CData',[152,173,176]./255); scatter([-45,0,20,-60]+300,[110,30,140,-160]+400,50,'filled','CData',[109,133,212]./255); scatter([-30,70]+300,[90,-180]+400,50,'filled','CData',[203,202,181]./255); end % ------------------------------------------------------------------------- % 山峰绘制函数
    function drawMountains() % 绘制山峰- - - - - - - - - - - - - - - - - - - - - - - - - - - - - colFurther=[240 34 75]; % 远山的颜色
        colCloser=[219 43 60]; % 近山的颜色
        
        a=800;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 8层山 for i=1:8 % 每次都生成一次二维柏林式噪声，并取其中一行的数据
            h=abs(ifftn(exp(5i*rand(a))./r.^1.05)).*(c*c').*10; nh=(8-i)*30+h(400,:)+180; hm=ceil(max(nh)); CMesh=zeros([hm,800,3]); % 颜色矩阵构造，
            tcol=colFurther+(colCloser-colFurther)./8.*(i); tcol=hsv2rgb(tcol./[360,100,100]); CMesh(:,:,1)=tcol(1); CMesh(:,:,2)=tcol(2); CMesh(:,:,3)=tcol(3); % 用nan数值框出山的轮廓
            alp=ones(hm,800); alp((1:hm)'>nh)=nan; % 绘制山峰 image([-50,650],[0,hm],CMesh,'AlphaData',alp.*0.98); end % 画个雪地- - - - - - - - - - - - - - - - - - - - - - - - - - - - - nh=h(800,:); fill([linspace(0,600,length(nh)),600,0],[nh./1.5+140,0,0],[218,222,234]./255,... 'EdgeColor','none') % 雪地附近雾气- - - - - - - - - - - - - - - - - - - - - - - - - - - colMist=[155,172,198]./255.*1.2; a=500;b=(.5:a)'/a; c=(-cos(b*2*pi)+1).^.2; f=b-.5;r=f'.^2+f.^2; % 二维柏林式噪声生成
        alp=abs(ifftn(exp(3i*rand(a))./r.^.8)).*(c*c'); alp=alp./max(alp,[],'all'); CMesh=zeros([size(alp),3]); CMesh(:,:,1)=colMist(1); CMesh(:,:,2)=colMist(2); CMesh(:,:,3)=colMist(3); % 越向上、雾气越透明
        dy=(500:-1:1)./500.*.6+.2; image([0,600],[0,300],CMesh,'AlphaData',alp.*(dy')); end
end