| 查看: 2351 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
三星2010木虫 (小有名气)
|
[求助]
请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
|
|||
|
主要工作就是求解隐函数方程，参数及方程如下： taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=250000; t1=5.12; b2=400; E2=10000; t2=20; L=1200; 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; 参数a和h满足一下关系l amda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h)-1/beta/sin(lamda2*h)=0 我能肯定的是a从0单调增大至892.806169，h的范围大致是[135,150]我想要的结果是将a的变化区间若干等分，比如1000等分，对每一个a去求对应的h并将这两个参数输出. 我自己编写的程序如下： for i=1:1:892.9/0.001 a(i)=i*0.001-0.001; error=1; for j=135/0.001:1:150/0.001 h1(j)=j*0.001; error1(j)=lamda2/lamda1*sinh(lamda1*(L-a(i)-h1(j)))- ... cot(lamda2*h1(j))*cosh(lamda1*(L-a(i)-h1(j)))- ... 1/beta/sin(lamda2*h1(j)); if error>abs(error1(j)) error=abs(error1(j)); h(i)=h1(j); end end if L end ltpu(i)=lamda2/lamda1*tanh(lamda1*(L-h(i)-a(i)))*cos(lamda2*h(i))+sin(lamda2*h(i)); ltpb(i)=lamda2/((lamda1)^2*delta1*b1*E1*t1)*(1+cos(lamda2*h(i))/(beta*cosh(lamda1*(L-a(i)-h(i))))); ltp(i)=ltpu(i)/ltpb(i); lts(i)=deltaf+ltp(i)*a(i)/(b1*E1*t1); end plot(lts,ltp,'b'); 程序运行后总提示h未定义，不知道是哪里出了问题，有没有更好的办法求解啊，我都被这个问题折磨晕了  ，请大家出手相助！[ Last edited by 三星2010 on 2012-12-26 at 16:36 ] |
回复此楼» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有10人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab编程中遇到的问题
已经有4人回复
- 我想请问一个matlab编程的问题，sin（v/2）=3*v/14/pi
已经有10人回复
- 请问如何用Matlab编程做这个图？谢谢大家！
已经有5人回复
- 求高人指点 MATLAB 编程求解 洛伦兹的最大李雅普诺夫指数程序！！！！
已经有5人回复
- 100BB急求matlab编程中的问题解决，如何修改这个程序呢？
已经有4人回复
- 【求助】vb与matlab混合编程时图形显示窗口中，怎样使图形比窗口小一些？
已经有3人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】请问matlab编程语言如何打包成为独立的应用程序啊？谢谢
已经有8人回复
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
三星20103楼2012-12-28 00:49:29
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
2楼2012-12-26 23:57:16
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
>> taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=250000; t1=5.12; b2=400; E2=10000; t2=20; L=1200; 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; syms a h f=lamda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h))-1/beta/sin(lamda2*h); ezplot(f,[0 900]) |
» 本帖已获得的红花（最新10朵）
4楼2012-12-28 09:05:17
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
5楼2012-12-28 13:29:32
