matlab 7.0三维心形函数 求matlab的三维心形方程的代码

\u7528MATLAB\u7ed8\u5236\u5fc3\u5f62\u51fd\u6570

\u672c\u89c6\u9891\u5c55\u793a\u5982\u4f55\u7528matlab\u7ed8\u5236\u6563\u5706\u72b6\u6001\u56fe\uff0c\u53ef\u7528\u4e8e\u76f8\u5173\u79d1\u7814\u6570\u636e\u7ed8\u56fe\uff01

\u7ed9\u4f60\u4e00\u6bb5\u5f88\u5168\u7684\u4ee3\u7801\uff0c\u81ea\u5df1\u8fd0\u884c\u770b\u770b
clf;
clc;
clear;
close all;
figure(1);
%\u7528\u6781\u5750\u6807\u7684\u5f62\u5f0f\u5b9e\u73b0
%r=a(1-sin\u03b8)
x=-pi:pi/1000:pi;
y=2*(1-sin(x));
polar(x,y,'r-');
title('\fontsize{16}To SmallWhite');
xlabel('\fontsize{16}\copyrightDingQian');
figure(2);
%\u7528\u9690\u51fd\u6570\u7684\u5f62\u5f0f\u5b9e\u73b0
%(x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0 -3<=x,y,z<=3
[x,y,z]=meshgrid(linspace(-3,3)); %\u505a\u51fa\u7f51\u683cmeshgrid
p=(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80)*y.^2.*z.^3; %\u5b9e\u73b0\u7ed3\u679c\u7684\u8868\u8fbe
isosurface(x,y,z,p,0); %\u4f7f\u7528\u51fd\u6570\u505a\u56femesh surf
axis equal;
axis off;
view(0,45); %\u89c6\u89d2\u7684\u63a7\u5236
colormap([1 0 0]); %\u7ed8\u56fe\u989c\u8272\u7ea2\u8272
brighten(0.5); %\u589e\u4eae
camlight right; %\u5149\u6e90\u4f4d\u7f6e
lighting phong; %\u5149\u7167\u6a21\u5f0f
title('\fontsize{16}To SmallWhite');
xlabel('\fontsize{16}\copyrightDingQian');
text(-0.37,0.75,1,'\fontsize{16}\copyrightDingQian');
figure(3);
isosurface(x,y,z,p,0);
view(0,45);
colormap([1 0.2 0.2]);
axis off;
axis equal;
title('\fontsize{16}To SmallWhite');
text(-0.37,0.75,1.20,'\fontsize{16}\copyrightDingQian');
figure(4);
isosurface(x,y,z,p,0);
view(-10,24);
colormap([1 0.2 0.2]);
axis off;
axis equal;
title('\fontsize{16}To SmallWhite');
text(-0.37,0.75,1.20,'\fontsize{16}\copyrightDingQian');
figure(5);
isosurface(x,y,z,p,0);
view(-10,24);
colormap([1 0 0]);
axis off;
axis equal;
title('\fontsize{16}To SmallWhite');
text(-0.37,0.75,1.20,'\fontsize{16}\copyrightDingQian');
figure(6);
isosurface(x,y,z,p,0);
view(1,11);
colormap([1 0 0]);
axis off;
axis equal;
title('\fontsize{16}To SmallWhite');
\u8fd9\u662f\u5f53\u5e74\u6211\u5ba4\u53cb\u5199\u7ed9\u4ed6\u5973\u670b\u53cb\u7684\u3002\u503c\u5f97\u4e00\u770b

% 二维心形图
figure(1)
ezplot('(x^2-1+y^2)^3-x^2*y^3 = 0')
axis([-2.5 2.5 -2 2])

% 三维心形图
figure(2)
n=100;
x=linspace(-3,3,n);
y=linspace(-3,3,n);
z=linspace(-3,3,n);
[X,Y,Z]=ndgrid(x,y,z);
F=((-(X.^2) .* (Z.^3) -(9/80).*(Y.^2).*(Z.^3)) + ((X.^2) + (9/4).* (Y.^2) + (Z.^2)-1).^3);
isosurface(F,0)
lighting phong
caxis
axis equal
colormap('flag');
view([55 34]);

扩展阅读：matlab7.0安装包下载 ... matlab7.0中文版界面 ... mathematica心形3d图 ... 如何安装matlab7.0版本 ... matlab mesh ... matlab7.0破解版安装教程 ... 用matlab画跳动的爱心三维 ... matlab笛卡尔心形图 ... matlab7.0安装详细教程 ...