MATLAB绘制爱心

2022-01-17 05:31:06 阅读：1584 来源： 互联网

标签：val title 0.1 爱心 grid MATLAB 绘制 sin axis

以下为爱心的MATLAB源代码与对应图片

代码开源无误，皆可复现

love1

 clear;clc
 b=0:0.01*pi:2*pi;
 a=2;
 r=a*(1-sin(b));
 h=polarplot(b,r,'r');
 box on
 set(h,'LineWidth',3);
 title({'公式1';'笛卡尔表白：p=a(1-sin(\theta))'});

love2

 clear;clc
 a=ezplot('x^2-abs(x)*y+y^2-15');
 axis square
 grid on
 set(a,'color','r','LineWidth',3);
 title({'公式2';'x^2-|x|y+y^2-15=0'});

love3

 clear;clc
 f=@(x,y,z)(x.^2+2.*y.^2+z.^2-1).^3-x.^2.*z.^3-0.1.*y.^2.*z.^3;
 [x,y,z]=meshgrid(linspace(-3,3));
 val=f(x,y,z);
 [p,v]=isosurface(x,y,z,val,0);
 patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');
 view(3);
 title({'公式3';'(x^2+2y^2+z^2-1)^3-x^2z^3-0.1y^2z^3=0'});
 box on;
 grid on;
 axis equal;

love4

 clear;clc
 [x,y]=meshgrid(-5:0.01:5,-4:0.01:5);
 v=x.^2-abs(x).*y+y.^2+10./abs(5*x+sin(5*y))<15;
 contourf(x,y,v);
 axis square;
 title({'公式4';'x^2-|x|y+y^2+10/|5x+sin(5y)|<15'});
 colormap(jet);

love5

 clear;clc
 b=ezplot('x^2+(y-(x^2)^(1/3))^2=13');
 set(b,'color','g','LineWidth',5);
 title({'公式5';'(x^2+(y-(x^2)\^(1/3))^2=9'});
 axis square;
 grid on
 axis([-5,5,-5,6])

love6

 clear;clc
 [x,y,z]=meshgrid(linspace(-1.3,1.3));
 val=(x.^2+2*y.^2+z.^2-1).^3-x.^2.*z.^3-0.1*y.^2.*z.^3;
 isosurface(x,y,z,val,0)
 title({'公式6';'(x^2+2y^2+z^2-1)^3-x^2z^3-0.1y^2z^3=0'});
 axis equal
 grid on
 box on;
 view(-10,24);

love合集

 clear;clc
 subplot(2,3,1)
 b=0:0.01*pi:2*pi;
 a=2;
 r=a*(1-sin(b));
 h=polarplot(b,r,'r');
 box on
 set(h,'LineWidth',3);
 title({'公式1';'p=a(1-sin(\theta))'});
 
 subplot(2,3,2)
 a=ezplot('x^2-abs(x)*y+y^2-15');
 axis square
 grid on
 set(a,'color','r','LineWidth',3);
 title({'公式2';'x^2-|x|y+y^2-15=0'});
 
 subplot(2,3,3)
 f=@(x,y,z)(x.^2+2.*y.^2+z.^2-1).^3-x.^2.*z.^3-0.1.*y.^2.*z.^3;
 [x,y,z]=meshgrid(linspace(-3,3));
 val=f(x,y,z);
 [p,v]=isosurface(x,y,z,val,0);
 patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');
 view(3);
 title({'公式3';'(x^2+2y^2+z^2-1)^3-x^2z^3-0.1y^2z^3=0'});
 box on;
 grid on;
 axis equal;
 
 subplot(2,3,4)
 [x,y]=meshgrid(-5:0.01:5,-4:0.01:5);
 v=x.^2-abs(x).*y+y.^2+10./abs(5*x+sin(5*y))<15;
 contourf(x,y,v);
 axis square;
 title({'公式4';'x^2-|x|y+y^2+10/|5x+sin(5y)|<15'});
 colormap(jet);
 
 subplot(2,3,5)
 b=ezplot('x^2+(y-(x^2)^(1/3))^2=13');
 set(b,'color','g','LineWidth',5);
 title({'公式5';'(x^2+(y-(x^2)\^(1/3))^2=9'});
 axis square;
 grid on
 axis([-5,5,-5,6])
 
 subplot(2,3,6)
 [x,y,z]=meshgrid(linspace(-1.3,1.3));
 val=(x.^2+2*y.^2+z.^2-1).^3-x.^2.*z.^3-0.1*y.^2.*z.^3;
 isosurface(x,y,z,val,0)
 title({'公式6';'(x^2+2y^2+z^2-1)^3-x^2z^3-0.1y^2z^3=0'});
 axis equal
 grid on
 box on;
 view(-10,24);

标签：val,title,0.1,爱心,grid,MATLAB,绘制,sin,axis
来源： https://www.cnblogs.com/bloghjx/p/15811992.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9