24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2011-11-15 09:35:30

已阅关注TA 给TA发消息送TA红花 TA的回帖

yqx1985

木虫 (著名写手)

云中仙

应助: 34 (小学生)
贵宾: 0.032
金币: 1608.8
散金: 1180
红花: 33
帖子: 1544
在线: 620.3小时
虫号: 550008
注册: 2008-04-22
性别: GG
专业: 药物设计与药物信息

【答案】应助回帖

★ ★ ★
余泽成(金币+3): 2011-11-19 12:19:58
余泽成: 谢谢参与应助！ 2011-11-19 12:20:07

简单的来讲，这种问题可以采用如下方式估计，不过具体的函数参数设置还要你自己摸索。
（1）建立一个函数，名为fun
R=fun(p,data) %
x=data(:,1);y1=data(:,2);y2=data(:,3);
n0 = 1; % air
n1 = p(1)+p(2)./x.^2+p(3)./x.^2; % oxide
n2 = 3.4693+130149./x.^2+2.851e10./x.^4 + 1i*(0.52686-320950./x.^2+6.2225e10./x.^4); % silicon
d1 = p(4);
%Refracted angles
tet0 = pi*60/180; % incident beam angle
tet1 = asin((n0./n1)*sin(tet0));
tet2 = asin((n1./n2).*sin(tet1));

%Boundary reflection coefficients
r01s = (n0.*cos(tet0)-n1.*cos(tet1))./(n0.*cos(tet0)+n1.*cos(tet1));
r12s = (n1.*cos(tet1)-n2.*cos(tet2))./(n1.*cos(tet1)+n2.*cos(tet2));
r01p = (n1.*cos(tet0)-n0.*cos(tet1))./(n1.*cos(tet0)+n0.*cos(tet1));
r12p = (n2.*cos(tet1)-n1.*cos(tet2))./(n2.*cos(tet1)+n1.*cos(tet2));
dt1 = (2*pi./x).*n1.*d1.*cos(tet1);
%Reflection coefficients total
rs = (r01s+r12s.*exp(2*1i*dt1))./(1+r01s.*r12s.*exp(2*1i*dt1));
rp = -(r01p+r12p.*exp(2*1i*dt1))./(1+r01p.*r12p.*exp(2*1i*dt1));
F = rp./rs;
psi=atan(abs(F))*180/pi;
delta = angle(F)*180/pi;
m=find(delta<0);
delta(m(1:end))=delta(m(1:end))+360;

% y1=psi;
% y2=delta;

R1=(y1-psi).^2;
R2=(y2-delta).^2;
R=R1+R2;
% R1，R2和R的计算很重要，会影响到最终的参数估计，你可以%更具你的需要选一个最稳健的表达方法；这样，如果能选取到合适的参数的话，得到的R应该是一个尽可能接近于0的向量

options=statset; %根据需要设定参数，很重要
R=zeros(legth(data),1);
[beta,r,J,Sigma,mse] = nlinfit(data,R, @fun,p,options); % p是初始估计值

整体思路就是这个样子的，利用非线性回归的方法解决。中间的很多参数你得好好调整下。