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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 请问matlab解方程组结果如何保存
41/1返回列表
查看: 902  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhaoxyah

铁杆木虫 (著名写手)

[求助] 请问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);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-28 13:50:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +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);

请问matlab解方程组结果如何保存
附图1.png
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2015-01-28 14:23:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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的回帖

zhaoxyah

铁杆木虫 (著名写手)
引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 zhaoxyah 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