Z={(a*x+b*y)*(x*y)/[(0.5x+0.5y)^2]}^2)如何进行系统结构画图仿真（使用matlab)

nlm0402 发表于 2011-8-10 06:27
谢谢您的帮助，如果能帮忙写出MATLAB或vensim的系统结构图。
将非常感谢，倾我所有给予您报酬。

但是vensim貌似做不出三维图呀。

nkunku 发表于 2011-8-10 15:34
但是vensim貌似做不出三维图呀。

我实际上需要的系统仿真图，不要三维图。

tmdxyz 发表于 2011-8-8 16:57
这个问题Matlab我不会。但是如果用Mathematica来做，就很简单了。

其命令为（建议版本在7.0以上）：

如果能够滑动x、y就好了，就像设置a、b那样。

nlm0402 发表于 2011-8-10 15:52
如果能够滑动x、y就好了，就像设置a、b那样。

那就只能变通一下。如：
Animate[Plot3D[((a*(x + c) +
       b*(y + d))*((x + c)*(y + d))/(0.5*(x + c) +
         0.5*(y + d))^2)^2, {c, 0, 1}, {d, 0, 1},
  PlotRange -> Automatic], {a, 0, 1}, {b, 0, 1}, {x, 0, 1}, {y, 0, 1},
  AnimationRunning -> False]
不知是否合你的意？

nlm0402 发表于 2011-8-10 15:46
我实际上需要的系统仿真图，不要三维图。

附件是假设的模型ppp.mdl，其中设
a，b，x，y均为【0,1】之间的均匀分布随机数。即
a=RANDOM UNIFORM(0,1,12345)
b=RANDOM UNIFORM(0,1,12345)
x=RANDOM UNIFORM(0,1,12345)
y=RANDOM UNIFORM(0,1,12345)
z=={(a*x+b*y)*(x*y)/[(0.5x+0.5y)^2]}^2)

于是得到z的变化情况，如附图所示。

上述回答有些牵强，纯粹消遣，哈哈

clear
a=0.24;
b=0.54;
[X Y]=meshgrid(0:0.01:1, 0:0.01:1);
Z=((a*X + b*Y)*(X*Y)/(0.5*X + 0.5*Y)^2)^2;
mesh(X, Y, Z)

Animate[Plot[((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2, {x, 0, 1},
  PlotRange -> Automatic], {y, 0, 1}, {a, 0, 1}, {b, 0, 1},
AnimationRunning -> False]

Animate[Plot[((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2, {y, 0, 1},
  PlotRange -> Automatic], {x, 0, 1}, {a, 0, 1}, {b, 0, 1},
AnimationRunning -> False]

Animate[Plot3D[((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2, {a, 0, 1}, {b,
    0, 1}, PlotRange -> Automatic], {x, 0, 1}, {y, 0, 1},
AnimationRunning -> False]

nkunku 发表于 2011-8-11 05:02
Animate[Plot[((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2, {x, 0, 1},
  PlotRange -> Automatic], {y, 0, ...

好的，非常感谢您，如有什么需要，说一声。

nkunku 发表于 2011-8-10 20:03
附件是假设的模型ppp.mdl，其中设
a，b，x，y均为【0,1】之间的正态分布随机数。即
a=RANDOM UNIFORM(0 ...

这个结构图，也挺怪的，不能进行“reality  checks"

nkunku 发表于 2011-8-10 20:03
附件是假设的模型ppp.mdl，其中设
a，b，x，y均为【0,1】之间的正态分布随机数。即
a=RANDOM UNIFORM(0 ...

这里的12345是什么意思？
另外如何设计当z小于1时，x、y不断增加，当z等于1时，x、y停止。
不知道是否有这样的例子？
谢谢！