请大家帮忙,谢谢
最近做关于吸声材料含粒子算等效参数的问题,可是根据表达式编程序求解怎么也解不出来,请大家帮小弟看一下程序并给予指点,困扰了我几天的问题了.谢谢大家,顺便说名下:我要求解的可能都是复数量;程序原代码:
clc;
clear all;
close all;
ra=0.5;
dertaa=0.001;
fI=0.005;
EM=6e8+j*0.45*6e8;
EI=5.4e8+j*0.32*5.4e8;
Ec=5e8+j*0.7*5e8;
xigmaM=0.45;
xigmaI=0.32;
xigmac=0.37;
lamdaM=EM*xigmaM/((1+xigmaM)*(1-2*xigmaM));
lamdaI=EI*xigmaI/((1+xigmaI)*(1-2*xigmaI));
lamdac=Ec*xigmac/((1+xigmac)*(1-2*xigmac));
miuM=EM/(2*(1+xigmaM));
miuI=EI/(2*(1+xigmaI));
miuc=Ec/(2*(1+xigmac));
KM=EM/(3*(1-2*xigmaM));
KI=EI/(3*(1-2*xigmaI));
Kc=Ec/(3*(1-2*xigmac));
alphaM=1/3*(1+xigmaM)/(1-xigmaM);%%%%%,可算
alphac=1/3*(1+xigmac)/(1-xigmac);%%%%%,可算
alphaI=1/3*(1+xigmaI)/(1-xigmaI);%%%%%,可算
beitaM=2/15*(4-5*xigmaM)/(1-xigmaM);%%%%%,可算
beitaI=2/15*(4-5*xigmaI)/(1-xigmaI);%%%%%,可算
beitac=2/15*(4-5*xigmac)/(1-xigmac);%%%%%,可算
faic=(Kc+alphac*(KI-Kc))/Kc; %%%%%,可算
fsaic=(miuc+beitac*(miuI-miuc))/miuc;%%%%%,可算
J2=3*dertaa/ra*(miuI-miuc)/miuc;
J3=3*dertaa/ra*(KI-Kc)/Kc;
J4=Kc*(1-alphac);
J5=miuc*(1-beitac);
J6=fI*(lamdaI-lamdaM);
J7=fI*2/3*(miuI-miuM);
J8=fI*3*dertaa/ra*(lamdac-lamdaM)*faic;
J9=fI*3*dertaa/ra*2/3*(miuc-miuM);
J10=fI*((miuI-miuM)+3*dertaa/ra*(miuc-miuM)*fsaic);
K1=J6+J7+J8+J9;
K2=J7+J9*fsaic;
K3=1+J3*alphac;
K4=KI-J3*J4;
K5=1+J2*beitac;
K6=miuI-J2*J5;
eq1='miue=Ee/(2*(1+xigmae))';
eq2='Ke=Ee/(3*(1-2*xigmae))';
eq3='lamdae=Ee*xigmae/((1+xigmae)*(1-2*xigmae))';
eq4='alphae=1/3*(1+xigmae)/(1-xigmae)';
eq5='beitae=2/15*(4-5*xigmae)/(1-xigmae)';
eq6='A=Ke/(Ke*(1-alphae)*K3+alphae*K4)';
eq7='B=miue/(miue*(1-beitae)*K5+beitae*K6)';
eq8='lamdae=lamdaM+K1*A-K2*B';
eq9='miue=miuM+B*J10';
=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,'Ee,xigmae,miue,Ke,beitae,lamdae,alphae,A,B');
页:
[1]