求解多元积分

shuddkk

clc;%总润滑方程求解
, |& `* w' z0 U8 m! F7 X4 rclear;
close all
4 @6 I6 R0 k: b2 Dsyms a b c d e x L O P U t o p
K=5;%常量
h=2;%常量 多孔储液介质厚度
D=0.012;%半径
H=0.2;%压缩后的高度
" m- o$ p* `" AL=@(x)0.2.*x;
% % % % 计算与孔隙率相关的brinkman数
c=@(t)((t.^3.*D.^2)./(36.*(1-t).^2.*K)).^2;%渗透率求解
t1=0.0033;%孔隙率
! u/ q+ z3 V1 ^6 Z* e0 wa=L./c(t1);%a为brinkman数
7 w+ F% R$ y# }8 {b=H./c(t1);
1 N, s8 d3 u1 h% % % % 计算方程中的f,控制此处f不与x一起变化
f=(cos(h.*b.*H)-1)/(b.*sin(h.*b.*H));
% % % % 计算方程中常数余项C,此处x只与a有关
o=@(x)(a(x).^2.*(H-f)./(2.*f-H));
" ~3 v  ~7 x+ H0 L, W  Ip=@(x)(a(x).^2./(2.*f-H));
O=int(o(x),0,1);
( K9 j) {% f" O; p$ }/ n& {P=int(p,0,1);
4 u1 x( P& {, i% L$ fc=O./P;
%%%%计算润滑升力
d=@(x)(1-x).*L./c(t1).^2.*(H-f+c)./(2.*f-H);
+ `! W" y# D# Z5 @9 G! M+ Y5 MU=int(d,0,1);
+ X4 D0 n, A. V' B: l- h我的代码一直在报错，请求高手帮忙看看哪里出了问题。
错误提示：未定义与 'function_handle' 类型的输入参数相对应的函数 'rdivide'。

出错 runhuashengli (line 13)
, B/ p: Y; C3 F" ~. ]) V& za=L./c(t1);%a为brinkman数

llbnmo

错的地方太多，你需要先搞清楚MATLAB匿名函数的定义方法和调用的语法规则
0 |3 `6 N+ t+ H: b& E

kekek

前面既然使用了匿名函数 L=@(x)0.2.*x; 那么后面除非是使用函数句柄时可以只写L，当你需要使用含有某个自变量的式子的时候，就需要写成 L(x)。
, q+ n7 `% I6 R2 v, S" t( j/ W( {" I
: B6 m0 L# b' X- M! T9 s相应地，a的赋值右侧包含L函数，所以a也需要写成函数形式，否则后面对o和p的计算，你写的这个a(x)就没法算了。

再者，避免变量与函数重复命名，你既写了匿名函数
0 U  C2 t: t$ K$ m8 B& F5 t( Kc=@(t)((t.^3.*D.^2)./(36.*(1-t).^2.*K)).^2;%渗透率求解
后面却又写了
c=O./P; %这算出来要么是个数值要么是个符号表达式

可是MATLAB系统里面只能有一个被命名为c的量，它不能同时表示这两者，这直接导致你最后对d的算式中同时出现的
6 v+ p! @! `! B3 E" a: [# q: [d = @( x )( 1 - x ) .* L( x ) ./ c( t1 ).^2 .* ( H - f + c )./( 2 .* f - H );
会出错
7 p$ n6 g9 s3 C" W' U6 f4 n+ K" P

kekek

syms x
K = 5;%常量
h = 2;%常量 多孔储液介质厚度
: {' B0 L+ }' f0 d+ y5 }D = 0.012;%半径
H = 0.2;%压缩后的高度
- q( m7 j$ P# r+ S3 g; U6 L. q9 g" O+ AL = @( x ) 0.2.*x;
% % % % 计算与孔隙率相关的brinkman数
c = @( t )( ( t.^3.*D.^2 )./( 36.*( 1-t ).^2.*K ) ).^2;%渗透率求解
t1 = 0.0033;%孔隙率
* C0 E9 t& f, ia =@( x ) L( x )./c( t1 );%a为brinkman数
b = H./c( t1 );
4 u, v( k3 R. j! s5 z; d; [% % % % 计算方程中的f,控制此处f不与x一起变化
  k8 s$ R8 k& l  R' Df = ( cos( h.*b.*H )-1 )/( b.*sin( h.*b.*H ) );
4 Z  W% H& s8 g' v$ Y0 B. d) R# S% % % % 计算方程中常数余项C,此处x只与a有关
o = @( x )( a( x ).^2.*( H-f )./( 2.*f-H ) );
p = @( x )( a( x ).^2./( 2.*f-H ) );
O = int( o( x ), 0, 1 );
P = int( p( x ), 0, 1 );