版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

cs24451

鐵蟲 (初入文壇)

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

[求助] matlab使用自定義函數(shù)文件擬合出現(xiàn)錯誤已有1人參與

代碼如下，x,y是想要擬合的曲線對應的x，y值；L1，L2是想要擬合的參數(shù)；
下面就是一系列的矩陣計算，通過x求得y
使用ft=fittype('c(x,，L1，L2)');
出現(xiàn)錯誤：
錯誤使用 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 此處顯示有關此函數(shù)的摘要
% 此處顯示詳細說明
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

回復此樓

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

學習方法

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調劑已經(jīng)有6人回復
0702一志愿吉大B區(qū)求調劑已經(jīng)有5人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有273人回復
求調劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調劑信息已經(jīng)有31人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有65人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復
調劑已經(jīng)有0人回復

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

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

cs24451

鐵蟲 (初入文壇)

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

引用回帖:

6樓: Originally posted by 512095518 at 2016-09-01 13:31:34
L1 = -9.02438565124979e-09
L2 = -8.31832684073256e-09
你驗證下對不對...

額樓主怎么實現(xiàn)的能說下么我這個L1 L2代表長度，希望在我給出的值附近浮動，不能為負。。。

贊一下

回復此樓

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

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

查看全部 13 個回答

512095518

專家顧問 (正式寫手)

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

【答案】應助回帖

感謝參與，應助指數(shù) +1

你先給出你的擬合數(shù)據(jù)，x，y，然后給出你估計的L1 L2的值
你的這個函數(shù)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
你先給出你的擬合數(shù)據(jù)，x，y，然后給出你估計的L1 L2的值
你的這個函數(shù)function = c(x,L1，L2)，是表示用x，L1，L2就可以計算出y值對吧

是的，我已經(jīng)設置了x,y值，想通過我的函數(shù)c（x,L1,L2）來求得合適的L1與L2值

贊一下

回復此樓

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

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

512095518

專家顧問 (正式寫手)

專家經(jīng)驗: +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
是的，我已經(jīng)設置了x,y值，想通過我的函數(shù)c（x,L1,L2）來求得合適的L1與L2值...

我想幫你，就必須只要你的x，y數(shù)值，還有待定的L1 L2，然后我再嘗試看看怎么編寫代碼

贊一下

回復此樓

沒有極限，只有超越！

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

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

查看全部 13 個回答

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

[考研] 350求調劑 +4	weudhdk 2026-03-19	4/200	2026-03-20 21:06 by 無際的草原
[考研] 296求調劑 +5	www_q 2026-03-18	9/450	2026-03-20 20:52 by zhukairuo
[考研] 085700資源與環(huán)境308求調劑 +11	墨墨漠 2026-03-18	12/600	2026-03-20 19:43 by 丁丁*
[考研] 環(huán)境工程調劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風里
[考研] 317求調劑 +4	申子申申 2026-03-19	8/400	2026-03-20 11:20 by 申子申申
[考研] 求調劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學校友
[考研] 081700化工學碩調劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 0703化學調劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 材料工程專碩調劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 328求調劑，英語六級551，有科研經(jīng)歷 +3	生物工程調劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 295求調劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 311求調劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 070300化學319求調劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 312求調劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 280求調劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 考研求調劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 0856專碩279求調劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 070305求調劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

cs24451

[求助] matlab使用自定義函數(shù)文件擬合出現(xiàn)錯誤 已有1人參與

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

» 猜你喜歡

cs24451

512095518

【答案】應助回帖

cs24451

512095518

【答案】應助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] matlab使用自定義函數(shù)文件擬合出現(xiàn)錯誤已有1人參與