|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
我的计算中使用了最小二乘法拟合函数,因此自定义了一个函数,为了能把主函数中定义的变量直接拿过来在自定义函数中使用我就使用了global全局变量,第一次使用global变量,为什么在自定义函数中一使用global,就显示“此处使用的FUNCTION关键字无效”,求助2 A& {$ f3 L5 t" G) D
%主程序部分
4 |- K2 F6 x' h' ^0 Qclc;close all;clear all;
( ~$ z) ]0 N6 L/ T. d& \8 k0 aglobal wavelength_0 intensity OPD;% unit: m
; v7 \2 G% `" s/ \% Wwavelength_0=532*10^-9;$ a/ t! Q! v: `$ k
data=xlsread('F:\扫描结果\aa.xlsx');%读数据5 w0 {' N& ^' q) I; y
intensity = data(:,1);
% w2 U7 n7 b2 y" \8 T2 ^intensity = intensity'./max(intensity);
* M2 ^- ~ D6 B. Vintensity_number = length(intensity);4 |- h X9 A2 T. v
OPD= linspace(0,intensity_number./30,intensity_number);
4 B% m" c% K% q% Hfigure(1)2 h: \4 A( l2 R! c/ i
plot(OPD,intensity);
& R1 A$ u) G# [0 ca0 = [0.133,0,2]; %初始估计值6 {4 x6 ^: P6 |8 b P7 i
a = lsqcurvefit(@myfun,a0,OPD,intensity);
& j/ b9 {0 A2 uF = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));! l* g3 h& K4 Z* r' r
coherence = (max(newtype)-min(newtype))./(max(newtype)+min(newtype));2 C7 u" y6 B" O! A9 Q
figure(2)
7 x8 a: X% c+ k* v4 cplot(OPD,coherence);
2 b0 A" r w( K$ p5 v& d4 a! |
' N% z* P1 i# `9 l9 M1 B# J8 _) R4 a5 G* ]# K
%下面是我新建了一个.m文件的自定义函数
I5 Y3 R6 L+ f# I e- R. }global wavelength_0 OPD;
+ a1 I; B- I' ~9 k) Nfunction F = myfun(a,OPD)%就是这里出问题的# d; c# ~8 _; \* j/ A# J3 w
F = a(1)+a(2)*cos(2*pi.*OPD./wavelength_0+a(3));
$ |; g+ Z0 N$ P- ] a2 j9 Fend
8 j' ^& Y! C; _/ P H9 |6 Y9 i$ {$ P d1 F1 w
' A( t: B1 r$ `; H1 l0 a& B |
|