24小时热门版块排行榜    

查看: 1036  |  回复: 3

geegee

铜虫 (小有名气)

[求助] matlab 解方程组并作图 -----望高人指导。。。

%-------------------赋值---------------------
aa2=1.0*10^(-8);
aa3=3.0*10^(-6);
aa7=5.0*10^(3);
aa10=1.0;
aa11=1.0*10^(3);
A=5.00;
B=0.41;

%------求解非线性方程组-------------
0.5*aa2*A+0.5*aa3*(B-C)*A-0.5*aa11*x0*y0+aa7*A^2*x0^2=0
-aa10*(B-C)*y0+aa8*A^2*x0-aa11*x0*y0=0
%结果要满足 x0〉0且,y0〉0

%------将方程组的x0,y0带入a11,a12,a21,a22中-------------
a11=-0.5*aa11*y0+2*aa7*A^2*x0;
a12=-0.5*aa11*x0;
a21=aa8*A^2-aa11*y0;
a22=-aa10*(B-C)-aa11*x0;

T0=a11+a22;
M=a11*a22-a12*a21;
D0=T0*T0-4*M ;

%-------------------参数假定----------------------
%C的范围为0-0.41;其中每0.01取一个值
[0:0.01:0.41];

%-----------------作图----------------------
%求解T0=0时的C与aa8的关系,并将其作图;此处 作图可能不止一条曲线,有几条画几条 。
%求解D0=0时的C与aa8的关系,并将其作图;此处 作图可能不止一条曲线,有几条画几条 。
%将上述画在同一图中 。
        

万分感谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

拥有梦想只是一种智力,实现梦想才是一种能力。
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
dbb627: 金币+1, 谢谢指导 2013-09-24 09:41:42
geegee: 金币+10, ★★★很有帮助 2013-09-24 10:31:37
1:
T0=0:

循环常数 c        x0        y0        aa8
0        6.91945363771624E-6        0.00262588791158269        6.32872946755882
0.01        6.80873021834862E-6        0.00259074764873762        6.19169053546243
0.02        6.69735697882898E-6        0.00255528377545683        6.05417524112168
0.03        6.58530651556797E-6        0.00251948264475285        5.91617090672901
0.04        6.47254956743838E-6        0.00248332968458431        5.77766426061978
0.05        6.3590548511697E-6        0.00244680931588251        5.63864139780377
0.06        6.24478886324493E-6        0.00240990485389599        5.49908773319192
0.07        6.12971564254713E-6        0.00237259838998847        5.35898794626507
0.08        6.01379654337403E-6        0.00233487067860027        5.21832593102055
0.09        5.89698991999917E-6        0.00229670098015959        5.07708472515348
0.1        5.77925081038936E-6        0.0022580669035739        4.93524644272839
0.11        5.66053055668418E-6        0.00221894421722872        4.79279219445677
0.12        5.540776372095E-6        0.00217930663330331        4.64970199918568
0.13        5.41993083691825E-6        0.00213912555678529        4.50595468348301
0.14        5.29793131851213E-6        0.00209836979661904        4.36152776939003
0.15        5.17470929202234E-6        0.00205700522742713        4.21639734650247
0.16        5.05018954062078E-6        0.00201499439122915        4.070537925538
0.17        4.92428921662681E-6        0.00197229602988015        3.92392227176284
0.18        4.79691672338061E-6        0.00192886452824355        3.77652121247993
0.19        4.66797037730129E-6        0.00188464924789604        3.62830341389354
0.2        4.53733679851163E-6        0.00183959372565879        3.47923512156776
0.21        4.40488895438558E-6        0.00179363469928402        3.32927985543737
0.22        4.27048376333247E-6        0.00174670091413957        3.17839804935821
0.23        4.13395913049182E-6        0.00169871164698493        3.02654662135717
0.24        3.99513024304008E-6        0.00164957486103396        2.87367845682913
0.25        3.85378488130079E-6        0.00159918487088779        2.71974177999722
0.26        3.70967741939141E-6        0.00154741935485692        2.5646793828947
0.27        3.56252101537885E-6        0.00149413546565867        2.4084276641454
0.28        3.41197730159543E-6        0.00143916469619452        2.2509154171945
0.29        3.25764250356821E-6        0.00138230596677697        2.0920622750583
0.3        3.0990283590322E-6        0.00132331612279804        1.93177667786227
0.31        2.93553523680556E-6        0.00126189654792917        1.76995315806057
0.32        2.76641316782822E-6        0.00119767375757846        1.60646862049202
0.33        2.59070340491196E-6        0.00113017029564616        1.4411770873932
0.34        2.4071471171134E-6        0.00105875926432247        1.2739019940793
0.35        2.21403536468675E-6        0.000982589611614001        1.10442436143361
0.36        2.00894635744411E-6        0.000900455286007166        0.93246353729042
0.37        1.78824524208166E-6        0.000810546130556664        0.757643296046054
0.38        1.54601589310102E-6        0.000709915914764306        0.579425317135365
0.39        1.27135669842591E-6        0.000593135635816103        0.39695546326565
0.4        9.39253903613955E-7        0.00044774844399975        0.208592539036411
0.41        4.49013255114144E-7        0.000223608601046844        0.00894434403876471

2: D0=0时,多解,不稳定。
2楼2013-09-23 22:15:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

geegee

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-09-23 22:15:22
1:
T0=0:

循环常数 c        x0        y0        aa8
0        6.91945363771624E-6        0.00262588791158269        6.32872946755882
0.01        6.80873021834862E-6        0.00259074764873762        6.19169053546243
0.02        6.69735697882898E-6        0.00255528377545 ...

把运行的源程序给我吧。
1:T0=0:求出C,aa8。
2:D0=0时,多解,不稳定。此处求出同一个C时,不同的aa8。应该就是至少两条曲线。


非常感谢。
拥有梦想只是一种智力,实现梦想才是一种能力。
3楼2013-09-24 08:34:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

geegee

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-09-23 22:15:22
1:
T0=0:

循环常数 c        x0        y0        aa8
0        6.91945363771624E-6        0.00262588791158269        6.32872946755882
0.01        6.80873021834862E-6        0.00259074764873762        6.19169053546243
0.02        6.69735697882898E-6        0.00255528377545 ...

请把源程序给我吧。。。。。谢谢了
拥有梦想只是一种智力,实现梦想才是一种能力。
4楼2013-09-24 10:32:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 geegee 的主题更新
信息提示
请填处理意见