版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

汕头大学海洋科学接受调剂

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhaoxyah

铁杆木虫 (著名写手)

应助: 3 (幼儿园)
金币: 5329.7
散金: 2779
红花: 5
帖子: 1089
在线: 974.6小时
虫号: 533803
注册: 2008-03-27
专业: 等离子体物理

[求助] 请问matlab解方程组结果如何保存已有1人参与

matlab程序如下，想把运行结果保存至一.txt文件，便于用ORIGIN画图分析，请问该如何设置如下运算结果保存，期待大侠指导
function solve_noneqs2
clear all;clc
format long
global thita
x0=[0.8 -5 0.1];
for i=1:1:88
thita=i*3.14/180.0;
[x]=fsolve(@objfun,x0)
end
function f=objfun(x)
global thita
gam=0.2;
gama=gam/(cos(thita)+gam*(1-cos(thita)));
nn=0;

mie=1836.0*4.0;
miu=1.0/mie;
tie=0.01;
rou=0.2;
uiy0=-0.01;
e0=-0.01;
f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));
f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));
f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...
+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);

回复此楼

» 猜你喜欢

恳请有学校收留已经有3人回复
335求调剂已经有21人回复
材料专业344求调剂已经有22人回复
求调剂学校已经有9人回复
复试调剂已经有20人回复
085404 22408 309分求调剂已经有6人回复
22408 312求调剂已经有13人回复
0854调剂已经有14人回复
一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药已经有8人回复
材料299专硕求调剂已经有15人回复

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

Matlab如何解同余方程组已经有3人回复
MATLAB解非线性方程组已经有9人回复
为什么matlab算出的数不正确，在线等已经有5人回复
matlab可以求解差分方程系统的解吗？已经有6人回复
matlab解方程组求助攻已经有4人回复
matlab 求解线性方程组Ax=b 已经有10人回复
Matlab求解方程‘pi*x=1’的结果为什么是1/pi? 已经有9人回复
用matlab制作的M文件保存在哪里才能直接调用？已经有9人回复
Matlab求解多元多次方程组问题已经有7人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
求助，用matlab解微分方程组，希望高手能给予指点已经有17人回复
如何用Matlab求时变微分方程的解？已经有5人回复
matlab求解高阶方程组已经有3人回复
求Matlab解非线性方程代码已经有12人回复
用matlab如何解方程？已经有6人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
求助 MATLAB解方程组-fslove 已经有7人回复
搞地学的朋友，请问怎么用matlab做Mann-Kendall趋势检验啊？已经有10人回复
求助matlab---fsolve解非线性方程组已经有6人回复
matlab计算面积求助！！难已经有8人回复
matlab中怎么取方程组的一组解并画图已经有5人回复
Matlab求解偏微分方程组已经有22人回复
matlab 计算过程中工作空间变量保存及读取求解决内存不足的方法已经有5人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2015-01-28 13:50:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhaoxyah: 金币+30, ★★★★★最佳答案 2015-01-28 14:34:46

或者更简单点，把x结果显示在主程序窗口中，复制粘贴到Origin即可：

CODE:

function solve_noneqs2

clear all;clc

format long

global thita

x0=[0.8 -5 0.1];

for i=1:1:88

    thita=i*3.14/180.0;

x(i,1:3)=fsolve(@objfun,x0);

end

x

function f=objfun(x)

global thita

gam=0.2;

gama=gam/(cos(thita)+gam*(1-cos(thita)));  

nn=0;

mie=1836.0*4.0;

miu=1.0/mie;

tie=0.01;

rou=0.2;

uiy0=-0.01;

e0=-0.01;

f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));

f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));

f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...

+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);

主程序窗口显示数据：
x =

