|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本人的程序如下:
%uniform cylindrical tube terminated with a impedence
clc; clear all; close all;
ear_length=0.02; ear_radius=3.75e-3; ear_area=pi*ear_radius^2;
freq=4000; c=340; k=2*pi*freq/c; rho=1.41;
Zsp=402; Sd=12.44e-6; Zd=Zsp/Sd; omega=2*pi*freq;
syms z
pressure_z=dsolve('D2p=-(k)^2*p','p(0)=2','Dp(0)=0','z');
pressure_z_10=subs(pressure_z);
pressure_z_1=subs(pressure_z_10,z,0.003)
pressure_z_2=diff(pressure_z_10,'z');
pressure_z_3=subs(pressure_z_2,z,0.003);
pressure_z_4=dsolve('D2p=-(73.92)^2*p','Dp(0.003)=pressure_z_3+j*1.41*omega*pressure_z_1/(4.4179e-5*3.2315e7)','p(0.003)=pressure_z_1','z');
pressure_z_5=subs(pressure_z_4);
pressure_z_6=vpa(subs(pressure_z_5,z,0.003),5)
我想问一下:为什么两个的dsolve接下来的结果一个数值解,一个是符号解。
而且第二个dsolve后面的subs好像是不起作用的。求高手解答。
小弟在次拜谢。
[ 本帖最后由 ChaChing 于 2010-6-15 15:32 编辑 ] |
|