负载阻抗:Zl=a+j*b(a=50,b=100)
C1:电容C1  C1=0:0.0001:3
C2:电容C2  C2=0:0.0001:5
L:电感L   L=5
Z0:传输线特征阻抗Z0 50

w:工作频率  13560000Hz
S1:负载到匹配网络之间传输线长度  0.1
S2:检测模块到匹配网络之间传输线路长  0.5
B:相速率
j:虚数
负载点反射系数:Fl=(1/(1/(a+j*b+1/(j*w*C1)+j*w*L)+j*w*C2)-Z0)/(1/(1/(a+j*b+1/(j*w*C1)+j*w*L)+j*w*C2)+Z0)

测量点反射系数:Fsk=(1/(j*w*C1/(Z0*(1+exp(-2*j*B*S1)*Fl)/(1-exp(-2*j*B*S1)*Fl)*j*w*C1+1-w*w*L*C1)+j*w*C2)-Z0)/(1/(j*w*C1/(Z0*(1+exp(-2*j*B*S1)*Fl)/(1-exp(-2*j*B*S1)*Fl)*j*w*C1+1-w*w*L*C1)+j*w*C2)+Z0)*exp(-2*j*B*S2)
请大侠帮忙给个matlab的程序,是反射系数与两个电容之间的三维图形,反射系数为y轴的,里面的参数不重要,随便设定,关键是过程,我想知道怎么画图,2个电容都在变化,而不是只有一个变化,即两个变量~~~

小菜在此谢过~~~