0.000975910082976  -4.456299390614189 0.997581989816602
0.000974491887485  -4.457933083051565 0.995861494999457
0.000973373594347  -4.459510833495430 0.994138700333996
0.000972555181984  -4.461031939273196 0.992414123062480
0.000972037100394  -4.462495608478736 0.990688283093975
0.000971820034219  -4.463901127867802 0.988961492335664
0.000971905879388  -4.465247179177315 0.987234706347581
0.000972296276528  -4.466532872934121 0.985508231684999
0.000972993613533  -4.467757048737792 0.983782597604362
0.000974000785205  -4.468918449757659 0.982058335328984
0.000975321206273  -4.470015720425237 0.980335977889230
0.000976958826632  -4.471047403729782 0.978616059975027
0.000978918148859  -4.472011938243347 0.976899117752067
0.000981204248167  -4.472907654760871 0.975185688683845
0.000983822795018  -4.473732772516348 0.973476311357398
0.000986780080491  -4.474485395033108 0.971771525292660
0.000990083044619  -4.475163505485457 0.970071870734375
0.000993739308027  -4.475764961634439 0.968377888460546
0.000997757207001  -4.476287490185926 0.966690119559326
0.001002145832424  -4.476728680662956 0.965009105219735
0.001006915072800  -4.477085978604960 0.963335386490364
0.001012075661864  -4.477356678171161 0.961669504064674
0.001017639231092  -4.477537914015268 0.960011998006846
0.001023618367741  -4.477626652358884 0.958363407518746
0.001030026678853  -4.477619681242085 0.956724270669641
0.001036878861917  -4.477513599847341 0.955095124115890
0.001044190782928  -4.477304806759326 0.953476502823156
0.001051979562526  -4.476989487155612 0.951868939769982
0.001060263671332  -4.476563598704270 0.950272965626457
0.001069063035314  -4.476022856124226 0.948689108451285
0.001078399152550  -4.475362714215790 0.947117893336593
0.001088295222606  -4.474578349190693 0.945559842071138
0.001098776290204  -4.473664638164040 0.944015472763928
0.001109869404846  -4.472616136541590 0.942485299463856
0.001121603798543  -4.471427053083234 0.940969831750880
0.001134011083889  -4.470091222386274 0.939469574311788
0.001147125475312  -4.468602074436102 0.937985026491452
0.001160984036531  -4.466952600924484 0.936516681815986
0.001175626957915  -4.465135317885419 0.935065027485895
0.001191097867948  -4.463142224192278 0.933630543839829
0.001207444183753  -4.460964755409359 0.932213703782042
0.001224717506395  -4.458593732330876 0.930814972167297
0.001242974067723  -4.456019303561788 0.929434805130577
0.001262275236763  -4.453230881257064 0.928073649382751
0.001282688094944  -4.450217069102645 0.926731941425302
0.001304286091341  -4.446965581400276 0.925410106706889
0.001327149791125  -4.443463151990092 0.924108558694452
0.001351367732971  -4.439695431438106 0.922827697847479
0.001377037414339  -4.435646870739886 0.921567910492499
0.001404266427297  -4.431300589416395 0.920329567552576
0.001433173772413  -4.426638225469351 0.919113023142749
0.001463891383994  -4.421639764261887 0.917918612983278
0.001496565907280  -4.416283342727678 0.916746652588892
0.001531360777478  -4.410545024707535 0.915597435219628
0.001568458661971  -4.404398542185730 0.914471229529014
0.001608064341769  -4.397814996324893 0.913368276841582
0.001650408126980  -4.390762510675060 0.912288787997128
0.001695749925212  -4.383205827368949 0.911232939652124
0.001744384112993  -4.375105834947984 0.910200869936291
0.001796645400960  -4.366419013804112 0.909192673290152
0.001852915937186  -4.357096781816254 0.908208394309342
0.001913633963974  -4.347084718326792 0.907248020339510
0.001979304438591  -4.336321638956804 0.906311472489330
0.002050512157033  -4.324738486258682 0.905398594672314
0.002127938095470  -4.312256991442435 0.904509140081428
0.002212379926894  -4.298788049302101 0.903642754405397
0.002304778009748  -4.284229730930023 0.902798954763590
0.002406248626099  -4.268464834869572 0.901977103030645
0.002518126937618  -4.251357844521312 0.901176371702644
0.002642023134824  -4.232751113689719 0.900395699709655
0.002779896748912  -4.212460037414608 0.899633734516190
0.002934156350967  -4.190266871974685 0.898888755266820
0.003107795337292  -4.165912732114099 0.898158569228656
0.003304579969084  -4.139087090964472 0.897440370064432
0.003529314656287  -4.109413800794020 0.896730540332491
0.003788224095930  -4.076432174511963 0.896024370817481
0.004089516851567  -4.039570905127388 0.895315652616919
0.004444239122095  -3.998111347125019 0.894596069133892
0.004867608652992  -3.951134554796789 0.893854263105887
0.005381174970585  -3.897442714627276 0.893074356036064
0.006016468959968  -3.835438671220115 0.892233504032521
0.006821488777965  -3.762933744132360 0.891297668693953
0.007872960954398  -3.676826032918474 0.890213869991026
0.009301380227083  -3.572528644488063 0.888894953575232
0.011347505935047  -3.442872097900130 0.887186803547526
0.014508161408529  -3.275768993328194 0.884788750793815
0.019992655528039  -3.048465485954458 0.881024739323846
0.031669557869812  -2.709785321526165 0.873989668808664

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2015-01-28 14:28:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 4 个回答

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

