用matlab画椭圆??

用matlab画椭圆??,第1张

可以用椭圆的参数方程嘛。

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这个词条。

欢迎分享,转载请注明来源:内存溢出

原文地址:https://54852.com/langs/11676256.html

(0)
打赏 微信扫一扫微信扫一扫 支付宝扫一扫支付宝扫一扫
上一篇 2023-05-17
下一篇2023-05-17

发表评论

登录后才能评论

评论列表(0条)

    保存