24小时热门版块排行榜    

查看: 353  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 竿竿儿 的主题更新
信息提示
请填处理意见