如何按下面给定的方程建模

hrbustmqc

注：a角为0到90度范围的某一已知量，t在一定范围内变动
要求：在nasys中建立该曲线模型

tjqhnu

3# hrbustmqc
finish
/clear      
/title,curve_model
/prep7
*set,R,10                             ! 相当于赋值R=10
*set,pi,acos(-1)                    ! 给pi附上常数3.1415926...
i=1                                        ! 起始关键点号
afa=1                                   ! *do循环的afa起始值（自己设定，
                                               但该曲线方程中afa不能等于0）
t=-10                                     ! *do循环的afa起始值（自己设定）
*do,afa,1,90,5                       !*do开始
*do,t,-10,10,2
x=R*sqrt(1+t*t)*sin(afa*pi/180)*cos(1/sin(afa*pi/180)*(t-atan(t)))
y=R*sqrt(1+t*t)*sin(afa*pi/180)*sin(1/sin(afa*pi/180)*(t-atan(t)))
z=R*sqrt(1+t*t)*cos(afa*pi/180)
k,i,x,y,z                                   !创建关键点
i=i+1
*enddo
*enddo                                 !循环结束
bspline,all                             !将关键点生成样条曲线

tjqhnu

finish
/clear
/title,curve_model
/prep7
*set,R,10
*set,pi,acos(-1)
i=1
afa=1
t=-10
*do,afa,1,90,5
*do,t,-10,10,2
x=R*sqrt(1+t*t)*sin(afa*pi/180)*cos(1/sin(afa*pi/180)*(t-atan(t)))
y=R*sqrt(1+t*t)*sin(afa*pi/180)*sin(1/sin(afa*pi/180)*(t-atan(t)))
z=R*sqrt(1+t*t)*cos(afa*pi/180)
k,i,x,y,z
i=i+1
*enddo
*enddo
bspline,all

hrbustmqc

2# tjqhnu 能详细解释下每部的含义么？

cccjob699

学习了，谢谢

livelife

如果要画的更精确些，可以将
*do,afa,1,90,5                       !*do开始
*do,t,-10,10,2
改为：
*do,afa,1,90,1                       !*do开始
*do,t,-10,10,1

hrbustmqc

谢谢大家的解答，！
我还想请教个的问题
既然给出了曲线精确的表达式，为什么不能直接精确的画出该曲线呢？非要采用关键点拟合的形式。
这样也做不到精确的建模啊，会出现误差。
本人新手，请好心人赐教

saintseraph

民间自有高手啊