| 查看: 1988 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料专业344求调剂
已经有19人回复
-
085404 298分求调剂
已经有11人回复
-
296求调剂
已经有7人回复
-
初试261
已经有6人回复
-
284求调剂
已经有9人回复
-
材料085601调剂
已经有25人回复
-
272分材料子求调剂
已经有23人回复
-
22408调剂求助
已经有7人回复
-
调剂
已经有17人回复
-
材料类284调剂
已经有43人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 专业: 理论和计算化学
2楼2012-12-26 23:57:16
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
三星20103楼2012-12-28 00:49:29
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
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
6楼2012-12-28 13:39:03
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
三星2010: 金币+30, ★★★很有帮助, 百度中也有个账号是dbb627，莫非就是阁下，真乃高人也，你的代码我看了一阵才明白，不过还是明白了，非常感谢！ 2012-12-28 19:48:07
三星2010: 金币+30, ★★★很有帮助, 百度中也有个账号是dbb627，莫非就是阁下，真乃高人也，你的代码我看了一阵才明白，不过还是明白了，非常感谢！ 2012-12-28 19:48:07
|
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=@(a,h)lamda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h))-1/beta/sin(lamda2*h); x=0:5:890; h=zeros(numel(x),1); for i=1:numel(x) opt=optimset('Display','off'); a=x(i); h(i)=fsolve(@(h)f(a,h),135,opt); end plot(x,h) |
7楼2012-12-28 14:39:01
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
8楼2012-12-28 20:10:29
★
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2012-12-30 01:55:12
jjdg: 编辑内容 2012-12-30 01:55
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2012-12-30 01:55:12
jjdg: 编辑内容 2012-12-30 01:55
|
本帖内容被屏蔽 |
9楼2012-12-28 20:47:05
三星2010
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1394.5
- 散金: 70
- 红花: 1
- 帖子: 197
- 在线: 95.4小时
- 虫号: 1210684
- 注册: 2011-02-23
- 性别: GG
- 专业: 力学中的基本问题和方法
10楼2012-12-29 09:43:57
