?? trajectory.asv
字號:
clear;
clc;
global y;
global AMa;
global Aa;
global ACX;
global ACY;
global AP;
global APT;
global AXG;
global AXGT;
global AMZWZ1;
global AMZWZ2;
global AMZAF;
global AJZ;
global AJZT;
load -ascii ACX.txt; %讀入數據
load -ascii ACY.txt;
load -ascii AP.txt;
load -ascii APT.txt;
load -ascii AJZ.txt;
load -ascii AJZT.txt;
load -ascii AXG.txt;
load -ascii AXGT.txt;
load -ascii Aa.txt;
load -ascii AMa.txt;
load -ascii AMZAF.txt;
load -ascii AMZWZ1.txt;
load -ascii AMZWZ2.txt;
AE=0.01745329; %角度化弧度
h=0.005;
fo1=fopen('result.dat','wt');
fprintf(fo1,'y8 y1 y2 y3 y4 y5 y6 y7 ALPHA \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%賦初值
y(1)=20.0; %初速
y(2)=18*AE; %初始彈道傾角
y(3)=0; %初始轉速
y(4)=18*AE; %初始俯仰角
y(5)=0; %初始x坐標
y(6)=20; %初始y坐標
y(7)=52.38; %起飛質量
y(8)=0 %初始時間
fprintf(fo1,'%10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f\n',...
y(8),y(1),y(2)*57.3,y(3)*57.3/360,y(4)*57.3,y(5),y(6),y(7),y(4)-y(2));
while y(6)>=0
rk(8,h);
fprintf(fo1,'%10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f\n',y(8),y(1),y(2)*57.3,y(3)*57.3/360,y(4)*57.3,y(5),y(6),y(7),y(4)-y(2));
end
fclose(fo1);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -