| 4 | 1/1 | 返回列表 |
| 查看: 598 | 回復(fù): 3 | ||
阿德金蟲 (正式寫手)
|
[求助]
擬合的matlab程序,幫忙看看有什么問(wèn)題,運(yùn)行出錯(cuò) 已有1人參與
|
|
大家?guī)臀铱纯次业倪@個(gè)matlab程序有什么問(wèn)題沒有? 我是做擬合,I1,Q1為數(shù)據(jù)點(diǎn),Ip為關(guān)于q的函數(shù) clc clear all; syms q fai beita d N det_p det_pu; Q1=[0.0253 0.026 0.0267 0.0275 0.0286 0.0298 0.0311 0.0326 0.0346 0.0364 0.0395 0.0423 0.0457 0.0472 0.0492 0.052 0.0538 0.0554 0.0571 0.0577 0.0606 0.0633 0.0645 0.0649 0.0658 0.0671 0.0685 0.07 0.0715 ]; I1=[0.0275 0.0263 0.025 0.0238 0.0226 0.0215 0.0203 0.019 0.018 0.0169 0.0159 0.0151 0.0152 0.0153 0.0163 0.017 0.018 0.0191 0.0199 0.0201 0.0197 0.019 0.0183 0.0179 0.0171 0.016 0.0148 0.0137 0.0126 ]; Ip=ones(1,200);%預(yù)先分配空間 q=(0.02:0.001:0.08); deta_c=fai*beita*d; deta_a=(1-fai)*beita*d; xc=fai*d; xa=(1-fai)*d; Fc=exp(-1i.*xc.*q).*exp(-0.5*deta_c^2.*q.^2);%對(duì)fc的傅里葉變換的結(jié)果fc=(1/(2*pi*deta_c^2))*exp(-((x-xc)^2/(2*deta_c^2))) Fa=exp(-1i.*xa.*q).*exp(-0.5*deta_a^2.*q.^2); %對(duì)fa進(jìn)行傅里葉變換的結(jié)果fa=(1/(2*pi*deta_a^2))*exp(-((x-xa)^2/(2*deta_a^2))) R1=(2*det_p^2./(q.^2)).*real(((1-Fc).*(1-Fa))./(1-Fc.*Fa)+(1/N).*Fa.*((1-Fc).^2./(1-Fa.*Fc).^2).*(1-(Fa.*Fc).^N)); R2=(2*det_p^2)./(q.^2*N).*real((det_pu^2)/(det_p^2).*(1-(Fc.*Fa).^N)-(det_pu/det_p).*((1+Fa).*((1-Fc)./(1-Fc.*Fa).*(1-(Fc.*Fa).^N)))); Ip=(R1+R2)./(4*pi*q.^2*50000);; plot(q,Ip,'-r',Q1,I1,':b') xlim([0 0.1]) ylim([0.01,0.05]) grid |
金蟲 (正式寫手)
木蟲 (正式寫手)
金蟲 (正式寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|