24小时热门版块排行榜    

查看: 450  |  回复: 2

竿竿儿

新虫 (小有名气)

[求助] 急!毕业论文需要!求助:Undefined function or variable 'Etheta_34'. 已有1人参与

程序如下: 运行之后出现:Undefined function or variable 'Etheta_34'. 但是就算是我把Etheta_34赋值为0,都是同样的问题。求大神解决一下!毕业论文需要啊!急!
%2-Ca 3-Cl 4-SO4
a1=2;
a11=1.4;a22=12;
m2=(0.6435+30.4521)*0.001;m3=2*0.6435;m4=30.4521;
b=1.2;
A_fy=0.392;
beta0_23=0.3159;
beta1_23=1.614;
beta0_24=0.2;
beta1_24=3.1973;
beta2_24=-54.24;
Cfy_23=-0.00034;
Cfy_24=0;
pusy_234=-0.018;
theta_34=0.02;
I=0.5*(4*m2+m3+4*m4);
x34=6*(-1)*(-2)*A_fy*I^0.5;
x33=6*(-1)*(-1)*A_fy*I^0.5;
x44=6*(-2)*(-2)*A_fy*I^0.5;
Z=2*m2;
f_r=-A_fy*(I^0.5/(1+b*I^0.5)+2*log(1+b*I^0.5)/b);
g1=2*(a1*I^0.5)^(-2)*(1-(1+a1*I^0.5)*exp(-a1*I^0.5));
dg1=-2*(a1*I^0.5)^(-2)*(1-(1+a1*I^0.5+0.5*(a1*I^0.5)^2)*exp(-a1*I^0.5));
g11=2*(a11*I^0.5)^(-2)*(1-(1+a11*I^0.5)*exp(-a11*I^0.5));
dg11=-2*(a11*I^0.5)^(-2)*(1-(1+a11*I^0.5+0.5*(a11*I^0.5)^2)*exp(-a11*I^0.5));
g22=2*(a22*I^0.5)^(-2)*(1-(1+a22*I^0.5)*exp(-a22*I^0.5));
dg22=-2*(a22*I^0.5)^(-2)*(1-(1+a22*I^0.5+0.5*(a22*I^0.5)^2)*exp(-a22*I^0.5));
B_23=beta0_23+beta1_23*g1;
B_24=beta0_24+beta1_24*g11+beta2_24*g22;
C_23=Cfy_23/(2*2^0.5);
C_24=0;
dB_23=beta1_23*dg1/I;
dB_24=(beta1_24*dg11+beta2_24*dg22)/I;
y1=m3*(B_23+Z*C_23+theta_34)+m4*(B_24+Z*C_24);
y2=m2*(B_24+Z*C_24);
y3=m2*m3*(4*dB_23+2*C_23+0.5*pusy_234)+m2*m4*4*dB_24;
y4=0.5*m3*m4*(0.5*pusy_234+4*Etheta_34);
J(3,4)=x34*(4+(4.581*(x34^(-0.7237)))*exp(-0.012*x34^0.528))^(-1);
J(3,3)=x33*(4+(4.581*(x33^(-0.7237)))*exp(-0.012*x33^0.528))^(-1);
J(4,4)=x44*(4+(4.581*(x44^(-0.7237)))*exp(-0.012*x44^0.528))^(-1);
Etheta_34=(0.5/I)*(J(3,4)-0.5*J(3,3)-0.5*J(4,4));
dJ(3,4)=(4+4.581*x34^(-0.7237)*exp(-0.012*x34^0.528))^(-1)+(4+4.581*x34^(-0.7237)*exp(-0.012*x34^0.528))^(-2)*(4.581*x34*exp(-0.012*x34^0.528)*(0.7237*x34^(-4.581-1)+0.012*0.528*x34^(0.528-1)*x34^(-0.7237)));
dJ(3,3)=(4+4.581*x33^(-0.7237)*exp(-0.012*x33^0.528))^(-1)+(4+4.581*x33^(-0.7237)*exp(-0.012*x33^0.528))^(-2)*(4.581*x33*exp(-0.012*x33^0.528)*(0.7237*x33^(-4.581-1)+0.012*0.528*x33^(0.528-1)*x33^(-0.7237)));
dJ(4,4)=(4+4.581*x44^(-0.7237)*exp(-0.012*x44^0.528))^(-1)+(4+4.581*x44^(-0.7237)*exp(-0.012*x44^0.528))^(-2)*(4.581*x44*exp(-0.012*x44^0.528)*(0.7237*x44^(-4.581-1)+0.012*0.528*x44^(0.528-1)*x44^(-0.7237)));
dtheta_34=-(Etheta_34/I)+(2/(8*I^2))*(x34*dJ(3,4)-x33*dJ(3,3)/2-x44*dJ(4,4)/2);
y=4*f_r+y1+y2+y3+y4;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

mylifeljy

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
竿竿儿: 金币+10, ★★★很有帮助 2015-05-02 09:54:04
本帖内容被屏蔽

2楼2015-05-02 08:40:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竿竿儿

新虫 (小有名气)

引用回帖:
2楼: Originally posted by mylifeljy at 2015-05-02 08:40:00
lz,报错是因为你定义y4前没有定义Etheta_34!将J和Etheta_34的定义放到y4前面即可,修改后程序如下:
clc;  clear all; close all;
a1=2;
a11=1.4;a22=12;
m2=(0.6435+30.4521)*0.001;m3=2*0.6435;m4=30.4521;
...

谢谢!这样就调出来了,但是我以前的程序也是这么编的啊,都没出问题呀
3楼2015-05-02 09:53:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 竿竿儿 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见