声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1330|回复: 0

[综合讨论] 求助,关于desolve函数

[复制链接]
发表于 2010-4-20 19:26 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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 编辑 ]
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 16:23 , Processed in 0.061776 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表