应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab模拟求助
52/1返回列表
查看: 1498  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

geegee

铜虫 (小有名气)

[求助] matlab模拟求助

各常量的取值如下：
a2=1.0*10^(-8);
a5=4.0*10^(4);
a7=1.0*10^(-7);
a8=1.0*10^(-2);
a10=1.0*10^(0);
a11=2.0*10^(-1);
A=5.00;
B=0.36;
三个方程组如下：
0=0.5*a2*A+0.5*a3*(B-C)*A-0.5*a11*x0*y0+a7*a5^2*A^2*x0^2;
0=-a10*(B-C)*y0+a8*a5*A^2*x0-a11*x0*y0;
0=(-1/2)*a11*y0+2*a7*a5^2*A^2*x0-a10*(B-C)-a11*x0;
变量为C,a3.
要求当C取值为[0,0.36）时，做出a3与C的关系图。
回复此楼

» 猜你喜欢
拥有梦想只是一种智力，实现梦想才是一种能力。
1楼 2013-07-05 11:33:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
17楼: Originally posted by geegee at 2013-07-05 20:43:40
抱歉，我才开始学习，有些不懂。请赐教。。。。
错误提醒如下：


??? Error using ==> optimset>checkfield at 318
Invalid value for OPTIONS parameter Algorithm:
must be 'active-set', 'trust- ...

也许版本问题吧，我是2010a，可能没你的新。
你把options=optimset('Algorithm','Levenberg-Marquardt','Display','on')中的马夸特算法即'Levenberg-Marquardt’，换成'active-set', 'trust-region-reflective', or 'interior-point‘中的任意一个好了。本来默认用trust-region-reflective算法的，后来程序警告说这个算法对方程不太适合，所以我换成了马夸特。
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
18楼2013-07-05 21:13:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

geegee

铜虫 (小有名气)
补充，x0，y0都是大于零的数值。
赞一下(1人) 回复此楼
拥有梦想只是一种智力，实现梦想才是一种能力。
2楼2013-07-05 11:35:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-05 17:25:39
geegee: 金币+20, ★★★★★最佳答案 2013-07-05 17:37:12
以下程序复制进一个m文件中，运行即可。我是按照C从0到0.36取20个点做的。

function solution_nonl_C
clear all;clc
format long

global C

%-------------根据计算结果，修改迭代初值------------------------------------
x0=[1 1 1];
C=0;
n=21;  %C的取值个数
for i=1:n
    xx=fsolve(@nonlin,x0);
   
    a3(i)=xx(1);
    x0(i)=xx(2);
    y0(i)=xx(3);
   
    C=C+0.36/(n-1);
end
[a3' x0' y0']
aa3=a3';
aaa3=[aa3(1:14);aa3(end)];
cspan=0:0.36/(n-1):0.36;
cspan1=[cspan(1:14) cspan(end)];
figure,plot(cspan,aa3,'ko',cspan1,aaa3,'b'),legend('全部点','满足x0和y0大于0的点','Location','best'),...
     xlabel('C'),ylabel('a3');

function f= nonlin(x)
global C
a2=1.0*10^(-8);
a5=4.0*10^(4);
a7=1.0*10^(-7);
a8=1.0*10^(-2);
a10=1.0*10^(0);
a11=2.0*10^(-1);
A=5.00;
B=0.36;

a3=x(1);
x0=x(2);
y0=x(3);

f(1)=0.5*a2*A+0.5*a3*(B-C)*A-0.5*a11*x0*y0+a7*a5^2*A^2*x0^2;
f(2)=-a10*(B-C)*y0+a8*a5*A^2*x0-a11*x0*y0;
f(3)=(-1/2)*a11*y0+2*a7*a5^2*A^2*x0-a10*(B-C)-a11*x0;
matlab模拟求助
附图23.jpg
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2013-07-05 14:30:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

geegee

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by 月只蓝 at 2013-07-05 14:30:54
以下程序复制进一个m文件中，运行即可。我是按照C从0到0.36取20个点做的。

function solution_nonl_C
clear all;clc
format long

global C

%-------------根据计算结果，修改迭代初值----------------- ...

能否将里面的C, '全部点', '满足x0和y0大于0的点' 输出成TXT格式。谢谢了哦。
赞一下 回复此楼
拥有梦想只是一种智力，实现梦想才是一种能力。
4楼2013-07-05 17:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