% 北太天元畫圖（紅心、金龍)

?clc

?clf

?close all

clear all

load_plugin("time");

fill([-12 -12 8 8],[-12 4 4 -12],[?0.5843???0.8157???0.9882]);

%axis equal;

%axis([-12 7 -12 3]);

hold on;

text(-10,2,'謝謝你', 'Color','red','FontSize',30)

%畫一個紅心

t=linspace(-2.9,2.9,1000);

X=16.*(sin(t)).^3;

Y=13.*cos(t)-5.*cos(2.*t)-2.*cos(3.*t)-cos(4.*t);

fill(X/10-7,Y/10-3,[180,39,45]./255,'EdgeColor',[180,39,45]./255,'LineWidth',2)

hold on

pause(1)

%畫一條金黃色的巨龍

w=[0.824074 0.281482 -0.212346 0.864198 -1.882290 -0.110607

???0.088272 0.520988 -0.463889 -0.377778 0.785360 8.095795];

N=10000;

x=zeros(1,N);

y=zeros(1,N);

x0=0;

y0=0;

for i=1:N

???if rand<0.8

???????j=1;

???else

???????j=2;

???end

???x01=w(j,[1 2 5])*[x0,y0,1]';

???y01=w(j,[3 4 6])*[x0,y0,1]';

???x(i)=x01;

???y(i)=y01;

???x0=x01;

???y0=y01;

end

%藍色的龍

%plot(x,y,'.', 'MarkerSize', 60, 'Color', [0.2,0.2,0.9])

%金黃色的龍

scatter(x,y-8,200,'o','filled','MarkerFaceColor',[255/255,233/255,0])

hold on

pause(1)

text(-10,-9,'燕園相逢多少秋，今秋又遇論胖瘦。', 'Color','red','FontSize',30)

text(-10,-11,'年年歲歲人不同，慢性鼻炎仍依舊。', 'Color','red','FontSize',30)

hold off

unload_plugin("time");