| 7 | 1/1 | 返回列表 |
| 查看: 3491 | 回復(fù): 6 | |||
shengjunjie木蟲 (正式寫手)
|
[交流]
如何在matlab中繪制含有三個參數(shù)方程的三維圖形 已有1人參與
|
|
如下面公式 對matlab不是很熟悉,想根據(jù)這個公式繪制如下三維圖,其中c是常數(shù),可以去1或者100、1000等,該如何編程,謝謝 [ Last edited by shengjunjie on 2011-4-15 at 09:33 ] |

|
【1】ezimplot3 存儲該M文件,并將文件路徑設(shè)置為matlab當(dāng)前路徑 function h = ezimplot3(fun,domain,n,color) % EZIMPLOT3 Easy to use 3D implicit plotter. % EZIMPLOT3(FUN) plots the inline function FUN(X,Y,Z) = 0 over the % default domain -2*PI < X < 2*PI, -2*PI < Y < 2*PI, -2*PI < Z < 2*PI. % % EZIMPLOT3(FUN,DOMAIN)plots FUN over the specified DOMAIN instead of the % default domain. DOMAIN can be the vector [XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX] % or the vector [A,B] (to plot over A < X < B, A < Y < B, A < Z < B). % % EZIMPLOT3(...,N) plots FUN over the default domain using an N-by-N grid. % The default value for N is 60. % % Example % Plot x^3+exp(y)-cosh(z)=4 % % via a string: f='x^3+exp(y)-cosh(z)-4' % ezimplot3(f) % % via a vectorized function handle: f = @(x,y,z) x.^3+exp(y)-cosh(z)-4 % ezimplot3(f) % % Note: this function do not use the ezgraph3 standard, like ezsurf, ezmesh % ,etc, does. Because of that, ezimplot3 only tries to imitate that % interface. A future work must be to modify the ezgraph3 to include a % routine for implicit surfaces based on this file % % Inspired by works of: A.Jutan UWO 02-02-98 ajutan@julian.uwo.ca % Made by: G.Morales UC 03-20-09 gmorales@uc.edu.ve % if nargin == 1 domain = [-2*pi, 2*pi]; % default domain: -2*pi < xi < 2*pi n = 60; % default grid size elseif nargin == 2 n = 60; % just default grid end if size(domain,2) == 2 domain = repmat(domain,1,3); %domain repeated in all variables end xm = linspace(domain(1), domain(2), n); % generating the volume data ym = linspace(domain(3), domain(4), n); zm = linspace(domain(5), domain(6), n); [x,y,z] = meshgrid(xm, ym, zm); if ischar(fun) fun = inline(vectorize(fun)); % making sure string "fun" is vectorized fvalues = feval(fun,x,y,z); % evaluating "fun" in domain elseif isa(fun,'function_handle') fvalues = fun(x,y,z); % evaluating "fun" in domain fun = char(fun); fun = fun(9:end); % pre-formatting of graph title end h = patch(isosurface(x,y,z, fvalues, 0)); %"patch" handles the structure %sent by "isosurface" isonormals(x,y,z,fvalues,h) %Recalculating the isosurface normals based on %the volume data set(h,'FaceColor','red','EdgeColor','none'); xlabel('x');ylabel('y');zlabel('z');% naming the axis alpha(0.8) % adjusting for some transparency grid on view(3) axis equal camlight lighting gouraud title([strrep(char(fun),'.','') ' = 0']); % graph title without "." 【2】在matlab命令窗口中輸入 f='x-y*(y+1)*(z^(-2)-z^(-8))*x^2+(y+1)^2*(z^2-z^(-4))^3/100' 這里c我給取了10,所以c^2=100 |


木蟲 (正式寫手)

木蟲 (正式寫手)


木蟲 (正式寫手)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +4 | yy要上岸呀 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +3 | beefly 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +5 | 為學(xué)666 2026-03-16 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|