QQ在线咨询
售前咨询热线
15821572819
售后咨询热线
15821572819

matlab心形图大全,几个经典的函数图像,有趣的函数图像,matlab画图


收藏几个经典的函数,用来当作数学优化方法求解的例子。再收藏几个有趣的函数图像。



1 平面心形图像

经过 网上搜索,平面坐标系中,下面函数最像心形图像,而且函数形式简单,不需要分段。

 

x2+(y(x2)1/3)2=9


改变右面的数值可以变化心形图像的大小。


[plain] view plain copy

  1. > ezplot('x^2+(y-(x^2)^(1/3))^2=9')  




还有一个函数:

                         

x2y3+(x2+y21)3=0




[plain] view plain copy

  1. ezplot('-x^2*y^3+(x^2+y^2-1)^3=0',[-1.5,1.5])  




另外一个函数:


17x216|x|y+17y2=200




[plain] view plain copy

  1. ezplot('17*x.^2-16*abs(x).*y+17*y.^2=200'  



另外一个分段函数:

   

f(x)=2x2x2g(x)=2.142|x|



[plain] view plain copy

  1. x=linspace(-2,2,1000);  

  2. y1=sqrt(2*sqrt(x.^2)-x.^2);  

  3. y2=-2.14*sqrt(sqrt(2)-sqrt(abs(x)));  

  4. plot(x,y1,'b',x,y2,'b');  

  5. axis([-2.5,2.5,-3,1.5]);  


另外一个:


x=16(sint)3y=13cost5cos(2t)2cos(3t)cos(4t)





[plain] view plain copy

  1. t=linspace(-6,6,1000);  

  2. x=16*(sin(t)).^3;  

  3. y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);  

  4. plot(x,y);  





2. 立体心形图像

立体心形貌似只能由笛卡尔方程得到:

 

(x2+9y2/4+z21)3x2z39y2z3/80=0


画图时用 isosurface 和 patch  函数:

图形1:


[plain] view plain copy

  1. f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;  

  2. [x,y,z]=meshgrid(linspace(-3,3));  

  3. val=f(x,y,z);  

  4. [p,v]=isosurface(x,y,z,val,0);  

  5. patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');  

  6. view(3);grid on;axis equal;  




图形2:


[plain] view plain copy

  1. f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;  

  2. [x,y,z]=meshgrid(linspace(-1.5,1.5));  

  3. val=f(x,y,z);  

  4. isosurface(x,y,z,val,0);   

  5. axis equal;view(3);colormap([1 0.2 0.2])  







3. 轴对称振荡器函数

   

f(x)=xsin(x)




[plain] view plain copy

  1. ezplot('x*sin(x)')  




4.  刚好有一个极大点,一个极小点的二元函数

   

f(x,y)=xex2y2


    


[plain] view plain copy

  1. ezmesh(@(x,y) x.*exp(-x.^2-y.^2))  




5. peaks 函数(二元高斯分布的概率密度函数)(三个极大点,三个极小点)

   

f(x,y)=3(1x)2ex2(y+1)210(15xx3y5)ex2y213e(x+1)2y2




[plain] view plain copy

  1. f=@(x,y)3*(1-x).^2.*exp(-(x.^2) - (y+1).^2)- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2)- 1/3*exp(-(x+1).^2 - y.^2);  

  2. ezmesh(f);  






注:本文转载自https://blog.csdn.net/robert_chen1988/article/details/52865490,转载目的在于传递更多信息,并不代表本网赞同其观点和对其真实性负责。如有侵权行为,请联系我们,我们会及时删除。
上一篇 下一篇

分享