matlab如何给分式加限制条件,数据拟合加限定条件

本帖最后由 dabinlee 于 2018-1-8 22:41 编辑

本人用lsqcurvefit对阻抗模型进行数据拟合，拟合公式Z=a(1)+(2000.*pi.*xdata).*a(2).*1i+(((((((a(5).^(-1))+(((2000.*pi.*xdata).*a(6).*1i).^(-1))).^(-1))+(2000.*pi.*xdata).*a(4).*1i).^(-1))+(a(3).^(-1))).^(-1))，其中xdata为自变量，Z为因变量。a(1)~a(6)应均为正实数，但拟合出来的结果a(1)~a(6)都为复数。恳请问一下如何加限制条件？或者事先定义a(1)~a(6)为正实数？十分感谢！

代码如下

%输入选定频率值xdata/kHz

xdata=[0.02 1 2 3 4 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100];

wL=[0.021296228

0.825987526

1.635136116

2.428765239

3.201157196

3.947096943

7.27215855

10.04304322

12.51484828

14.85659141

17.14744073

19.41818386

21.67950221

23.93705066

26.18831591

28.43047055

30.66445705

32.89310284

35.10666899

37.31269531