查看: 4196|回复: 2|关注: 0
MATLAB怎么求解有贝塞尔函数的问题，求高手帮帮忙，谢谢
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
[code]%%声波声速、声衰减计算公式程序%%%%%
syms kb u w n p q m kesi x
coef={kb,u,w,n,p,q,m,kesi};
%u为骨架剪切模量,kb为骨架体变模量,p=(1-n)ps固相多孔隙介质的等效密度,q=npf流体的等效密度,w圆
频率,kr颗粒体变模量,kf流体体变模量,nta粘滞系数,kao渗透率,alfa弯曲度,a孔隙尺寸,n孔隙度%%；
kr=3.6e10;
kf=2.25e9;
nta=0.001;
kao=1.0e-10;
alfa=1.25;
D=kr*(1+(kr/kf-1)*n);
H=(kr-kb)^2/(D-kb)+kb+4*u/3;
C=kr*(kr-kb)/(D-kb);
M=kr^2/(D-kb);
m=1.25*1023/n;
for f=1:1200
kesi=a*(w*1023/nta).^1/2;
T=besselj(1,kesi)./besselj(0,kesi);
F=kesi.*T/4/(1-2*T/(i.*kesi));
A=C^2-H*M;
B=m*H*w^2+1999.2*w^2*M-i*H*w*F*nta/kao-2*C*1023*w^2;
C=-1999.2*m*w^4+i*1999.2*F*nta*w^3/kao+1023^2*w^4;
equa=A*x^4+B*x^2+C;
g=solve(equa,'x');
lr=subs(g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n**1023/n,5.0e-5*
(2*pi*f*)^1/2});
E=real(lr);
vp=[1 1 1 1]'./E;
li=subs(-w*g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n**1023/n,5.0e
-5*(2*pi*f*)^1/2});
F1=imag(li);
v=[vp(1) 0];
& & a=[ap(1) 0];
& & for j=2:4
& && &&&if&&vp(j)&vp(1)
& && && && &v(1)=vp(j);
& && && && &a(1)=ap(j);
& && &&&end
& & v1(f)=v(1);
& & a1(f)=a(1);
subplot(1,2,1);plot(v1,'-','linewidth',2);xlabel('f/Hz');ylabel('Vp/m/s');title('沉积物声速-
频率计算曲线');
subplot(1,2,2);plot(a1,'-','linewidth',2);xlabel('f/Hz');ylabel('ap/dB/m');title('沉积物声衰
减-频率计算曲线');
Function 'gt' is not implemented for MuPAD symbolic objects.[
不知道程序怎么弄了，有两个问题1.引入了贝塞尔函数，2.有复数解情况/code]
哪位高手帮我瞧瞧，非常感谢！
<h1 style="color:# 麦片财富积分
关注者： 3
clc
clear all
syms kb u w n p q m kesi x
coef={kb,u,w,n,p,q,m,kesi};
%u为骨架剪切模量,kb为骨架体变模量,p=(1-n)ps固相多孔隙介质的等效密度,q=npf流体的等效密度,w圆
%频率,kr颗粒体变模量,kf流体体变模量,nta粘滞系数,kao渗透率,alfa弯曲度,a孔隙尺寸,n孔隙度%%；
a=5e-5;
kr=3.6e10;
kf=2.25e9;
nta=0.001;
kao=1.0e-10;
alfa=1.25;
D=kr*(1+(kr/kf-1)*n);
H=(kr-kb)^2/(D-kb)+kb+4*u/3;
C=kr*(kr-kb)/(D-kb);
M=kr^2/(D-kb);
m=1.25*1023/n;
for f=1:1200
tt(f)=f;
w=2*pi*f;
kesi=a*(w*1023/nta).^1/2;
T=besselj(1,kesi)./besselj(0,kesi);
F=kesi.*T/4/(1-2*T/(i.*kesi));
A=C^2-H*M;
B=m*H*w^2+1999.2*w^2*M-i*H*w*F*nta/kao-2*C*1023*w^2;
CC=-1999.2*m*w^4+i*1999.2*F*nta*w^3/kao+1023^2*w^4; % 此处改C为CC，因与上面定义的参数重合
g=solve('A*x^4+B*x^2+CC',x);& & % 此处g有4个解，你取哪个？已经将equa=A*x^4+B*x^2+C;g=solve(equa,'x');改为现在形式
n=0.5;
lr=subs(g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n*1023,...
& & 1.25*1023/n,5.0e-5*(2*pi*f*)^1/2});
E=real(lr);
vp=[1;1;1;1]./E; % vp里面含有A,B,CC，没有给出代入计算值，这直接导致下面的for循环由于含有字符而无法比较
li=subs(-w*g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n*1023,...
& & 1.25*1023/n,5.0e-5*(2*pi*f*)^1/2});
F1=imag(li);
ap=F1;
v=[vp(1) 0];
& & a=[ap(1) 0]
& & for j=2:4
& && &&&if&&vp(j)&vp(1)
& && && && &v(1)=vp(j)
& && && && &a(1)=ap(j);
& && &&&end
& & end
& & v1(f)=v(1);
& & a1(f)=a(1)
& & end
subplot(1,2,1);plot(v1,'-','linewidth',2);xlabel('f/Hz');ylabel('Vp/m/s');
title('沉积物声速频率计算曲线');
subplot(1,2,2);plot(a1,'-','linewidth',2);xlabel('f/Hz');ylabel('ap/dB/m');
title('沉积物声衰减-频率计算曲线');
复制代码遗留的问题望LZ自己解决
[ 本帖最后由 xuwq909412 于
22:22 编辑 ]
<h1 style="color:# 麦片财富积分
(104 KB, 下载次数: 57)
22:48 上传
点击文件名下载附件
%%声波声速、声衰减计算公式程序%%%%%
syms kb u w n p q m kesi x
coef={kb,u,w,n,p,q,m,kesi};
%u为骨架剪切模量,kb为骨架体变模量,p=(1-n)ps固相多孔隙介质的等效密度,q=npf流体的等效密度,w圆频率,kr颗粒体变模量,kf流体体变模量,nta粘滞系数,kao渗透率,alfa弯曲度,a孔隙尺寸,n孔隙度%%；
a=5e-5;
kr=3.6e10;
kf=2.25e9;
nta=0.001;
kao=1.0e-10;
alfa=1.25;
D=kr*(1+(kr/kf-1)*n);
H=(kr-kb)^2/(D-kb)+kb+4*u/3;
C=kr*(kr-kb)/(D-kb);
M=kr^2/(D-kb);
m=1.25*1023/n;
for f=1:120
tt(f)=f;
w=2*pi*f;
kesi=a*(w*1023/nta)^1/2;
T=besselj(1,kesi)/besselj(0,kesi); 这里引入了一个贝塞尔函数，对它不懂，不知道怎么表达在f循环下矩阵的形式
F=kesi.*T./4./(1-2*T./(i.*kesi));
A=C.^2-H.*M;
B=m*H*w^2+1999.2*w^2*M-i*H*w*F*nta/kao-2*C*1023*w^2;
CC=-1999.2*m*w^4+i*1999.2*F*nta*w^3/kao+1023^2*w^4;这里A,B,CC表达也应该是矩阵形式吧，有跟频率相关的角频率w
g=solve('A*x^4+B*x^2+CC','x')
lr=subs(g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n**1023/n,5.0e-5*(2*pi*f*)^1/2});
E=real(lr);
vp=[1 1 1 1]'./E;这里有四个解，我取正解的最大值的，后面的比较由于A,B,CC都是含有f的变量，所以毛病很多
li=subs(-w*g,coef,{4.4e7-i*2.0e6,2.6e7-i*1.25e6,2*pi*f,n,(1-n)*2650,n**1023/n,5.0e-5*(2*pi*f*)^1/2});
F1=imag(li);
ap=F1;
v=[vp(1) 0];
a=[ap(1) 0];
for j=2:4
if vp(j)&vp(1)
v(1)=vp(j);
a(1)=ap(j);
end
end
v1(f)=v(1);
a1(f)=a(1);
end
subplot(1,2,1);plot(eval(v1),'-','linewidth',2);xlabel('f/HZ');ylabel('Vp/m/s');title('沉积物声速-频率计算曲线');
subplot(1,2,2);plot(eval(a1),'-','linewidth',2);xlabel('f/HZ');ylabel('ap/dB/m');title('沉积物声衰减-频率计算曲线');
复制代码
求解这样的四次解，分为实部和频率图和虚部和频率图，比较难啊，阁下能不能帮帮忙啊，先谢谢了
站长推荐 /2
Powered by已收藏本页面
帮帮忙，解一下
最讨厌老坟
妈逼，叫他俩开房好了
男人苦的是连襟，女人是死者的女儿，男人是女人的姨夫
蒙逼了
[流鼻血]
这不重要，重要的是他们开房了没有，，？
在我家这面男人和女婿是一担挑 女人是大姨子或者小姨子
这男的大女的一个辈分，这女的是男人的晚辈，极可能没有关系
女的是他小姨子
这是个多选题
楼下什么结论就想上神评？这个只有一种可能，这个男的俩姨娶了自己丈母娘！
姐夫和小姨子
是姐夫和小姨子
绝逼是丈母娘和女婿的关系
你妈和隔壁老王
丈母娘和女婿
母子关系，女儿没来！
女人应该是男人的小姨子或者大姑姐！
死人和活人的关系
丈母娘和女婿啊
回复 47楼：你能去给你老婆的前任上坟啊，佩服佩服，你感情真丰富。
回复 137楼：的丈人
姨夫和外甥女，姨夫哭的是自己的连襟，就是自己丈人的另一女婿，女的哭的是自己的爹，也是自己女婿（老公）
姐夫和小姨子
姨夫跟外甥女的关系
姐夫小姨子，或者妹夫和大姨姐
好恐怖，哭的是自己墓呢
回复 130楼：的亲爹，岂不就是自己的老公。呵呵。
回复 112楼：不错不错，但绝不能叫混了，这个死人是自己老丈人的女婿，也就是自己的连襟，或姐夫或妹夫。但女人却是指着自己的女儿叫的，自己女婿的老丈人，女儿
回复 41楼：这智商，绝对外星人
回复 9楼：应该是女人死了老公，然后再找了一个，那个男的和她一起去上坟，女人哭时因为想亡夫，男的哭是因为后悔娶这个女人了，为什么后悔呢，因为，靠实在编不下去了
男的是女子的姐夫，或者妹夫
那个男的应该是丈人的女婿的亲生兄弟类吧
女的不能上坟
回复 115楼：对
回复 113楼：他大姨夫跟他妈去给他爸上坟
姐夫和小姨子
姨夫跟外甥女
男人是女人的姐或妹的老公
另一个女婿和丈母娘
姐夫和小姨子
丈母娘和姑爷
男人就是丈人的女婿，所以这个丈人有两女儿；女人女婿的丈人就是她自己老公；结论女人是男人的小姨子，并且都是中年人。带我上神评～
又TM是姐夫和小姨子干的
^_^要么妻姐或者妻妹
哭的神贴的
外婆外孙关系
我诧异的是那男的：为什么会哭他丈人的女婿，这还有点信息量……
回复 100楼：女的是男的的丈母娘，也是死者的老婆，这么简单都能搞错。
丈母娘和女婿关系。
姐夫和小姨子
回复 98楼：正解
男的和死者是连襟，女的是死者的女儿，这两个人关系是姨夫个妻子的外甥女吧！
哪个男的会给自己丈人的女婿哭坟？除非是自己亲弟或亲哥吧
男女夫妻关系，女的哭的是前夫
女婿和丈母娘
女人是男的丈母娘
乱！！！！！
姐夫（妹夫）和小姨子的关系 就这么简单！
认不到！喊老表！
这么简单，没关系啊，各哭各的
是楼主和楼主的丈母娘吧？？
[笑着流泪]
[笑着流泪]
姐夫跟小姨子
姐夫与小姨子
姐夫和小姨子
男女关系，太简单了，整个五岁的问题来
贱人就是矫情。。。。。
男的是楼主的姐夫（妹夫），女的是楼主的老婆。这就是姐夫和小姨子（大姨姐）的关系。
跟我有啥关系
卧槽，跟我有啥关系！
你家上个坟都tm事多
丈母娘和女婿！准的
Sb，你和你妈去哭你爹
外甥女和姨夫去她爸坟上哭了
丈母娘和女婿
两口子，再婚关系，而且这女人嫁的是前夫的亲人
正在拿笔算的。请举手
亲兄妹，为亲爸上坟
姐夫和小姨子的关系！
姑爷和丈母娘
姐夫和小姨子
好像是舅舅和外甥女

