版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

[求助] matlab使用自定义函数文件拟合出现错误已有1人参与

代码如下，x,y是想要拟合的曲线对应的x，y值；L1，L2是想要拟合的参数；
下面就是一系列的矩阵计算，通过x求得y
使用ft=fittype('c(x,，L1，L2)');
出现错误：
错误使用 fittype/testCustomModelEvaluation (line 12)
Expression c(x,L1,L2) is not a valid MATLAB expression, has non-scalar coefficients, or cannot be evaluated:
Error in fittype expression ==> c(x,L1,L2)
??? 索引超出矩阵维度。
出错 fittype>iCreateFittype (line 371)
testCustomModelEvaluation( obj );
出错 fittype (line 328)
obj = iCreateFittype( obj, varargin{:} );
原因:
错误使用 fittype/evaluate (line 102)
Error in fittype expression ==> c(x,L1,L2)
??? 索引超出矩阵维度。
怎么解决这个问题？还有其他拟合办法吗？求助好心的坛友~~~~~~~~~~~~~~~~~~~

function [ y ] = c(x,L1，L2)
%UNTITLED2 此处显示有关此函数的摘要
% 此处显示详细说明
m=20; % number of section
n1 = 2.8924; %AlAs
n2 = 3.3702; %GaAS
na=1;
ns=n2;
ra=(na-n2)/(na+n2);
r=(n2-n1)/(n2+n1);
ta=sqrt(1-ra^2);
t=sqrt(1-r^2);
y=zeros(1,2500);
ab1=0;
ab2=0;
for kl=1:2500

beta1 = 2*pi*n1/x(kl)-1i*ab1/2;
beta2 = 2*pi*n2/x(kl)-1i*ab2/2;

T11_a = exp(1i*beta2*L2)/ta;
T12_a = ra*exp(-1i*beta2*L2)/ta;
T21_a = ra*exp(1i*beta2*L2)/ta;
T22_a = exp(-1i*beta2*L2)/ta;
T=[T11_a,T12_a;T21_a,T22_a];

      for k=1:m-1
         T11_n2 = exp(1i*beta1*L1)/t;
         T12_n2 = r*exp(-1i*beta1*L1)/t;
         T21_n2 = r*exp(1i*beta1*L1)/t;
         T22_n2 = exp(-1i*beta1*L1)/t;
         T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];
         T=T*T_n2;

         T11_n1 = exp(1i*beta2*L2)/t;
         T12_n1 = -r*exp(-1i*beta2*L2)/t;
         T21_n1 = -r*exp(1i*beta2*L2)/t;
         T22_n1 = exp(-1i*beta2*L2)/t;
         T_n1=[T11_n1,T12_n1;T21_n1,T22_n1];
         T=T*T_n1;

      end
      T11_n2 = exp(1i*beta1*L1)/t;
      T12_n2 = r*exp(-1i*beta1*L1)/t;
      T21_n2 = r*exp(1i*beta1*L1)/t;
      T22_n2 = exp(-1i*beta1*L1)/t;
      T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];
      T=T*T_n2;

      T_12=[1/t,-r/t;-r/t,1/t]; % from n1 to n2, just cross the interface
      T = T*T_12;

      T11_c = exp(1i*beta2*2*L2);
      T12_c = 0;
      T21_c = 0;
      T22_c = exp(-1i*beta2*2*L2);
      T_c = [T11_c,T12_c;T21_c,T22_c]; % n2 cavity
      T = T*T_c;

      for k=(1:m)
         T11_n2 = exp(1i*beta1*L1)/t;
         T12_n2 = r*exp(-1i*beta1*L1)/t;
         T21_n2 = r*exp(1i*beta1*L1)/t;
         T22_n2 = exp(-1i*beta1*L1)/t;
         T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];
         T=T*T_n2;
         T11_n1 = exp(1i*beta2*L2)/t;
         T12_n1 = -r*exp(-1i*beta2*L2)/t;
         T21_n1 = -r*exp(1i*beta2*L2)/t;
         T22_n1 = exp(-1i*beta2*L2)/t;
         T_n1=[T11_n1,T12_n1;T21_n1,T22_n1];
         T=T*T_n1;
      end


      T11_n2 = exp(1i*beta1*L1)/t;
      T12_n2 = r*exp(-1i*beta1*L1)/t;
      T21_n2 = r*exp(1i*beta1*L1)/t;
      T22_n2 = exp(-1i*beta1*L1)/t;
      T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];
      T=T*T_n2;

      T_12=[1/t,-r/t;-r/t,1/t];
      T = T*T_12;

      T11_c = exp(1i*beta2*9*L2);
      T12_c = 0;
      T21_c = 0;
      T22_c = exp(-1i*beta2*9*L2);
      T_c = [T11_c,T12_c;T21_c,T22_c]; % n2 cavity
      T = T*T_c;

      T11_n2 = exp(1i*beta1*L1)/t;
      T12_n2 = r*exp(-1i*beta1*L1)/t;
      T21_n2 = r*exp(1i*beta1*L1)/t;
      T22_n2 = exp(-1i*beta1*L1)/t;
      T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];
      T=T*T_n2;

      rs1=(n1-ns)/(ns+n1);
      ts1=sqrt(1-rs1^2);
      T11_s = 1/ts1;  % n1 to substrate
      T12_s = rs1/ts1;
      T21_s = rs1/ts1;
      T22_s = 1/ts1;
      T_s = [T11_s,T12_s;T21_s,T22_s];
      T=T*T_s;
      y(kl) = (abs(T(2,1)/T(1,1)))^2;
