[matlab] 用 matlab 做數(shù)據(jù)

[復(fù)制鏈接]

shouce

電梯直達(dá)

1^#

發(fā)表于 2016-2-1 21:01:08 | 只看該作者 |倒序瀏覽 |閱讀模式

本帖最后由 shouce 于 2016-2-1 21:19 編輯

function ABCk1 %圓弧包絡(luò)線  陰推陽 20160131
clc
clear
t=linspace(-0.2401727,-0.65911,200); %原參數(shù)范圍
%t=linspace(-1.608786,-1.35,100);
i=4/6;k=i+1;
a=95.937985;b=8.801204;A=100;r=60;
q=-(t - asin((b*k*r*cos(t) + a*k*r*sin(t))/(A*r)))/i;  %
x=A-(A*cos(q) - cos(k*q).*(a - r*cos(t)) - sin(k*q).*(b + r*sin(t)));
y=cos(k*q).*(b + r*sin(t)) - sin(k*q).*(a - r*cos(t)) + A*sin(q);
z=0*t;
plot(x,y,'-go') %方便查看點的間距
ABk1=[x',y',z'];
vpa(x',8)  % 輸出8位數(shù)值  20160201修改
save ABk1.txt ABk1 -ascii
% 用 matlab 做數(shù)據(jù)

function ABCk2 %圓弧包絡(luò)線
clc
clear
t=linspace(-0.6981317,0.65911,200); %原參數(shù)范圍
%t=linspace(-1.608786,-1.35,100);
i=4/6;k=i+1;
a=40.60035549;b=-34.06774331;A=100;r=10;
q=(t + asin((k*(b*cos(t) - a*sin(t)))/A))/i;
x=A-(A*cos(q) - cos(k*q).*(a + r*cos(t)) - sin(k*q).*(b + r*sin(t)));
y=cos(k*q).*(b + r*sin(t)) - sin(k*q).*(a + r*cos(t)) + A*sin(q);
z=0*t;
plot(x,y,'-go') %方便查看點的間距
ABk2=[x',y',z'];
save ABk2.txt ABk2 -ascii
% 用 matlab 做數(shù)據(jù)

function ABCk3 %圓弧包絡(luò)線
clc
clear
t=linspace(-1.59196897,0.30879537,500); %原參數(shù)范圍
%t=linspace(-1.608786,-1.35,100);
i=4/6;k=i+1;
a=57.06675156;b=18.20428153;A=100;r=3.1;
q=(t + asin((k*(b*cos(t) - a*sin(t)))/A))/i;
x=A-(A*cos(q) - cos(k*q).*(a + r*cos(t)) - sin(k*q).*(b + r*sin(t)));
y=cos(k*q).*(b + r*sin(t)) - sin(k*q).*(a + r*cos(t)) + A*sin(q);
z=0*t;
plot(x,y,'-go') %方便查看點的間距
ABk3=[x',y',z'];
save ABk3.txt ABk3 -ascii
% 用 matlab 做數(shù)據(jù)

function, matlab