level 1
k = 2.5*10^-8;
T = 293.2;
tao = 0.05;
p = 1000;
R = 287;
pv = 3000;
ci = 18.2*10^-3;
cs = 1.44*10^-3;
NDSolve[{dd'[t] ==
2 k R T (ci - cs)/(p - pv +
8*tao/3/dd[t])*((dd[t]/2)^-1 + (\[Pi]*k*t)^-0.5),
dd[10^-10] == 10^-5}, dd, {t, 10^-10, 0.014}];
Plot[(dd[t]/10^-6) /. %, {t, 10^-10, 0.014}]
x = (2*k*zeta*cs/roug/R0/R0*t)^(1/2)
我现在画了一个dd[t]和t的图像，然后x是t的函数，现在想画d是纵坐标x是横坐标的图像。
2021年05月02日 08点05分
1
level 1
roug = 14.53;
zeta = 0.14;
有两个参数没给
2021年05月02日 08点05分
2
level 9
roug = 14.53;
zeta = 0.14;
k = 2.5*10^-8;
T = 293.2;
tao = 0.05;
p = 1000;
R = 287;
pv = 3000;
ci = 18.2*10^-3;
cs = 1.44*10^-3;
s = NDSolve[{dd'[t] ==
2 k R T (ci - cs)/(p - pv +
8*tao/3/dd[t])*((dd[t]/2)^-1 + (\[Pi]*k*t)^-0.5),
dd[10^-10] == 10^-5}, dd, {t, 10^-10, 0.014}];
Plot[(dd[t]/10^-6) /. %, {t, 10^-10, 0.014}]
x = (2*k*zeta*cs/roug/R0/R0*t)^(1/2)
你的R0没有给，我瞎猜了一个
R0 = 1
2021年05月02日 15点05分
3
ParametricPlot[{x, First[(dd[t]/10^-6) /. s]}, {t, 10^-10, 0.014}] First函数是用来去除{},不然会让ParametricPlot函数误解
2021年05月02日 15点05分
@🌚黑黄🌝 您好 我改成了下面的样子 不知道为什么 还是没图像
2021年05月06日 14点05分
@🌚黑黄🌝 我把你的原封不动都复制了，还是没画出图来
2021年05月08日 02点05分
level 1
clear all; clc;
roug = 14.53;
zeta = 0.14;
k = 2.5*10^-8;
T = 293.2;
tao = 0.05;
p = 1000;
R = 287;
pv = 3000;
ci = 18.2*10^-3;
cs = 1.44*10^-3;
R0 = 10^-5;
s = NDSolve[{dd'[t] ==
k *cs/roug *zeta*((dd[t])^-1 + (\[Pi]*k*t)^-0.5),
dd[10^-8] == 10^-5}, dd, {t, 10^-10, 0.010}];
ParametricPlot[{((2*k*zeta*cs/roug/R0/R0*t)^(1/2)), First[s/R0]}, {t,
10^-10, 0.010}]
2021年05月06日 14点05分
4
我照例要问一句，你们为什么觉得在代码开头写上 clear all; clc; 能起到你们想象中的作用？是你们自己想当然的还是从哪个教程里学的？
2021年06月05日 01点06分
level 4
ClearAll["Global`*"]
roug = 14.53;
zeta = 0.14;
k = 2.5*10^-8;
T = 293.2;
tao = 0.05;
p = 1000;
R = 287;
pv = 3000;
ci = 18.2*10^-3;
cs = 1.44*10^-3;
R0 = 10^-5;
sol = NDSolve[{dd'[t] ==
k*cs/roug*zeta*((dd[t])^-1 + (\[Pi]*k*t)^-0.5),
dd[10^-8] == 10^-5}, dd, {t, 10^-10, 0.010}];
x[t_] := (2*k*zeta*cs/roug/R0/R0*t)^(1/2);
dd = dd /. First@sol;
ParametricPlot[Evaluate[{x[t], dd[t]}], {t, 10^-10, 0.010},
AspectRatio -> 0.75]
2021年05月08日 02点05分
5