| 5 | 1/1 | 返回列表 |
| 查看: 604 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
阿德金蟲(chóng) (正式寫(xiě)手)
|
[求助]
擬合的matlab程序,幫忙看看有什么問(wèn)題,運(yùn)行出錯(cuò) 已有1人參與
|
|
|
大家?guī)臀铱纯次业倪@個(gè)matlab程序有什么問(wèn)題沒(méi)有? 我是做擬合,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 |
金蟲(chóng) (正式寫(xiě)手)
金蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6 | 墨墨漠 2026-03-20 | 6/300 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專(zhuān)碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
[考研] 機(jī)械專(zhuān)碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|