end

end

回复此楼

» 收录本帖的淘帖专辑推荐

学习方法

» 猜你喜欢

重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。已经有0人回复
UJN物理学专业调剂已经有12人回复
物理学I论文润色/翻译怎么收费? 已经有138人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名已经有0人回复
山西大同大学物理学硕士研究生招收调剂生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
一志愿物理所凝聚态物理352求调剂已经有0人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有1人回复

1楼 2016-08-31 17:16:30

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

512095518

专家顾问 (正式写手)

专家经验: +52
应助: 118 (高中生)
金币: 7665.4
红花: 78
帖子: 712
在线: 186.1小时
虫号: 3195576
注册: 2014-05-10
性别: GG
专业: 大气环境与全球气候变化
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

你先给出你的拟合数据，x，y，然后给出你估计的L1 L2的值
你的这个函数function [ y ] = c(x,L1，L2)，是表示用x，L1，L2就可以计算出y值对吧

赞一下

回复此楼

没有极限，只有超越！

2楼2016-08-31 19:11:07

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

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

引用回帖:

2楼: Originally posted by 512095518 at 2016-08-31 19:11:07
你先给出你的拟合数据，x，y，然后给出你估计的L1 L2的值
你的这个函数function = c(x,L1，L2)，是表示用x，L1，L2就可以计算出y值对吧

是的，我已经设置了x,y值，想通过我的函数c（x,L1,L2）来求得合适的L1与L2值

赞一下

回复此楼

3楼2016-08-31 23:37:42

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

512095518

专家顾问 (正式写手)

专家经验: +52
应助: 118 (高中生)
金币: 7665.4
红花: 78
帖子: 712
在线: 186.1小时
虫号: 3195576
注册: 2014-05-10
性别: GG
专业: 大气环境与全球气候变化
管辖: 计算模拟

【答案】应助回帖

引用回帖:

3楼: Originally posted by cs24451 at 2016-08-31 23:37:42
是的，我已经设置了x,y值，想通过我的函数c（x,L1,L2）来求得合适的L1与L2值...

我想帮你，就必须只要你的x，y数值，还有待定的L1 L2，然后我再尝试看看怎么编写代码

赞一下

回复此楼

没有极限，只有超越！

4楼2016-08-31 23:48:58

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

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

引用回帖:

4楼: Originally posted by 512095518 at 2016-08-31 23:48:58
我想帮你，就必须只要你的x，y数值，还有待定的L1 L2，然后我再尝试看看怎么编写代码...

我把想x,y L1 L2的值都写在excel里面了

多谢！！！！

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 参数.xlsx

2016-09-01 00:44:54, 69.71 K

5楼2016-09-01 00:45:30

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

512095518

专家顾问 (正式写手)

专家经验: +52
应助: 118 (高中生)
金币: 7665.4
红花: 78
帖子: 712
在线: 186.1小时
虫号: 3195576
注册: 2014-05-10
性别: GG
专业: 大气环境与全球气候变化
管辖: 计算模拟

【答案】应助回帖

引用回帖:

5楼: Originally posted by cs24451 at 2016-09-01 00:45:30
我把想x,y L1 L2的值都写在excel里面了

多谢！！！！...

L1 = -9.02438565124979e-09
L2 = -8.31832684073256e-09
你验证下对不对

赞一下

回复此楼

没有极限，只有超越！

6楼2016-09-01 13:31:34

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

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

引用回帖:

6楼: Originally posted by 512095518 at 2016-09-01 13:31:34
L1 = -9.02438565124979e-09
L2 = -8.31832684073256e-09
你验证下对不对...

额楼主怎么实现的能说下么我这个L1 L2代表长度，希望在我给出的值附近浮动，不能为负。。。

赞一下

回复此楼

7楼2016-09-01 16:13:49

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

512095518

专家顾问 (正式写手)

专家经验: +52
应助: 118 (高中生)
金币: 7665.4
红花: 78
帖子: 712
在线: 186.1小时
虫号: 3195576
注册: 2014-05-10
性别: GG
专业: 大气环境与全球气候变化
管辖: 计算模拟

