| 2 | 1/1 | 返回列表 |
| 查看: 317 | 回復(fù): 1 | ||
wzycrfpxz新蟲 (初入文壇)
|
[求助]
求糾錯。。。
|
![]() function guiji clc;clear all; format long; global CL CD s m hs lo0 r0 mu E0 Ef rad CL=0.3515;CD=1.4595;s=15.9;m=2804;hs=9354;lo0=0.0145;r0=3395.76*10^3;mu=4.284e+013;E0=1.5600418*10^7;... Ef=1.0766379*10^5;rad=57.325; solinit=bvpinit(linspace(0,1,100),[1 1 1 1 1]); sol=bvp4c(@f,@bcfun,solinit); plot(sol.x(1, ,sol.x(2, );grid; end function dx=f(E1,x) global CL CD s m hs lo0 r0 mu E0 Ef rad L=lo0*(2.7183^(-(x(3)-r0)/hs))*s*CL*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0))/m; D=lo0*(2.7183^(-(x(3)-r0)/hs))*s*CD*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0))/m; g=mu/x(3)^2; n=length(E1); if E1>=0&&E1<=0.24 sigma=-10/rad; elseif 0.24 elseif E1>0.4&&E1<=0.61 sigma=-80/rad; elseif E1>0.61&&E1<=0.8 sigma=(E1^2+840.69*E1-593.193)/rad; elseif E1>0.8&&E1<=0.96 sigma=80/rad; elseif E1>0.96&&E1<0.97 sigma=23.9733/rad; elseif E1>0.97&&E1<=1 sigma=-23.9733/rad; end %dx=zeros(5,10); dx=[-(Ef-E0)/x(3)*cos(x(5)/rad)*cos(x(4)/rad)/cos(x(2)/rad)/D; -(Ef-E0)/x(3)*cos(x(5)/rad)*sin(x(4)/rad)/D; -(Ef-E0)*sin(x(5)/rad)/D; (Ef-E0)*[L*sin(sigma)/D/cos(x(5)/rad)/(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))... +cos(x(5)/rad)*cos(x(4)/rad)*tan(x(2)/rad)/x(3)/D]; -(Ef-E0)/(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))*[L*cos(sigma)/D-(g... -(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))/x(3))*cos(x(5)/rad)/D]]; end function res=bcfun(xa,xb) global CL CD s m hs lo0 r0 mu E0 Ef rad res=[xa(1)+90.072;xa(2)+43.898; xa(3)-3520.76*10^3;xa(4)-7.99; xa(5)+14.15;xb(1)+73.26; xb(2)+41.45;xb(3)-3403.76*10^3]; endSample Text |
銀蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +3 | vv迷 2026-03-21 | 4/200 |
|
|
[考研] 269專碩求調(diào)劑 +5 | 金恩貝 2026-03-21 | 5/250 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 一志愿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 |
|