24小时热门版块排行榜

返回列表

langdd007

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 37
在线: 8.6小时
虫号: 4073026
注册: 2015-09-15
专业: 核物理

[求助] matlab怎么带入一堆数值？

t=[0:0.1:2*pi];

x2= (320.*cos(t).^2 + 320.*sin(t).^2 - (sin(t).*(55.*sin(t) - 4.*cos(t).*(-(25.*(64.*cos(t).^2 + 64.*sin(t).^2 - 1).*(64.*cos(t).^2 + 64.*sin(t).^2 - 121))./16).^(1/2) + 320.*sin(t).^3 + 320.*cos(t).^2.*sin(t)))./(cos(t).^2 + sin(t).^2) + 55)./(8.*cos(t));

y2=(55.*sin(t) + 4.*cos(t).*(-((64.*cos(t).^2 + 64.*sin(t).^2 - 121).*(1600.*cos(t).^2 + 1600.*sin(t).^2 - 25))./16).^(0.5) + 320.*sin(t).^3 + 320.*cos(t).^2.*sin(t))./(8.*cos(t).^2 + 8.*sin(t).^2);

y3=(3925.*y2 - (12188224018053065.*x2)./2199023255552 + 24000.*(((6570214509731731.*x2.^3)./68719476736 - (3900231685776981.*x2.^3.*y2)./4398046511104 - 444.*x2.^4 + 68.*x2.^2.*y2.^2 + 142560.*x2.^2.*y2 - 6811425.*x2.^2 - (5179995207672553.*x2.*y2.^2)./549755813888 - (7498614386106867.*x2.*y2)./1073741824 + (661280917803823.*x2)./4194304 + 326400.*y2.^2 + 94464000.*y2 + 206280000)./(x2.^4 - 160.*3.^(1/2).*x2.^3 + 28800.*x2.^2 - (5713230008462375.*x2)./4294967296 + 23040000)).^(1/2) + 80.*x2.^2 + 5.*x2.^2.*(((6570214509731731.*x2.^3)./68719476736 - (3900231685776981.*x2.^3.*y2)./4398046511104 - 444.*x2.^4 + 68.*x2.^2.*y2.^2 + 142560.*x2.^2.*y2 - 6811425.*x2.^2 - (5179995207672553.*x2.*y2.^2)./549755813888 - (7498614386106867.*x2.*y2)./1073741824 + (661280917803823.*x2)./4194304 + 326400.*y2.^2 + 94464000.*y2 + 206280000)./(x2.^4 - 160.*3.^(1/2).*x2.^3 + 28800.*x2.^2 - (5713230008462375.*x2)./4294967296 + 23040000)).^(1/2) - 80.*3.^(1/2).*x2.*y2 - 400.*3.^(1/2).*x2.*(((6570214509731731.*x2.^3)./68719476736 - (3900231685776981.*x2.^3.*y2)./4398046511104 - 444.*x2.^4 + 68.*x2.^2.*y2.^2 + 142560.*x2.^2.*y2 - 6811425.*x2.^2 - (5179995207672553.*x2.*y2.^2)./549755813888 - (7498614386106867.*x2.*y2)./1073741824 + (661280917803823.*x2)./4194304 + 326400.*y2.^2 + 94464000.*y2 + 206280000)./(x2.^4 - 160.*3.^(1/2).*x2.^3 + 28800.*x2.^2 - (5713230008462375.*x2)./4294967296 + 23040000)).^(1/2) + 227000)/(2.*x2.^2 - 160.*3.^(1/2).*x2 + 2.*y2.^2 - 160.*y2 + 12800);

代码如上，我想先设置参数t，然后求出来x2，y2的值（从0到2*pi的对应值）
然后再把x2，y2的对应值导入y3的等式里的x2，y2项，
但是得到的答案只有一个虚数4.8593 +30.8495i
首先不说会得到虚数这么神经的答案，为什么答案只出来了一个？

这里y3是用一个方程((5675-y3*(80-2*y2))/(80*3^0.5-2*x2)-40*3^0.5)^2+(y3-40)^2=45^2来求解，但是为什么得到的y3解会含有y3？而不是只有x2，y2，t的解？
syms t x3 y3;
>> x2= (320*cos(t)^2 + 320*sin(t)^2 - (sin(t)*(55*sin(t) - 4*cos(t)*(-(25*(64*cos(t)^2 + 64*sin(t)^2 - 1)*(64*cos(t)^2 + 64*sin(t)^2 - 121))/16)^(1/2) + 320*sin(t)^3 + 320*cos(t)^2*sin(t)))/(cos(t)^2 + sin(t)^2) + 55)/(8*cos(t));
>> y2=(55*sin(t) + 4*cos(t)*(-(25*(64*cos(t)^2 + 64*sin(t)^2 - 1)*(64*cos(t)^2 + 64*sin(t)^2 - 121))/16)^(1/2) + 320*sin(t)^3 + 320*cos(t)^2*sin(t))/(8*(cos(t)^2 + sin(t)^2));

>> [y3]=solve('((5675-y3*(80-2*y2))/(80*3^0.5-2*x2)-40*3^0.5)^2+(y3-40)^2=45^2');

回复此楼

» 猜你喜欢

1楼 2016-06-14 16:33:13

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 langdd007 的主题更新

返回列表