仅供参考，运行完毕后，请到D盘找名为 shuju.xls 的excel文档

CODE:

function solve_noneqs2

clear all;clc

format long

global thita

x0=[0.8 -5 0.1];

for i=1:1:88

    thita=i*3.14/180.0;

x(1:3,i)=fsolve(@objfun,x0);

end

xlswrite('D:\shuju.xls',x)

function f=objfun(x)

global thita

gam=0.2;

gama=gam/(cos(thita)+gam*(1-cos(thita)));  

nn=0;

mie=1836.0*4.0;

miu=1.0/mie;

tie=0.01;

rou=0.2;

uiy0=-0.01;

e0=-0.01;

f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));

f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));

f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...

+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);

附图1.png

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2015-01-28 14:23:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaoxyah

铁杆木虫 (著名写手)

应助: 3 (幼儿园)
金币: 5329.7
散金: 2779
红花: 5
帖子: 1089
在线: 974.6小时
虫号: 533803
注册: 2008-03-27
专业: 等离子体物理

引用回帖:

3楼: Originally posted by 月只蓝 at 2015-01-28 14:28:26
或者更简单点，把x结果显示在主程序窗口中，复制粘贴到Origin即可：

function solve_noneqs2
clear all;clc
format long
global thita
x0=;
for i=1:1:88
thita=i*3.14/180.0;
x(i,1:3)=fsolve(@objf ...

非常感谢您的每一次指导！祝每天开心

赞一下

回复此楼

4楼2015-01-28 14:34:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 4 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 恳请有学校收留 +3	柯淮然 2026-04-12	3/150	2026-04-14 16:25 by 逆水乘风
[考研] 0854调剂 +11	长弓傲 2026-04-12	14/700	2026-04-14 16:07 by 逆水乘风
[考研] 297，工科调剂?河南农业大学本科 +9	河南农业大学-能 2026-04-14	9/450	2026-04-14 16:03 by Art1977
[考研] 求调剂 +20	MAX怅惘 2026-04-09	22/1100	2026-04-14 14:57 by 独醉梦孤城
[考研] 材料工程281还有调剂机会吗 +38	xaw. 2026-04-11	39/1950	2026-04-14 13:04 by wanganpin
[考研] 08工学 309分求调剂 +11	Yin DY 2026-04-08	11/550	2026-04-14 12:55 by fqwang
[考研] 人工智能320调剂08工类还有机会吗 +18	振—TZ 2026-04-10	19/950	2026-04-14 10:34 by screening
[考研] 332求调剂 +15	蕉蕉123 2026-04-10	15/750	2026-04-13 23:12 by pies112
[考研] 本科西工大 324求调剂 +5	wysyjs25 2026-04-10	5/250	2026-04-13 23:08 by pies112
[考研] B区0809 ，数一英一，290 求调剂 +3	泠潍1111 2026-04-12	4/200	2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿西交机械专硕求调剂 +9	求上岸的小王 2026-04-10	9/450	2026-04-13 16:08 by jiangguiquan11
[考研] 368化学求调剂 +14	wwwwabcde 2026-04-07	15/750	2026-04-13 08:36 by lhj2009
[考研] 339求调剂 +8	hanwudada 2026-04-11	9/450	2026-04-12 15:36 by laoshidan
[考研] 352 求调剂 +6	yzion 2026-04-11	8/400	2026-04-11 16:24 by 明月此时有
[考研] 生物学308求调剂（一志愿华东师大） +6	相信必会光芒万� 2026-04-10	6/300	2026-04-11 05:23 by zhuwenxu
[考研] 吉大计算机技术331分，英语六级，求调剂 +3	峰峰021116 2026-04-09	3/150	2026-04-10 20:01 by chemisry
[考研] 344求调剂 +7	丶风雪夜归人丶 2026-04-09	7/350	2026-04-10 12:05 by pengliang8036
[考研] 314求调剂 +14	weltZeng 2026-04-09	14/700	2026-04-09 23:14 by wolf97
[考研] 二次调剂求老师收留 +3	笑笑袁 2026-04-08	3/150	2026-04-08 23:50 by 醉在风里
[考研] 296求调剂 +3	汪！？！ 2026-04-08	3/150	2026-04-08 22:00 by zhouyuwinner

24小时热门版块排行榜

zhaoxyah

[求助] 请问matlab解方程组结果如何保存 已有1人参与

» 猜你喜欢

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

月只蓝

【答案】应助回帖

月只蓝

【答案】应助回帖

zhaoxyah

[求助] 请问matlab解方程组结果如何保存已有1人参与