| 查看: 1367 | 回复: 5 | ||
三星2010木虫 (小有名气)
|
[求助]
请高人指点一个matlab编程问题，谢谢
|
|
如题：主要是通过求解一个非线性方程组来求最值，一些基本的参数及方程如下： taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=500000; t1=5.12; b2=400; E2=10000; t2=20; eql=600; beta=b2*E2*t2/(b1*E1*t1); lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5; lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5; equation1：lamda2/lamda1*tanh(0.5*lamda1*(L-x(1)-x(2)))-(sin(lamda2*x(1))-beta*sin(lamda2*x(2)))/((beta*cos(lamda2*x(2)))-cos(lamda2*x(1)))=0; equation2:lamda2*sinh(L-x(1)-x(2))*tan(lamda2*x(1))+1-lamda1*(cosh(L-x(1)-x(2))+1)=0; 其中x(1)+x(2)<L, L从400变化到1200，L的步长为20，求每一L对应的x(1)和x(2)，最好通过循环输出吧，我用fsolve尝试了一下，并不理想，请高人给以指点，先谢谢啦！ |
回复此楼» 猜你喜欢
看来今天不会放榜了？
已经有9人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
时间戳又变了
已经有18人回复
-
时间戳又变了8-15
已经有27人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 中心差分求解二维热传导的matlab程序
已经有12人回复
- 关于matlab GUI打包生产exe文件，并使其在没有安装matlab的计算机上运行的问题。。
已经有4人回复
- 求解一元三次方程
已经有11人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 请教大家一个反卷积的问题~~帮忙推导或者matlab编程计算~~
已经有5人回复
- 求使用MATLAB语言将植被与土壤精确分割的程序
已经有8人回复
- 关于matlab中循环读取多个excel文件的问题？
已经有3人回复
- matlab 气泡边缘提取 求助高手！！
已经有5人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- 向各位matlab高手求助了，声音处理方面。
已经有3人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 求高人指点 MATLAB 编程求解 洛伦兹的最大李雅普诺夫指数程序！！！！
已经有5人回复
- matlab程序运行总是出现下面的错误，请各位高手指点一下，谢谢了！急 急急
已经有5人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】求助 matlab编程：可否实现清除文本文档数据
已经有5人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】用MATLAB写程序求解简单的正弦函数
已经有4人回复
- 【求助】请教把java程序转化成matlab程序
已经有9人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
2楼2013-12-01 21:33:23
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
3楼2013-12-01 21:59:02
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
4楼2013-12-01 23:34:45
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
xzhdty: 金币+1, 谢谢参与 2013-12-02 07:49:27
三星2010: 金币+50, ★★★很有帮助, 先给一半的金币，验证通过后补齐！ 2013-12-02 14:11:35
三星2010: 金币+50, ★★★★★最佳答案 2013-12-05 22:36:31
xzhdty: 金币+1, 谢谢参与 2013-12-02 07:49:27
三星2010: 金币+50, ★★★很有帮助, 先给一半的金币，验证通过后补齐！ 2013-12-02 14:11:35
三星2010: 金币+50, ★★★★★最佳答案 2013-12-05 22:36:31
|
楼主可尝试编写m文件如下： function xfsolve clear all;clc; global L; for L=400:20:1200 x0=[L/2,L/2]; k=1; xa(k,  =fsolve(@abc,x0)k=k+1; end function f=abc(x) global L; taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=500000; t1=5.12; b2=400; E2=10000; t2=20; eql=600; beta=b2*E2*t2/(b1*E1*t1); lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5; lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5; lamda2=(lamda^2*2*Gf/(deltaf-delta1)/taof)^0.5; f(1)=lamda2/lamda1*tanh(0.5*lamda1*(L-x(1)-x(2)))-(sin(lamda2*x(1))-beta*sin(lamda2*x(2)))/((beta*cos(lamda2*x(2)))-cos(lamda2*x(1))); f(2)=lamda2*sinh(L-x(1)-x(2))*tan(lamda2*x(1))+1-lamda1*(cosh(L-x(1)-x(2))+1); 其中x0是迭代初值，楼主可依情况适当更改，xa(k, 即为要求的x(1)，x(2)。 |
5楼2013-12-02 00:07:13
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
6楼2013-12-02 10:54:26
40