plot_ellipse.m

来自「利用ＭＡＴＬＡＢ实现了对测量中的平面网」· M 代码 · 共 24 行

function plot_ellipse(fit2,x,y,a,b,fi,scale)
BETA = pi/ 30;
fe=fi*pi/180;
for j = 1:60
	XE(j) = x + (b*cos(fe)*sin(BETA * j)+a *sin(fe)*cos (BETA * j))*scale;
	YE(j) = y - (b*sin(fe)*sin(BETA * j)-a *cos(fe)*cos (BETA * j))*scale;
	ZE(j) = 0;
end
for K = 1:59
	fprintf(fit2,' 0\n'); fprintf(fit2,'LINE\n');
	fprintf(fit2,' 8\n'); fprintf(fit2,'ELLIPSE\n');
	fprintf(fit2,' 10\n'); fprintf(fit2,'%15.3f\n',XE(K));
	fprintf(fit2,' 20\n'); fprintf(fit2,'%15.3f\n',YE(K));
	fprintf(fit2,' 11\n'); fprintf(fit2,'%15.3f\n',XE(K+1));
	fprintf(fit2,' 21\n'); fprintf(fit2,'%15.3f\n',YE(K+1));
end
fprintf(fit2,' 0\n'); fprintf(fit2,'LINE\n');
fprintf(fit2,' 8\n'); fprintf(fit2,'ELLIPSE\n');
fprintf(fit2,' 10\n'); fprintf(fit2,'%15.3f\n',XE(60));
fprintf(fit2,' 20\n'); fprintf(fit2,'%15.3f\n',YE(60));
fprintf(fit2,' 11\n'); fprintf(fit2,'%15.3f\n',XE(1));
fprintf(fit2,' 21\n'); fprintf(fit2,'%15.3f\n',YE(1));
return