MATLAB--绘图4

x=[0 1 1 0 0];
y=[0 0 1 1 0];
t=[1 1 1 1 1];
plot(x,y)
h=patch(x,y,'b');
h.FaceAlpha=0.5;
hold on
% axis([-4 4 -4 4])
% axis equal
% M=[1 0 2;
% 0 1 3;
% 0 0 1];
% M1=power(M,1/50);
axis([-6 6 -6 6])
axis equal
M=[cos(pi) -sin(pi) 2;
sin(pi) cos(pi) 3;
0 0 1];
M1=M^(1/50);
for i=1:50
T=M1^i*[x;y;t];
plot(T(1,:),T(2,:))
h1=patch(T(1,:),T(2,:),'r');
h1.FaceAlpha=0.5;
hold on
% axis([-4 4 -4 4])
axis([-6 6 -6 6])
axis equal
pause(0.1)
end

%% 旋转的五角星
clear
clc
point1=[0 1];
point2=[cos(18/180*pi),sin(18/180*pi)];
point3=[cos(54/180*pi),-sin(54/180*pi)];
point4=[-cos(54/180*pi),-sin(54/180*pi)];
point5=[-cos(18/180*pi),sin(18/180*pi)];
x=[point1(1),point3(1),point5(1),point2(1),point4(1),point1(1)];
y=[point1(2),point3(2),point5(2),point2(2),point4(2),point1(2)];
h1=plot(x,y)
hold on
axis equal
axis([-2 2 -2 2])
plot([0 0],[-1 1.5],':b','Linewidth',1)
%delete(h1)
%%
% for i=i:100
% M=[cos(10*pi/100*i) 0;
% 0 1;];
% wx=M*[x;y];
% h=plot(wx(1,:),wx(2,:))
% hold on
% axis([-2 2 -2 2])
% pause(0.1)
% delete(h)%防止原来的图形停留
% end
%% 循环前不能删除h1，否则就没有图形了
for i=i:100
M=[cos(10*pi/100*i) 0;
0 1;];
wx=M*[x;y];
h1.XData=(wx(1,:));%属性里有XData，YData
h1.YData=(wx(2,:));
hold on
axis([-2 2 -2 2])
pause(0.1)
end

%% 小石头
t=0:0.01:0.3;
v=10;
g=-9.8;
x=v*t;
y=1/2*g*t.^2;
comet(x,y,0);%只写y则按计算机给的坐标值。若写x，y则按实际的x的范围绘图,写0后无小尾巴，【0 1】之间，越大尾巴越延后

%% comet3()
close
t=0:0.001:1;
r=8*t;
sita=2*pi*t*4;%仰角
fai=-2*pi*t*8;%方位角
[x,y,z]=sph2cart(fai,pi/2-sita,r);%在matlab中球坐标的参数与数学中的坐标系有些区别，sita为其余角
comet3(x,y,z,0.2)

%% animatedline
clear
clc
ax1=subplot(2,1,1);
axis([0 2*pi -1 1])
ax2=subplot(2,1,2);
axis([0 2*pi -1 1])
h=animatedline(ax1);
%h=animatedline(ax1,x(1:50000),y(1:50000));%先绘制前50000点
h1=animatedline(ax2);
h1.MaximumNumPoints=100;
x=linspace(0,2*pi,100000);
y=sin(x);
for i=1:100000
%for i=50001:100000
addpoints(h,x(i),y(i));
drawnow limitrate
end
x1=linspace(0,2*pi,1000);
y1=cos(x1);
for i=1:1000
addpoints(h1,x1(i),y1(i));
drawnow
end
[xx,yy]=getpoints(h);%得到数据
clearpoints(h1)%清除数据

%% drawnow
figure
hold on
for i=2:1000
x(i)=i;
y(i)=i^2;
plot([x(i-1),x(i)],[y(i-1),y(i)])
drawnow limitrate
end

%% 画图
r=3;
t=linspace(0,2*pi,200);
x=r*cos(t);
y=r*sin(t);
plot(x,y);
%% 产生线
l=8;%发生线的长度
fx=[r r];
fy=[0 l];
hold on
h=plot(fx,fy,'r','linewidth',1);
h1=plot(r,0,'o','markersize',3);
axis([-4*r 4*r -4*r 4*r])
axis equal
fx=[0;0;1];
fy=[0;1;1];
v=VideoWriter('jiankaixian,avi');
open(v)
for i=linspace(0,pi,500)
M1=[1 0 r*cos(i);
0 1 r*sin(i);
0 0 1;];
M2=[cos(i) -sin(i) 0;
sin(i) cos(i) 0;
0 0 1;];
l1=i*r;
detax=l1*sin(i);
detay=-l1*cos(i);
M3=[1 0 detax;
0 1 detay;
0 0 1;];
fxx=M3*M1*M2*fx;
fyy=M3*M1*M2*fy;
axis([-4*r 4*r -4*r 4*r])
axis equal
plot(fxx(1),fxx(2),'o','markersize',3)
h1.XData=fxx(1);
h1.YData=fxx(2);
h.XData=[fxx(1),fyy(1)];
h.YData=[fxx(2),fyy(2)];
M=getframe;
M.cdata=imresize(M.cdata,[343,434]);
pause(0.01)
writeVideo(v,M);
end
close(v)

猜你喜欢