版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

cs24451

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時(shí)
蟲(chóng)號(hào): 1692275
注冊(cè): 2012-03-15
專業(yè): 光子與光電子器件

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

代碼如下，x,y是想要擬合的曲線對(duì)應(yīng)的x，y值；L1，L2是想要擬合的參數(shù)；
下面就是一系列的矩陣計(jì)算，通過(guò)x求得y
使用ft=fittype('c(x,，L1，L2)');
出現(xiàn)錯(cuò)誤：
錯(cuò)誤使用 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)
??? 索引超出矩陣維度。
出錯(cuò) fittype>iCreateFittype (line 371)
testCustomModelEvaluation( obj );
出錯(cuò) fittype (line 328)
obj = iCreateFittype( obj, varargin{:} );
原因:
錯(cuò)誤使用 fittype/evaluate (line 102)
Error in fittype expression ==> c(x,L1,L2)
??? 索引超出矩陣維度。
怎么解決這個(gè)問(wèn)題？還有其他擬合辦法嗎？求助好心的壇友~~~~~~~~~~~~~~~~~~~

function [ y ] = c(x,L1，L2)
%UNTITLED2 此處顯示有關(guān)此函數(shù)的摘要
% 此處顯示詳細(xì)說(shuō)明
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

回復(fù)此樓

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

學(xué)習(xí)方法

» 猜你喜歡

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

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

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

cs24451

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時(shí)
蟲(chóng)號(hào): 1692275
注冊(cè): 2012-03-15
專業(yè): 光子與光電子器件

引用回帖:

8樓: Originally posted by 512095518 at 2016-09-01 21:20:02
正數(shù)也一樣，把負(fù)號(hào)去了就行，代碼如下，我是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; ...

我?guī)肓艘幌?nbsp; 不太對(duì)啊這是我想要擬合的曲線

untitled.jpg

贊一下

回復(fù)此樓

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

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

查看全部 13 個(gè)回答

512095518

專家顧問(wèn) (正式寫(xiě)手)

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

【答案】應(yīng)助回帖

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

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

贊一下

回復(fù)此樓

沒(méi)有極限，只有超越！

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

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

cs24451

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 49
在線: 11.4小時(shí)
蟲(chóng)號(hào): 1692275
注冊(cè): 2012-03-15
專業(yè): 光子與光電子器件

引用回帖:

2樓: Originally posted by 512095518 at 2016-08-31 19:11:07
你先給出你的擬合數(shù)據(jù)，x，y，然后給出你估計(jì)的L1 L2的值
你的這個(gè)函數(shù)function = c(x,L1，L2)，是表示用x，L1，L2就可以計(jì)算出y值對(duì)吧

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

贊一下

回復(fù)此樓

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

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

512095518

專家顧問(wèn) (正式寫(xiě)手)

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

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by cs24451 at 2016-08-31 23:37:42
是的，我已經(jīng)設(shè)置了x,y值，想通過(guò)我的函數(shù)c（x,L1,L2）來(lái)求得合適的L1與L2值...

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

贊一下

回復(fù)此樓

沒(méi)有極限，只有超越！

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

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

查看全部 13 個(gè)回答

普通表情龍兔虎貓

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

[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 299求調(diào)劑 +4	shxchem 2026-03-20	5/250	2026-03-21 00:24 by 星空星月
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	13/650	2026-03-20 21:01 by 無(wú)際的草原
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明