[doge]
[doge]
[doge]亲戚关系
姐夫和小姨子
回复 9楼：瞎扯
二百年前的烂梗
[doge]
莫非是楼主！
死了的那个人，是那个女人的前夫
[嘿嘿]
[嘿嘿]
没有特定是同一个丈人或者同一个女婿吧？
假设坟子里是A，根据男人的描述，男人是A的连襟，换句话说，男人和A的妻子，都是同一个老爸。再说说女人，女人的女婿，女婿的丈人，就是女人和A是夫妻。所以结论是男人和女人是姐夫和小姨子。
女婿与丈母娘的故事
苍井空和安培
唉～竟然没有一个说对的！男女的关系是外公和外孙女
不假思索，男女关系！
女婿和丈母娘的关系呗
姐夫和小姨子
回复 34楼：他不只一个女婿不啊
好诡异！！丈人的女婿是他自己！女婿的丈人是他老公！男的是鬼！
男的是女的的姨夫
什么关系不重要！重要的是“这么巧！你也来嘚瑟”
[大笑]
[大笑]
[大笑]
大姨姐和妹夫
大姨子和妹婿，或者小姨子和姐夫
男人哭连襟，女人哭老公。
丈人是个什么鬼
哇，好老的，我96年听过，又遇到。
女婿丈母娘
狗男女，哭完不忘野战
小姨子与姐夫，或者大姨子与妹夫！
回复 18楼：正解
回复13楼：
[大哭]
[嘿嘿]
女婿丈母娘
女的和死者是夫妻关系，男的是女的
姐夫与小姨子
姐夫跟小姨子，死者是女人的丈夫，男人的连襟
大嫂和小叔
大姨子和妹夫；姐夫和小姨子；再婚的丈夫和未亡人
丈母娘，和女婿
男人是死者的姨夫，女人的老公，由此得出结论，女人是男人的姨姐或小姨子
回复5楼：不是了
纯洁的男女关系
回复 3楼：男的在哭自己的女婿，女的在哭自己的老公
回复 1楼：我也不知道啊
后才能发表评论
你可能喜欢的图片笑话
最受欢迎的爆笑笑话查看: 7256|回复: 8
最后登录在线时间2 小时CPS积分58 分注册时间阅读权限20帖子精华0UID475083
安防小学生, 积分 58, 距离下一级还需 42 积分
请问一下通过走云服务器产生的延迟问题你有没有好的解决方案啊谢谢，看你很有经验的感觉啊。或者给点你的看法建议什么的谢谢，镜头是海康威视的DS-2CD4026Fwd,200万的像素，延迟比较厉害。
最后登录在线时间43 小时CPS积分516 分注册时间阅读权限90帖子精华0UID474713
安防精英, 积分 516, 距离下一级还需 84 积分
试一下千兆交换机试试
北京安防一级资质,安防资质,安防公司,www.afzz.ne t 有需要咨询的朋友 联系我 QQ
欢迎加入安防交流群：
最后登录在线时间49 小时CPS积分135 分注册时间阅读权限30帖子精华0UID431436
安防中学生, 积分 135, 距离下一级还需 65 积分
局域网都卡，还上云服务？尽量减少数据流，比如调节帧率，码流
人活着本来就没什么意义，但只有活下去，才能找到有趣的事物……就像你找到了这朵花， ...
最后登录在线时间63 小时CPS积分986 分注册时间阅读权限30帖子精华0UID475194
安防大师, 积分 986, 距离下一级还需 14 积分
手动改码流，网络核心部分改成千兆交换机。传输部分用光传输，或者用光纤收发器。收发器的成本比模块低很多！估计你的前端设备少不了
最后登录在线时间353 小时CPS积分2253 分注册时间阅读权限255帖子精华1UID195683
安防偶像, 积分 2253, 距离下一级还需 747 积分
调节帧率，码流
深圳盼宝科技有限公司，专业监控产品生产厂家。电话：8、566038
最后登录在线时间59 小时CPS积分2129 分注册时间阅读权限255帖子精华0UID473259
安防偶像, 积分 2129, 距离下一级还需 871 积分
深圳邦亿视安防系统有限公司
Byvision security systems Limited
www-byvisionss-com
主营:CCTV camera,HD-CAM,IPC,DVR KIT,NVR KIT,VDP
最后登录在线时间44 小时CPS积分374 分注册时间阅读权限70帖子精华0UID473243
安防工作者, 积分 374, 距离下一级还需 26 积分
如果是百兆交换机，两端都强制100M，全双工试试，千兆也一样。
最后登录在线时间747 小时CPS积分7901 分注册时间阅读权限255帖子精华0UID463155
热心朋友还是多呢，看看，学习中
物联网摄像机研发制造商,提供方案模组二次开发（SDK包）微信热线：
最后登录在线时间81 小时CPS积分1063 分注册时间阅读权限255帖子精华0UID426808
安防知名大师, 积分 1063, 距离下一级还需 437 积分
个人觉得你得咨询你们的 厂家，首先排除产品是否有无这方面的缺陷，贷款OK的情况下延时多久是正常，内网多久延时是正常的 。然后再优化其他（成本允许的前提下）设备。当然你也可以和楼上说的 试着改变码流及分辨率等信息、不过这些都是必须建立在你客户接受的前提之下。希望可以帮到你
Powered byRe: 请问这是什么问题？
我是这方面的白痴，只是因为游戏进不去，想解决下。大哥们帮帮忙 - 高端调试 :: 论坛求解！！！急用，高手们帮帮忙，万分感谢！！！ -oracle- TryCatch
>> Content
求解！！！急用，高手们帮帮忙，万分感谢！！！
装oracle11g&发行版2&的时候&先决条件检查&全部失败下面是详细信息：Checks:物理内存&-&此先决条件将测试系统物理内存总量是否至少为&922MB&(KB)。预期值?:?N/A实际值?:?N/A?错误列表:&?-?PRVF-7531&:&无法在节点&"HP-"&上执行物理内存检查&?-&Cause:?&无法在指示的节点上执行物理内存检查。&?-&Action:?&确保可以访问指定的节点并可以查看内存信息。&可用物理内存&-&此先决条件将测试系统可用物理内存是否至少为&50MB&(51200.0KB)。预期值?:?N/A实际值?:?N/A?错误列表:&?-?PRVF-7563&:&无法在节点&"HP-"&上执行可用内存检查&?-&Cause:?&无法在指示的节点上执行可用内存检查。&?-&Action:?&确保可以访问指定的节点并可以查看内存信息。交换空间大小&-&此先决条件将测试系统是否具有足够的总交换空间。预期值?:?N/A实际值?:?N/A?错误列表:&?-?PRVF-7574&:&无法在节点&"HP-"&上执行交换空间大小检查&?-&Cause:?&无法在指示的节点上执行交换空间检查。&?-&Action:?&确保可以访问指定的节点并可以查看交换空间信息。&?-?PRVF-7531&:&无法在节点&"HP-"&上执行物理内存检查&?-&Cause:?&无法在指示的节点上执行物理内存检查。&?-&Action:?&确保可以访问指定的节点并可以查看内存信息。空闲空间:&HP-:C:\DOCUME~1\ADMINI~1\LOCALS~1\Temp&-&此先决条件将测试文件系统是否具有足够的空闲空间。预期值?:?N/A实际值?:?N/A体系结构&-&此先决条件将测试系统是否具有认证的体系结构。预期值?:?N/A实际值?:?N/A?错误列表:&?-?PRVF-7536&:&无法在节点&"HP-"&上执行体系结构检查&?-&Cause:?&无法确定系统体系结构。&?-&Action:?&确保正在使用正确的软件包Environment&variable:&"PATH"&-&This&test&checks&whether&the&length&of&the&environment&variable&"PATH"&does&not&exceed&the&recommended&length.预期值?:?N/A实际值?:?N/A以上&先决条件检查全部失败了。求高手们帮帮忙怎么改才能通过？物理内存我现在知道可以忽略，其他的也可以忽略吗？会不会对程序的运行产生影响？[img=http://hi.csdn.net/space--do-album-picid-961072.html][/img]不知道图能不能看到，我空间里有截图，恳求高手们帮忙解决问题。解决问题者，100分&送上。万分感谢！
------Solutions------
打开C盘的默认共享，再重新检查即可！
------Solutions------
看下你服务中的&"TCP/IP&NetBIOS&Helper"&服务有没有启动
------Solutions------
引用&1&楼&funfenffun&的回复:打开C盘的默认共享，再重新检查即可！谢谢啦，其他的都通过了，还剩下最后一个问题，请帮忙，万分感谢！！Environment&variable:&"PATH"&-&This&test&checks&whether&the&length&of&the&environment&variable&"PATH"&does&not&exceed&the&recommended&length.预期值?:?N/A实际值?:?N/A
------Solutions------
引用&2&楼&opps_zhou&的回复:看下你服务中的&"TCP/IP&NetBIOS&Helper"&服务有没有启动服务启动了，还剩下一个问题没解决，请帮忙，万分感谢！！！
------Solutions------
环境变量Path的长度超长，右键点击我的电脑-&属性-&高级-&环境变量，找到Path，把一些不用的部分去掉就可以了
------Solutions------
引用&5&楼&funfenffun&的回复:环境变量Path的长度超长，右键点击我的电脑-&属性-&高级-&环境变量，找到Path，把一些不用的部分去掉就可以了谢谢了，已经解决啦，万分感谢啊
------Solutions------
别都删掉,要保留几个重要的。把原来的path值备份下来(一定要备份),换成:&%SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem再看看"用户变量“里头有没有path，如果有，将其删除。然后注销，重新登录，再安装，应该没有问题。
------Solutions------
引用&5&楼&funfenffun&的回复:环境变量Path的长度超长，右键点击我的电脑-&属性-&高级-&环境变量，找到Path，把一些不用的部分去掉就可以了安装完成后C盘默认共享能不能关闭呀？
------Solutions------
net&share&c$&/delete
------Solutions------
引用&9&楼&funfenffun&的回复:net&share&c$&/delete谢谢啦，万分感谢！
------Solutions------
是怎么解决的呀&我遇到同样问题&一直解决不了都疯了&帮帮我吧
------Solutions------
引用&11&楼&beauling&的回复:是怎么解决的呀&我遇到同样问题&一直解决不了都疯了&帮帮我吧1&楼&正解，网上查一下怎么打开C盘默认共享，打开后再检查就ok了，希望能帮到你}