caizi2008 发表于 2007-7-13 17:06

请大家帮忙,谢谢

最近做关于吸声材料含粒子算等效参数的问题,可是根据表达式编程序求解怎么也解不出来,请大家帮小弟看一下程序并给予指点,困扰了我几天的问题了.谢谢大家,顺便说名下:我要求解的可能都是复数量;
      程序原代码:
    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]
查看完整版本: 请大家帮忙,谢谢