Question

matlab符号函数画图

求助~~~符号函数画图问题
我想画出一个t(R)的二维函数图~~~R的区间是（0,50^（-9）），这是我编的程序~~~但是结果老不对~~请教各位大侠！！
syms R n1 n2 a h r k u1 u2 e1 e2
f=6*R*(n1+h*(n2-n1)/(2*R))
g=2*a*R-2*r*h-k/R-(h^1.5)*(R^0.5)/(0.75*((1-u1^2)/e1+(1-u2^2)/e2))
t=int(f/g,h,0,2*R)
n1=10^(-3)
n2=2240
a=0.08*1.38*300*10^(-5)
k=20*1.38*300*10^(-23)
r=0.005*1.38*300*10^(-5)
e1=3000
e2=1.8*10^9
u1=0.38
u2=0.3
ezplot(t,0,50*10^(-9))

Answer 1

程序可以更改如下：但是执行结果表明找不到解析的结果，说明你的积分函数太复杂。以下满足语法规则但是找不到满意的结果。
% global n1 n2 a r k u1 u2 e1 e2 %这句可要可不要
n1=10^(-3);
n2=2240;
a=0.08*1.38*300*10^(-5);
k=20*1.38*300*10^(-23);
r=0.005*1.38*300*10^(-5);
e1=3000;
e2=1.8*10^9;
u1=0.38;
u2=0.3;

syms R h
f=6*R*(n1+h*(n2-n1)/(2*R));
g=2*a*R-2*r*h-k/R-(h^1.5)*(R^0.5)/(0.75*((1-u1^2)/e1+(1-u2^2)/e2));
t=int(f/g,h,0,2*R)

% ezplot(t,0,50*10^(-9))

建议尝试使用数值解。至于如何数值化需要看你的参数如何取值。

Answer 2

h没有定义

符号替换subs