| 4 | 1/1 | 返回列表 |
| 查看: 599 | 回復: 3 | ||
阿德金蟲 (正式寫手)
|
[求助]
擬合的matlab程序,幫忙看看有什么問題,運行出錯 已有1人參與
|
|
大家?guī)臀铱纯次业倪@個matlab程序有什么問題沒有? 我是做擬合,I1,Q1為數(shù)據(jù)點,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);%對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); %對fa進行傅里葉變換的結(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 5/250 |
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 收復試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 0817化學工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|