求解多元积分

shuddkk

clc;%总润滑方程求解
clear;
( C3 ]1 I, H3 ]& j  Rclose all
syms a b c d e x L O P U t o p
" F3 ~+ l# Y' e  F9 t$ xK=5;%常量
h=2;%常量 多孔储液介质厚度
D=0.012;%半径
H=0.2;%压缩后的高度
; S5 ?: q9 }/ f. \. OL=@(x)0.2.*x;
5 X# f9 i2 ^2 W( h! k2 l" z" S1 M% % % % 计算与孔隙率相关的brinkman数
& h* s" e% U( h+ ?) fc=@(t)((t.^3.*D.^2)./(36.*(1-t).^2.*K)).^2;%渗透率求解
; ~( J  E! E. T% ^* ^2 at1=0.0033;%孔隙率
6 K" K5 @8 N% j& }0 Wa=L./c(t1);%a为brinkman数
b=H./c(t1);
% % % % 计算方程中的f,控制此处f不与x一起变化
9 ]& ]" q" y: Tf=(cos(h.*b.*H)-1)/(b.*sin(h.*b.*H));
% % % % 计算方程中常数余项C,此处x只与a有关
1 h. ]! X0 r  |# b; H3 @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,0,1);
0 F" v) J+ K1 P$ |" n7 U& `: @: Kc=O./P;
%%%%计算润滑升力
d=@(x)(1-x).*L./c(t1).^2.*(H-f+c)./(2.*f-H);
U=int(d,0,1);
" C' g, n! U% h' y6 r6 p我的代码一直在报错，请求高手帮忙看看哪里出了问题。
错误提示：未定义与 'function_handle' 类型的输入参数相对应的函数 'rdivide'。
4 ~( a2 i0 p. y2 P7 A
/ _; x9 }- `) d4 [. O0 ]出错 runhuashengli (line 13)
* W+ G! B3 z$ Ca=L./c(t1);%a为brinkman数

llbnmo

错的地方太多，你需要先搞清楚MATLAB匿名函数的定义方法和调用的语法规则
  G# q3 s# G, y& k. k. p( D1 l" n

kekek

前面既然使用了匿名函数 L=@(x)0.2.*x; 那么后面除非是使用函数句柄时可以只写L，当你需要使用含有某个自变量的式子的时候，就需要写成 L(x)。

相应地，a的赋值右侧包含L函数，所以a也需要写成函数形式，否则后面对o和p的计算，你写的这个a(x)就没法算了。
9 u$ Z& u* a" a! r
再者，避免变量与函数重复命名，你既写了匿名函数
c=@(t)((t.^3.*D.^2)./(36.*(1-t).^2.*K)).^2;%渗透率求解
后面却又写了
% y3 ?. q) j  m" V6 Z* Qc=O./P; %这算出来要么是个数值要么是个符号表达式

可是MATLAB系统里面只能有一个被命名为c的量，它不能同时表示这两者，这直接导致你最后对d的算式中同时出现的
& N3 K/ c! h/ S+ gd = @( x )( 1 - x ) .* L( x ) ./ c( t1 ).^2 .* ( H - f + c )./( 2 .* f - H );
会出错

kekek

syms x
  D( E9 X2 n( i+ T, gK = 5;%常量
& @5 ^' K, L' Z! _9 s! I- ^h = 2;%常量 多孔储液介质厚度
4 D3 |' Z2 s3 U% D/ k) Z/ R# hD = 0.012;%半径
7 ?0 O5 d3 @1 X) Z! ~( IH = 0.2;%压缩后的高度
7 x: u6 C, `& {! U7 {" G/ ML = @( x ) 0.2.*x;
6 i! O; \$ f- P9 X+ q+ \% ?0 o% % % % 计算与孔隙率相关的brinkman数
c = @( t )( ( t.^3.*D.^2 )./( 36.*( 1-t ).^2.*K ) ).^2;%渗透率求解
t1 = 0.0033;%孔隙率
a =@( x ) L( x )./c( t1 );%a为brinkman数
7 Y1 Q; o, b- W6 `. ]0 R; ?b = H./c( t1 );
% % % % 计算方程中的f,控制此处f不与x一起变化
f = ( cos( h.*b.*H )-1 )/( b.*sin( h.*b.*H ) );
3 L; l" [5 }% v% % % % 计算方程中常数余项C,此处x只与a有关
o = @( x )( a( x ).^2.*( H-f )./( 2.*f-H ) );
% F: l. Y: x, [" {6 @p = @( x )( a( x ).^2./( 2.*f-H ) );
7 G( O: N2 a9 k% ?O = int( o( x ), 0, 1 );
2 ^, a$ q5 V) t! x& V( SP = int( p( x ), 0, 1 );