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

nlm0402 发表于 2011-8-11 09:53
这里的12345是什么意思？
另外如何设计当z小于1时，x、y不断增加，当z等于1时，x、y停止。
不知道是否有 ...

您真是太客气了，真不敢当啊。不要再给论坛币了。

1、函数中的12345用以确定均匀分布随机数的序列。不同的值就会产生不同的随机数序列。

2、关于“如何设计当z小于1时，x、y不断增加，当z等于1时，x、y停止”，见附件模型new.rar。
将z的表达式改为：
如果((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2小于1，那么z=((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2；
若否，则z=0。

即在模型中z被改变为：
IF THEN ELSE(((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2<1, ((a*x + b*y)*(x*y)/(0.5*x + 0.5*y)^2)^2 , 0 )

根据楼主对模型的调试过程，现将有关情况说明如下：

约束条件：
（1）0<a<1;
（2）0<b<1;
（3）a+b=1;
（4）0.1<x<1（单调递增）;
（5）0.1<y<1（单调递增）;
（6）当z<=1时，z=((a*x + b*y)^(1/2)*(x*y)/(0.5*x + 0.5*y)^2);否则z=0

那么
a=RANDOM UNIFORM(0,1,12345)
b=1-a
x=y=RAMP((1-0.1)/(FINAL TIME-INITIAL TIME),INITIAL TIME,FINAL TIME)+0.1
z=IF THEN ELSE(((a*x + b*y)^(1/2)*(x*y)/(0.5*x + 0.5*y)^2)<1, ((a*x + b*y)^(1/2)*(x*y)/(0.5*x + 0.5*y)^2) , 0 )
FINAL TIME=100
INITIAL TIME=0

于是，模型结构如图形vvv所示。

z的模拟结果如图形zzz所示。

nkunku 发表于 2011-8-10 22:59
clear
a=0.24;
b=0.54;

Warning: Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 2.672028e-020.
出现了这样的警告，不知该如何修改，谢谢！