请高人帮我看一下Matlab一段程序：p12怎么改？我怎么调不出来。

clear%熵波耗散模型
clc
w=0:0.01:400*2*pi;%RD-170
mg=930;%燃气发生器质量流量
rm=54.5;%混合比
mx=rm/(1+rm)*mg;mf=1/(1+rm)*mg;%氧化剂和燃料质量流量
tc=0.006;%燃气停留时间（取0.004~0.008）
mt=tc*mg;%燃气发生器容腔气体质量
pg=53.5*10^6;%燃气发生器压力
kco=pg/mg;kcf=pg/mg;%相对于氧化剂秒流量、燃料秒流量的压力放大系数
Tg=760;%温度
tt=0.003;%燃烧时滞
a=340;%声速
l=1.5;%气体管路长度
f=-1.5;%燃烧产物温度与推进剂组元混合比的关系曲线的无量纲斜率
k=1.2;%比热比
Ma=0.15;
kt=0.8;%导热系数
Pea=a*l/kt;
aa=k*Ma;
%subplot(211)
p11=0.5.*(exp(i*w*l/a/(1-Ma))+exp(-i*w*l/a/(1+Ma))).*exp(-i*w*tt).*kco./(1+i*w*tc)*mx/pg+(-aa/2.*(exp(i*w*l/a/(1-Ma))-exp(-i*w*l/a/(1+Ma))))*mx/mg.*exp(-i*w*tt);
p12=f*(k.*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea).^2./((i*w*l/a+Ma.*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea)).^2-(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea).^2)/Pea*(0.5.*(exp(-i*w*l/a/(1+Ma))-exp(i*w*l/a/(1-Ma))).*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea)-(i*w+Ma.*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea)).*(0.5.*(exp(i*w*l/a/(1-Ma))+exp(-i*w*l/a/(1+Ma)))-exp(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea))));
p1=p11+p12
%plot(w/2/pi,abs(p1))
grid on
title('系统频率特性')
xlabel('频率/Hz')
ylabel('|\deltaP/\deltaG_x|')
subplot(212)
p2=0.5.*(exp(i*w*l/a/(1-Ma))+exp(-i*w*l/a/(1+Ma))).*exp(-i*w*tt).*kcf./(1+i*w*tc)*mf/pg+(-aa/2.*(exp(i*w*l/a/(1-Ma))-exp(-i*w*l/a/(1+Ma))))*mf/mg.*exp(-i*w*tt);
plot(w/2/pi,abs(p2))
grid on
title('系统频率特性')
xlabel('频率/Hz)')
ylabel('|\deltaP/\deltaG_\Gamma|')

答案：擦，不知道变量是干啥用的么，多写几个变量会死啊，你见谁写这么长的公式能不出错的？？
p12a = k*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea).^2./((i*w*l/a+Ma*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea)).^2-(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea).^2)/Pea;
p12b = 0.5*(exp(-i*w*l/a/(1+Ma))-exp(i*w*l/a/(1-Ma))).*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea)-(i*w+Ma.*(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea));
p12c = (0.5.*(exp(i*w*l/a/(1-Ma))+exp(-i*w*l/a/(1+Ma)))-exp(-i*w*l/a/Ma-(w*l/a).^2/Ma^3/Pea));
p12 = f*(p12a.*p12b.*p12c);
主要错误就是p12a后面少了一个点

上一个：求轮系结构的智能移动机器人轨迹跟踪matlab仿真程序，最好是用simulink做的,邮箱nba_86@qq.com，谢谢
下一个：matlab程序出错，Inner matrix dimensions must agree. 希望高手能帮忙修改，绝对给分！