
x=-2+2cos(t)
y=6+3sin(t)
t从0到2pi。
具体程序就是:
>> t=linspace(0,2pi,1000);
>> x=-2+2cos(t);
>> y=6+3sin(t);
>> plot(x,y),grid on
画出来的就是这样:
你的第一行是要干嘛
function ececcmaking(a,b)
哪有这么写函数的
function output = myfunctionName(input1,input2)
end
问题比较复杂,值能给你一些参考
syms x1 y1 z1
x0=1;y0=2;z0=3;%椭圆中心坐标
a=1;b=2;c=3;%椭圆轴
f1=solve('(x1-x0)^2/(a^2)+(y1-y0)^2/(b^2)+(z1-z0)^2/(c^2)=1')
subs(f1)
帮助里面的solve里面有个解符号变量方程的例子,可以指定未知数,感觉这个里面你能用到,
说面说的是变形、
但是,,好像画图不是用这个啊,,椭圆参数方程式,x=acos(t);y=bsin(t);用这个画椭圆很方便,至于椭圆柱体,这需要在上面的方程上加上对应的z就行了
比如t=0:01:2pi;
x=acos(t);y=bsin(t);
好吧下面你是程序
t=0:01:2pi;
x=2cos(t);y=3sin(t);z1=linspace(0,5,length(x));
z=[meshgrid(z1)]';
x=meshgrid(x);
y=meshgrid(y);
surf(x,y,z)
axis equal
for k=1:4
view(-375,10k)
pause %按任意键
end
[x,y]=meshgrid(x,y);
z=0:01:5;
surf(x,y,z);
K(b)是第一类完全椭圆积分可以用ellipticCK(b)函数求解,E(b)是第二类完全椭圆积分可以用ellipticCE(b)函数求解。那么含有完全椭圆积分的方程,就可以用vpasolve函数得到其b解。解决主要代码:
syms b
b=vpasolve(ellipticCE(b)-(1-4b^2)ellipticCK(b)/2==08)
其结果为
你可以换个元,y^4=a(cosp)^2,方程两边开根号
于是方程变为a^(1/4)1/2(cosp)^(-1/2)(-sinp)(p')=a^(1/2)sinp
从而,p'=C(一个常数)(cosp)^(1/2)。
这是个单摆方程,解是Jacobi椭圆函数来表示的。
你可以wiki百科看一下pendulum equation这个词条。
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)