最小二乘法拟合函数

kekek · 发表于 2021-2-1 16:33

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

我的计算中使用了最小二乘法拟合函数，因此自定义了一个函数，为了能把主函数中定义的变量直接拿过来在自定义函数中使用我就使用了global全局变量，第一次使用global变量，为什么在自定义函数中一使用global，就显示“此处使用的FUNCTION关键字无效”，求助
%主程序部分
clc;close all;clear all;
global wavelength_0 intensity OPD;% unit: m
wavelength_0=532*10^-9;
data=xlsread('F:\扫描结果\aa.xlsx');%读数据
intensity = data(:,1);
intensity = intensity'./max(intensity);
intensity_number = length(intensity);
OPD= linspace(0,intensity_number./30,intensity_number);
figure(1)
plot(OPD,intensity);
a0 = [0.133,0,2]; %初始估计值
a = lsqcurvefit(@myfun,a0,OPD,intensity);
F = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
coherence = (max(newtype)-min(newtype))./(max(newtype)+min(newtype));
figure(2)
plot(OPD,coherence);

%下面是我新建了一个.m文件的自定义函数
global wavelength_0 OPD;
function F = myfun(a,OPD)%就是这里出问题的
F = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
end

Uifhjvv · 发表于 2021-2-1 17:03

顺序有问题
%下面是我新建了一个.m文件的自定义函数

function F = myfun(a,OPD)%就是这里出问题的
global wavelength_0
F = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
end

zaiyiaaaa · 发表于 2021-2-2 11:20

来学习一下

小白的白 · 发表于 2021-2-2 11:23

帐号		自动登录	找回密码
密码			注册

最小二乘法拟合函数

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1