最小二乘法拟合函数

kekek

我的计算中使用了最小二乘法拟合函数，因此自定义了一个函数，为了能把主函数中定义的变量直接拿过来在自定义函数中使用我就使用了global全局变量，第一次使用global变量，为什么在自定义函数中一使用global，就显示“此处使用的FUNCTION关键字无效”，求助
%主程序部分
clc;close all;clear all;
6 k1 i# r* {2 H1 ^- C* z+ Zglobal wavelength_0 intensity OPD;% unit: m
wavelength_0=532*10^-9;
" G0 J( {, }. A% l4 vdata=xlsread('F:\扫描结果\aa.xlsx');%读数据
intensity = data(:,1);
9 ~: r- V/ O1 z1 m3 lintensity = intensity'./max(intensity);
intensity_number = length(intensity);
( F8 ?2 C* d  m  n* p! q3 S) JOPD= linspace(0,intensity_number./30,intensity_number);
figure(1)
5 M# I  J8 V3 }3 B# E: Oplot(OPD,intensity);
  S% a2 P5 p5 O0 ]a0 = [0.133,0,2];    %初始估计值
a = lsqcurvefit(@myfun,a0,OPD,intensity);
9 `/ p! W* i# F7 y5 @2 JF = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
. G" C1 P0 V8 U/ A  ^* R  qcoherence = (max(newtype)-min(newtype))./(max(newtype)+min(newtype));
figure(2)
3 Y2 m# B2 i( @2 P. ~3 _' V5 vplot(OPD,coherence);
9 k% D% x! `! ^% N3 d

%下面是我新建了一个.m文件的自定义函数
global wavelength_0 OPD;
function F = myfun(a,OPD)%就是这里出问题的
! D3 J" ]$ P% a% R2 _1 s# xF = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
9 m8 q6 `3 ^% w! Lend
, x; O" K% P- p0 e/ H! j1 ^
5 ?. a  q$ W' a3 g9 ^# p

Uifhjvv

顺序有问题
1 p6 D, o- y% D5 v7 y( |%下面是我新建了一个.m文件的自定义函数
) F5 Y- h! p' c/ E- G3 @  W
1 N% O  T+ Q& N& }( G9 P9 M$ Ifunction F = myfun(a,OPD)%就是这里出问题的
global wavelength_0  
F = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
4 Y& s# g/ e: ]/ [4 mend

zaiyiaaaa

来学习一下

小白的白