Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 基础科学 > 物理学
物理学 A discussion board for Physics.
回复
 
主题工具 显示模式
旧 2009-04-11, 05:24 PM   #1
初级会员
 
注册日期: 2009-03-22
年龄: 35
帖子: 8
感谢他人: 2
有 0 帖获得 0 感谢
声望力: 0
Lixiney 正向着好的方向发展
警告 请MATLAB高手帮我解这道隐函数方程

急!请MATLAB高手帮我解这道隐函数方程
已知角1相应度数,求角3相应度数,我解了很久都解不出来,好急!
麻烦高手解围!


其中L1==70cm,L2=130cm,L3=90cm,L4=150cm, W=70rad/s

此帖于 2009-04-11 05:37 PM 被 Lixiney 编辑。
Lixiney 当前离线   回复时引用此帖
旧 2009-04-11, 06:40 PM   #2
hrt
初级会员
 
注册日期: 2009-03-11
年龄: 35
帖子: 9
感谢他人: 0
有 5 帖获得 5 感谢
声望力: 0
hrt 正向着好的方向发展
默认 回复: 急!请MATLAB高手帮我解这道隐函数方程

可利用fsolve函数求解;
先编写函数M文件,如下:
function y=f(x,a) %x是数组,存储theta1、theta2,a是theta1的值
y(1)=x(1)-a;
l1=0.7;l2=1.3;l3=0.9;l4=1.5;
y(2)=l1*cos(x(1))+l2*sqrt(1-((l3*sin(x(2))-l1*sin(x(1)))/l2)^2)-l3*cos(x(2))-l4;

再在matlab命令行中输入:
>> x0=[0,0.2];
>> [x,fv,ef]=fsolve(@f,x0,[],0)
即可得到:
x =

0 1.4033


fv =

1.0e-007 *

0 0.5511

即:theta1=0时,theta2=1.4033(rad)
其它情况依此即可,要注意初值x0的选取,否则可能不收敛。
~ ~
hrt 当前离线   回复时引用此帖
旧 2009-04-12, 12:17 AM   #3
初级会员
 
注册日期: 2009-03-22
年龄: 35
帖子: 8
感谢他人: 2
有 0 帖获得 0 感谢
声望力: 0
Lixiney 正向着好的方向发展
默认 回复: 急!请MATLAB高手帮我解这道隐函数方程

非常的谢谢您!很好,很开心!
Lixiney 当前离线   回复时引用此帖
旧 2009-05-01, 09:52 AM   #4
高级会员
 
注册日期: 2008-11-07
住址: 湖南长沙
帖子: 233
感谢他人: 0
有 24 帖获得 25 感谢
声望力: 18
laosam280 正向着好的方向发展
默认 回复: 急!请MATLAB高手帮我解这道隐函数方程

我把周期一个周期内计算了101个点,并绘制了图像。应该比较直观。txt里有这些点的结果和误差。
上传的图像
文件类型: jpg map of this function.jpg (30.5 KB, 10 次查看)
上传的附件
文件类型: txt solution of this equation.txt (3.1 KB, 10 次查看)
__________________
坚持就是胜利,努力就有奇迹。
laosam280 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码
Trackbacks are 禁用
Pingbacks are 禁用
Refbacks are 禁用



所有时间均为北京时间。现在的时间是 02:49 AM


Powered by vBulletin
版权所有 ©2000 - 2022,Jelsoft Enterprises Ltd.

SEO by vBSEO ©2009, Crawlability, Inc.