【答案】应助回帖

引用回帖:

7楼: Originally posted by cs24451 at 2016-09-01 16:13:49
额楼主怎么实现的能说下么我这个L1 L2代表长度，希望在我给出的值附近浮动，不能为负。。。...

正数也一样，把负号去了就行，代码如下，我是MATLAB R2016a

CODE:

function F = myfun(b)

data = load('a.txt');

x = data(:,1); yn = data(:,2);

m=20; % number of section

n1 = 2.8924;   %AlAs

n2 = 3.3702;   %GaAS

na=1;

ns=n2;

ra=(na-n2)/(na+n2);

r=(n2-n1)/(n2+n1);

ta=sqrt(1-ra^2);

t=sqrt(1-r^2);

y=zeros(1,2500);

ab1=0;

ab2=0;

for kl=1:2500

    beta1 = 2*pi*n1/x(kl)-1i*ab1/2;

    beta2 = 2*pi*n2/x(kl)-1i*ab2/2;

    T11_a = exp(1i*beta2*b(2))/ta;

    T12_a = ra*exp(-1i*beta2*b(2))/ta;

    T21_a = ra*exp(1i*beta2*b(2))/ta;

    T22_a = exp(-1i*beta2*b(2))/ta;

    T=[T11_a,T12_a;T21_a,T22_a];

        for k=1:m-1

            T11_n2 = exp(1i*beta1*b(1))/t;

            T12_n2 = r*exp(-1i*beta1*b(1))/t;

            T21_n2 = r*exp(1i*beta1*b(1))/t;

            T22_n2 = exp(-1i*beta1*b(1))/t; 

            T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];

            T=T*T_n2;

            T11_n1 = exp(1i*beta2*b(2))/t;

            T12_n1 = -r*exp(-1i*beta2*b(2))/t;

            T21_n1 = -r*exp(1i*beta2*b(2))/t;

            T22_n1 = exp(-1i*beta2*b(2))/t;  

            T_n1=[T11_n1,T12_n1;T21_n1,T22_n1];

            T=T*T_n1;

        end

        T11_n2 = exp(1i*beta1*b(1))/t;

        T12_n2 = r*exp(-1i*beta1*b(1))/t;

        T21_n2 = r*exp(1i*beta1*b(1))/t;

        T22_n2 = exp(-1i*beta1*b(1))/t; 

        T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];

        T=T*T_n2;

        T_12=[1/t,-r/t;-r/t,1/t]; % from n1 to n2, just cross the interface

        T = T*T_12;

        T11_c = exp(1i*beta2*2*b(2));

        T12_c = 0;

        T21_c = 0;

        T22_c = exp(-1i*beta2*2*b(2));

        T_c = [T11_c,T12_c;T21_c,T22_c]; % n2 cavity

        T = T*T_c;

        for k=(1:m)

            T11_n2 = exp(1i*beta1*b(1))/t;

            T12_n2 = r*exp(-1i*beta1*b(1))/t;

            T21_n2 = r*exp(1i*beta1*b(1))/t;

            T22_n2 = exp(-1i*beta1*b(1))/t; 

            T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];

            T=T*T_n2;

            T11_n1 = exp(1i*beta2*b(2))/t;

            T12_n1 = -r*exp(-1i*beta2*b(2))/t;

            T21_n1 = -r*exp(1i*beta2*b(2))/t;

            T22_n1 = exp(-1i*beta2*b(2))/t;  

            T_n1=[T11_n1,T12_n1;T21_n1,T22_n1];

            T=T*T_n1;

        end   

        T11_n2 = exp(1i*beta1*b(1))/t;

        T12_n2 = r*exp(-1i*beta1*b(1))/t;

        T21_n2 = r*exp(1i*beta1*b(1))/t;

        T22_n2 = exp(-1i*beta1*b(1))/t; 

        T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];

        T=T*T_n2;

        T_12=[1/t,-r/t;-r/t,1/t]; 

        T = T*T_12;

        T11_c = exp(1i*beta2*9*b(2));

        T12_c = 0;

        T21_c = 0;

        T22_c = exp(-1i*beta2*9*b(2));

        T_c = [T11_c,T12_c;T21_c,T22_c]; % n2 cavity

        T = T*T_c;

        T11_n2 = exp(1i*beta1*b(1))/t;

        T12_n2 = r*exp(-1i*beta1*b(1))/t;

        T21_n2 = r*exp(1i*beta1*b(1))/t;

        T22_n2 = exp(-1i*beta1*b(1))/t; 

        T_n2=[T11_n2,T12_n2;T21_n2,T22_n2];

        T=T*T_n2;

        rs1=(n1-ns)/(ns+n1);

        ts1=sqrt(1-rs1^2);

        T11_s = 1/ts1;  % n1 to substrate

        T12_s = rs1/ts1;

        T21_s = rs1/ts1;

        T22_s = 1/ts1;

        T_s = [T11_s,T12_s;T21_s,T22_s];

        T=T*T_s;

        y(kl) = (abs(T(2,1)/T(1,1)))^2;

end

F = y'-yn;

CODE:

b0 = [-0.5e-9,-0.5e-9];

options = optimoptions(@lsqnonlin,'TolFun',1e-20,'TolX',1e-20,...

    'Algorithm','levenberg-marquardt','Plotfcns',@optimplotfval,...

    'OptimalityTolerance',1e20,'MaxIterations',1500,...

    'Display', 'iter-detailed','MaxFunctionEvaluations',2000);

K2 = lsqnonlin(@myfun,b0,[],[],options);

赞一下

回复此楼

没有极限，只有超越！

8楼2016-09-01 21:20:02

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

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

引用回帖:

8楼: Originally posted by 512095518 at 2016-09-01 21:20:02
正数也一样，把负号去了就行，代码如下，我是MATLAB R2016a
function F = myfun(b)
data = load('a.txt');
x = data(:,1); yn = data(:,2);
m=20; % number of section
n1 = 2.8924; %AlAs
n2 = 3.3702; ...

我带入了一下不太对啊这是我想要拟合的曲线

untitled.jpg

赞一下

回复此楼

9楼2016-09-02 01:04:37

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

cs24451

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 49
在线: 11.4小时
虫号: 1692275
注册: 2012-03-15
专业: 光子与光电子器件

引用回帖:

这是我原有的L1 L2带入之后的曲线

111.jpg

赞一下

回复此楼

10楼2016-09-02 01:08:01

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

相关版块跳转我要订阅楼主 cs24451 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿华南师范361分，化学求调剂 +5	Nicole88888 2026-04-01	5/250	2026-04-04 13:56 by 聪明的大松鼠
[考研] 材料383求调剂 +4	郭阳阳阳成 2026-04-04	4/200	2026-04-04 12:07 by 无际的草原
[考研] 085600材料与化工调剂 +26	kikiki7 2026-03-30	27/1350	2026-04-04 09:18 by qlm5820
[考研] 343求调剂 +8	爱羁绊 2026-03-29	8/400	2026-04-03 22:19 by 柟�?
[考研] 化工求调剂 +11	荔香芝士椰奶 2026-04-03	11/550	2026-04-03 22:06 by 啵啵啵0119
[考研] 材料专硕322分 +13	哈哈哈吼吼吼哈 2026-04-01	13/650	2026-04-03 16:08 by 哦哦123
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +5	redcom 2026-04-01	6/300	2026-04-03 15:15 by Howard28
[考研] 326分求调剂 +3	于是乎呢 2026-04-01	5/250	2026-04-03 14:23 by 于是乎呢
[考研] 一志愿生物与医药，296分，求调剂 +8	66鹿 2026-04-03	9/450	2026-04-03 14:22 by 化学化工硕士招�
[考研] 专硕085601求调剂 +7	suyifei 2026-04-03	8/400	2026-04-03 14:00 by 欣喜777
[考研] 土木水利328分求调剂 +6	疾风知劲草666 2026-04-02	6/300	2026-04-03 11:38 by znian
[考研] 能源动力调剂 +3	不破不立0 2026-04-02	3/150	2026-04-02 12:46 by ffffjjjj
[考研] 一志愿9初试366 本双非求调剂 +4	运气来得若有似� 2026-04-02	4/200	2026-04-02 09:56 by guanxin1001
[考研] 298求调剂 +4	什么是胖头鱼 2026-03-30	6/300	2026-04-01 22:06 by 客尔美德
[考研] 生物与医药考研调剂 +5	铁憨憨123425 2026-03-31	5/250	2026-04-01 18:01 by syh9288
[硕博家园] 博一被送出联培感觉不适应怎么办 +3	全村的狗 2026-03-31	3/150	2026-04-01 10:44 by 328838485
[考研] 一志愿西交大080500材料学硕349 +6	jqx1258 2026-03-31	7/350	2026-03-31 21:08 by yuq
[考研] 江苏苏北高校诚邀调剂同学 +3	zzll406 2026-03-31	3/150	2026-03-31 16:54 by 及时行乐fan
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu锂
[考研] 305求调剂 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito

24小时热门版块排行榜

[求助] matlab使用自定义函数文件拟合出现错误 已有1人参与

» 收录本帖的淘帖专辑推荐

» 猜你喜欢

【答案】应助回帖

【答案】应助回帖

» 本帖附件资源列表

【答案】应助回帖

【答案】应助回帖

[求助] matlab使用自定义函数文件拟合出现错误已有1人参与