谢谢celestialgod！
简单来说，我有三个公式想要记算，分别是S1, F2, A2的算式
但各自公式成立之前都有一些不一样的限制式需要被满足，
如果限制式有 > 0才计算数值，不然就显示0。
所以，
1. 在S1的前面应该要有：
if......(S1成立的条件式)
S1 = ......
else
S1 = 0
if......(F2成立的条件式)
F2 = ......
else
F2 = 0
if......(A2成立的条件式)
A2 = ......
else
A2 = 0
2. 这样一轮下来应该就会有S1, F2, A2的值 (有数字或是0)
接着去比较在b = 0, D = 0.1的状况下哪个是Max并且print name出来
然后循环后在b = 0.1, D = 0.1的状况下再比较一次
...经过很多轮
然后循环后在b = 1, D = 1的状况下再比较一次
3. 最后希望可以整理出在不同b跟D的组合下，什么时候S1, F2, A2会是max
然后想要画出一个图http://i.imgur.com/NNTHMT9.jpg
像右边图，在不同的组合下面max的值是哪一个
不好意思再麻烦各位帮我解惑了 > <
谢谢！
※ 引述《celestialgod (天)》之铭言：
有些问题看得不是很懂...
有问题推文问再回你吧~"~
q=0.9;
a=0.9;
c0=0.2;
c1=0.01;
X=2^(1/2);
[b, D] = meshgrid(0:0.1:1, 0.1:0.1:1);
S1 = (q*(a+b(:))-2*c0).^2/(16*(q-c1));
F2 = (q^2*(2*D(:) + a - b(:)) .* (a*q - 2*c0 + b(:)*q))/(16*(c1 - q).^2);
A2 = ((X*q*(q*(a + b(:))-2*c0)).^2)/((4*(q-c1)).^3);
Y = [S1, F2, A2];
Y(q - b(:) > 0 & q*(q + b(:)) - 2*c0 > 0) = 0;
if q - c1 <= 0
Y(:) = 0;
end
[Z, code] = max(Y, [], 2);
names = {'S1', 'F2', 'A2'};
maxNames = arrayfun(@(x) names{x}, code, 'uni', false);
储存的话，要同时存chr跟数字的话，就用open跟fprintf解决吧，网络上很多
存数字就用dlmwrite就好
※ 引述《kevinjean (Kevin)》之铭言：
: 各位版上前辈好，
: 我想要请教几个问题希望大家可以帮助我 > <
:

1. 计算后，我再取条件等于0，不是一样？如果你要个别条件，在个别加就好2. 我不知道你在写什么...

例如(b,D,max)=(0,0.1,S1);(0.1,0.1,S1);(0.2,0.1,S1)(0.3,0.1,S1);...;(0.5,0.5,A2);(0.6,0.5,A2);......;(1,0.9,F2);(1,1,F2) 这样的输出

所以我前面不就说可以用fprintf做？前面程式OUTPUT是错的，你也没发现...你到底有没有跑过我写的程式，再来问问题？

有我有跑过了，谢谢！目前先看到的状况就是有一连串S1, F2, A2。我还没有先对output的内容，我会再研究一下，谢谢celestialgod

所以output应该都给你了，我不知道你后面要干嘛就是

限制式只限制b , 解出来刚好是一个开区间。S1 F2 A2又刚好只是b D的二次式，应该可以手算才是