社区首页 >专栏 >matlab画心形曲线_笛卡尔心形曲线方程

matlab画心形曲线_笛卡尔心形曲线方程

全栈程序员站长

发布于 2022-11-02 01:58:22

1.2K00

代码可运行

文章被收录于专栏：全栈程序员必看全栈程序员必看

运行总次数：0

代码可运行

MATLAB 心形曲线

基本知识

clc；指令可以清除屏幕，所以你可以通过clc指令来清理屏幕

clc

hold on；指令可以将画的图连起来

hold on

clear；清除之前所留的定义

clear

笛卡尔爱心曲线

ezpolar('1-sin(t)')

a=1;
theta = 0 : 0.01 : 2 * pi;
r = a*(1 - sin(theta));
polar(theta, r, '-r');

特别定制

第一种实现方式

clear
t=-pi:pi/100:pi;

r=abs(t);
x=r.*sin(t);
y=r.*cos(t);
plot(x,y)
title('I love you.')
axis equal

方法二：grid on可以加上网格，可以通过删除下面代码中的grid on删除表格

clear
x=-2:0.01:2;
y=sqrt(2*sqrt(x.^2)-x.^2);
z=asin(abs(x)-1)-pi./2;plot(x,y);
grid on;
hold on;
plot(x,z);
axis equal;

fill语句填色

clear
x=-2:0.01:2;
y=sqrt(2*sqrt(x.^2)-x.^2);
z=asin(abs(x)-1)-pi./2;plot(x,y);
grid on;
hold on;
plot(x,z);
axis equal;
fill(x,y,'r')
fill(x,z,'r')

方案3

t=0:0.1:2*pi;
x=16*sin(t).^3;
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
figure;
axis([-10,10,-10,10]);
plot(x,y)

axis off 可以关闭坐标轴，进行动态绘图

clear
t=0:0.1:2*pi;
x=16*sin(t).^3;
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
figure;
h = plot(x(1),y(1),'r');
axis([-20,20,-20,20]);
axis off

for idx = 2:length(t)
    h.XData(idx) = x(idx);
    h.YData(idx) = y(idx);
    drawnow
end

如果你需要保存到本地的话，你需要进一步添加imwrite保存到本地，delayTime延时画图

clear
t=0:0.1:2*pi;
x=16*sin(t).^3;
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
figure;
h = plot(x(1),y(1),'r');
axis([-20,20,-20,20]);
axis off
[A,map] = rgb2ind(frame2im(getframe),256);
imwrite(A,map,'love.gif','LoopCount',inf,'DelayTime',0.1);
for idx = 2:length(t)
    h.XData(idx) = x(idx);
    h.YData(idx) = y(idx);
    drawnow
%这两句话是为了在本地保存的图片也能显示出来，保存地址默认
    [A,map] = rgb2ind(frame2im(getframe),256);
    imwrite(A,map,'love.gif','WriteMode','append','DelayTime',0.1);  
end

最终实现版本

进一步动态实现（借用movie函数）再进行本地保存改进

https://blog.csdn.net/lpsl1882/article/details/50806694?depth_1-utm_source=distribute.pc_relevant.none-task&utm_source=distribute.pc_relevant.none-task

首先添加函数保存

function ratio=HeartRate(i,heartSpeed)
i=mod(i,heartSpeed);
if i < heartSpeed / 2
    ratio=sin( i/heartSpeed*pi);
elseif i >= heartSpeed / 2 && i < heartSpeed * 3/4
    ratio=2*cos(i/heartSpeed*2/3*pi);
else
    ratio = 0;
end

通过Movie函数进行跳动Heartbeats

clear
gcp=figure;
M=moviein(100,gcp);
heartSpeed=15;
radiusRate=3;

for i = 1:100
    
    axis([-20 20 -20 20]);
axis off
[A,map] = rgb2ind(frame2im(getframe),256);
imwrite(A,map,'1.gif','LoopCount',inf,'DelayTime',0.1);
    cla(gcp);
    x=[];
    y=[];
    
    for t = 0:0.01:2*pi
        r=(sin(t)*sqrt(abs(cos(t))))/(sin(t)+7/5)-2*sin(t)+2;
        r = radiusRate * r;
        x=[x,r*cos(t)*(0.5+0.5*HeartRate(i,heartSpeed))];
        y=[y,r*sin(t)*(0.5+0.5*HeartRate(i,heartSpeed))];
    end
    hold on;
    fill(x,y,'r');
    
    M(:,i)=getframe;
    [A,map] = rgb2ind(frame2im(getframe),256);
    imwrite(A,map,'1.gif','WriteMode','append','DelayTime',0.1);  
end
movie(M,1);