Question

我正试图通过此函数传递I，V，H的三个cell2mat数组，并绘制下面来自nlinfit模型的参数结果。但是，当代码运行时，它什么也不做，只存储一个值。任何帮助表示赞赏：）

function [Icp] = Fraunhofer_Function(I,V,H)        

V1 = @(b,I)(b(1).*sign(I).*real(sqrt(I.^2 - (sign(I).*(   (b(2)+b(3)/2)  )).^2)) + b(4));
Vthresx = find(V<=1e-3 & V>=0);
Ithresvec = max(I(Vthresx));

Voffsetx = find(I<=0.1e-3 & I>=-.1e-3);
Voffset = max(V(Voffsetx));

Rn =  (max(V)-min(V))/(max(I)-min(I));

beta1 = [Rn; Ithresvec; -Ithresvec; Voffset]; %Init values  b1=Rn b2 = Icp, b3 = Icm, b4 = Voffset
opts = statset('MaxIter', 500000, 'MaxFunEvals', 100000, 'RobustWgtFun', 'andrews');
B1 = nlinfit(I, V, V1, beta1, opts ); %Fit 
Icp = V1(B1,V);


end



  files = dir('*.xlsx*');
  for k =1:length(files)
  filenames = files(k).name;

  txt = 'I,V,H';
  [num,txt,raw] = xlsread(filenames);

  %Put data into numerical columns
  Idata = num(:,1)'; Vdata = num(:,2)'; Hdata = num(:,3)';        
  [Hu,~,idx] = unique(Hdata);
  Isplit = splitapply(@(x) {x}, [Idata(:)],idx);
  Vsplit = splitapply(@(x) {x}, [Vdata(:)],idx);
  Hsplit = splitapply(@(x) {x}, [Hdata(:)],idx);    

  for l = 1:length(Isplit)

      I = (Isplit{l,1});
      V = (Vsplit{l,1});
      H = (Hsplit{l,1});

     %fit the data to the functional form
      Icp = Fraunhofer_Function(I,V,H);              

  end

结束

I，V，H数据示例如下：enter image description here

Answer 1

现在，您只需将Icp设置为第二个估计系数。在Matlab文档nlinfit中，输出beta：

beta —估计回归系数   向量

估计的回归系数，以向量形式返回。 beta中的元素数量等于beta0中的元素数量。

因此，要使用估计的参数，应调用modelfun，并将参数存储在B1中：

Icp = V1(B1,V);
plot(H,Icp);

Matlab：从nlinfit返回多个参数值（现在我只遍历循环的一个索引）

1 个答案: