cs24451

鐵蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時
蟲號: 1692275
注冊: 2012-03-15
專業(yè): 光子與光電子器件

[求助] 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

回復此樓

» 收錄本帖的淘帖專輯推薦

學習方法

» 猜你喜歡

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

512095518

專家顧問 (正式寫手)

專家經驗: +52
應助: 118 (高中生)
金幣: 7665.4
紅花: 78
帖子: 712
在線: 186.1小時
蟲號: 3195576
注冊: 2014-05-10
性別: GG
專業(yè): 大氣環(huán)境與全球氣候變化
管轄: 計算模擬

【答案】應助回帖

引用回帖:

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發(fā)消息送TA紅花 TA的回帖

查看全部 13 個回答

512095518

專家顧問 (正式寫手)

專家經驗: +52
應助: 118 (高中生)
金幣: 7665.4
紅花: 78
帖子: 712
在線: 186.1小時
蟲號: 3195576
注冊: 2014-05-10
性別: GG
專業(yè): 大氣環(huán)境與全球氣候變化
管轄: 計算模擬

【答案】應助回帖

感謝參與，應助指數 +1

你先給出你的擬合數據，x，y，然后給出你估計的L1 L2的值
你的這個函數function [ y ] = c(x,L1，L2)，是表示用x，L1，L2就可以計算出y值對吧

贊一下

回復此樓

沒有極限，只有超越！

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

cs24451

鐵蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時
蟲號: 1692275
注冊: 2012-03-15
專業(yè): 光子與光電子器件

引用回帖:

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發(fā)消息送TA紅花 TA的回帖

cs24451

鐵蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時
蟲號: 1692275
注冊: 2012-03-15
專業(yè): 光子與光電子器件

引用回帖:

4樓: Originally posted by 512095518 at 2016-08-31 23:48:58
我想幫你，就必須只要你的x，y數值，還有待定的L1 L2，然后我再嘗試看看怎么編寫代碼...

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

多謝�。。�！

贊一下

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯系郵箱：xiaomuchong@tal.com
附件 1 : 參數.xlsx

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

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿西南交通專碩材料355 本科雙非求調劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 求調劑 +5	Mqqqqqq 2026-03-19	5/250	2026-03-20 20:46 by zhukairuo
[考研] 一志愿西安交通大學學碩 354求調劑 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 材料080500調劑求收留 +8	一顆meteor 2026-03-13	8/400	2026-03-20 20:00 by zhukairuo
[考研] 一志愿吉林大學材料學碩321求調劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 261求B區(qū)調劑，科研經歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 招收調劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	14/700	2026-03-20 09:30 by kkcoco25
[考研] 085410人工智能專碩317求調劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 081700化工學碩調劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 0703化學調劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 324分 085600材料化工求調劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 346求調劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 0703化學調劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 301求調劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 290求調劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 中科院材料273求調劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 復試調劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen