[問題] 圓柱座標畫立體圖，邊緣鋸齒狀如何修正已回收
Matlab 如何更改座標以圓柱畫圖
小弟知道的 Matlab畫圖 surf等指令都是以直角座標下去畫圖，對於
圓柱座標要表示就比較困難些，而在圖的邊緣如果不是矩形的情況，
硬將邊緣消掉卻又會使邊緣不圓滑，小弟想要問的是如何畫出邊緣圓滑
的立體圖，附上程式碼請前輩指點，謝謝
% 根據 Siedel aberration 的方程式，畫出初階像差的分佈
% 初階像差方程式
% thi = -1/4*B*rho^4 -c*y0^2*rho^2*(cos(theta))^2 -1/2*D*y0^2*rho^2
% +E*yo^3*rho*cos(theta) +F*y0*rho^3*cos(theta)
% B,C,D,E,F,y0 are constants
% 像差簡稱根據codev的命名 SA(spherical abettation); TCO(tengtial coma);
% TAS(tangtial astigmatism); DST(distorsion); PTZ(field curvature)
% 清除參數
clc;
clear;
close all;
% 基本設定輸入
sampling=81;
% 產生基本直角座標
x = ones(sampling,1)*linspace(-1,1,sampling);
y = linspace(-1,1,sampling)'*ones(1,sampling);
[theta,rho] = cart2pol(x,y);
for i=1:sampling
for j=1:sampling
if rho(i,j) > 1
rho(i,j) = nan;
end
end
end
% 計算normalized的像差
SA = rho.^4;
TCO = rho.^3.*cos(theta);
TAS = rho.^2.*(cos(theta)).^2;
PTZ = rho.^2;
DST = rho;
x = rho.*cos(theta); % spherical coordinate equations
y = rho.*sin(theta);
t=[zeros(1,20) ones(1,20)];
[X,Y,Z] = cylinder(t);
surf(x,y,SA)
邊緣若切掉就會有鋸齒狀出現，要如何以改座標設定的方式修正呢？
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.46.89
討論串 (同標題文章)
完整討論串 (本文為第 1 之 2 篇)：
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章