2Â¥: Originally posted by sogooda at 2012-04-06 09:31:01:
¾ÝÎÒËù֪ûÓÐÕâÖÖÏֳɵÄÃüÁһ¸öÃüÁî¾ÍÄÜʵÏÖÂ¥Ö÷Ëù˵µÄ¹¦ÄÜ£©¡£
µ«ÊÇÓÃmatlab×Ô¼º±à³ÌµÄ»°
¡°ÔÚ»­³öÖʵãÔ²Ðι켣µÄÇ°ÌáÏ£¬»­³öÁâÐÎËùÖ¸³öµÄ·½Ïò¡£ÁâÐεÄλÖúÍÖ¸ÏòÊý¾Ý¶¼¿ÉÒԵõ½¡£ÁâÐεÄÑÕÉ«´ÓÉdz´ú±í×Å ...

function Untitled2
clear;clc;close all
R=10;
x0=linspace(-R,R,50);
x=[x0,fliplr(x0)];
y=[sqrt(R^2-x0.^2),-fliplr(sqrt(R^2-x0.^2))];


theta=atan(y./x);

figure1=figure('color','w');plot(x,y);hold on;axis equal
myr=linspace(0,1,length(x));
myg=linspace(0,0,length(x));
myb=linspace(1,0,length(x));
for jj=1:length(x)   
    plotSquare(x(jj),y(jj),theta(jj),1,[myr(jj),myg(jj),myb(jj)]);
end

function plotSquare(x0,y0,theta,a,mycolor)
xs=[-a/2,0,a/2,0,-a/2];
ys=[0,-a/4,0,a/4,0];
xy=[cos(theta),-sin(theta);
    sin(theta),cos(theta)]*[xs;ys] ;

fill(xy(1,:)+x0,xy(2,:)+y0,mycolor)