matlab画简单图

clc;clear all;close all;
%% 
m=0,n=0;%圆心
r=6;%半径
h=zeros(1,100);h2=h;h3=h;
w=zeros(1,100);w2=w;w3=w;
for i=0:300
    
    w(i+1)=m+r*cos(i/100*2*pi);
    h(i+1)=r*sin(i/100*2*pi);
    %画圆
    x=m-r:0.1:m+r;
    y1=sqrt(r.^2-(x-m).^2);subplot(121);plot(x,y1,'b');hold on;grid on;
    y2=-sqrt(r.^2-(x-m).^2);subplot(121);plot(x,y2,'b');
    set(gca,'XLim',[-15,15]);set(gca,'YLim',[-15,15]);
    %画点
    subplot(121);plot([0,w(i+1)],[0,h(i+1)]);%plot(w,h);
        
    m2=w(i+1);
    n2=h(i+1);%圆心
    r2=3;%半径
    w2(i+1)=m2+r2*cos(i*3/100*2*pi);
    h2(i+1)=n2+r2*sin(i*3/100*2*pi);
    x=m2-r2:0.1:m2+r2;
    y1=n2+sqrt(r2.^2-(x-m2).^2);subplot(121);plot(x,y1,'b');hold on;grid on;
    y2=n2-sqrt(r2.^2-(x-m2).^2);subplot(121);plot(x,y2,'b');
    subplot(121);plot([m2,w2(i+1)],[n2,h2(i+1)]);%plot(w2,h2);
    
    m3=w2(i+1);
    n3=h2(i+1);%圆心
    r3=1.5;%半径
    w3(i+1)=m3+r3*cos(i*5/100*2*pi);
    h3(i+1)=n3+r3*sin(i*5/100*2*pi);
    x=m3-r3:0.1:m3+r3;
    y1=n3+sqrt(r3.^2-(x-m3).^2);subplot(121);plot(x,y1,'b');hold on;grid on;
    y2=n3-sqrt(r3.^2-(x-m3).^2);subplot(121);plot(x,y2,'b');
    subplot(121);plot([m3,w3(i+1)],[n3,h3(i+1)]);plot(w3,h3,'r+');
    
    %画右边输出信号线
    subplot(122);plot(h(1:i+1)+h2(1:i+1)+h3(1:i+1),'r*-');grid on;
    set(gca,'XLim',[0,300]);set(gca,'YLim',[-30,30]);
    
    %暂停
    pause(0.1);
    clf;
end
close;


%%