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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？
51/1返回列表
查看: 2903  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cw277

木虫 (正式写手)

[求助] 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？

当i=600时，结果可以输出。
但当i>=700  这个char就卡住了，可能是没办法转换了。。
求高手指教啊。。  我应该怎么办？？

程序如下

clc
clear
tic;

S(1)=sym('k1');
S(2)=sym('k2');
S(3)=sym('k3');
S(4)=sym('k4');
S(5)=sym('k5');
S(6)=sym('k6');
S(7)=sym('k7');
S(8)=sym('k8');
S(9)=sym('k9');
S(10)=sym('k10');
S(11)=sym('k11');
S(12)=sym('k12');
S(13)=sym('k13');
S(14)=sym('k14');
S(15)=sym('k15');
S(16)=sym('k16');
S(17)=sym('k17');
S(18)=sym('k18');
S(19)=sym('k19');
S(20)=sym('k20');
S(21)=sym('k21');
S(22)=sym('k22');
S(23)=sym('k23');
S(24)=sym('k24');
S(25)=sym('k25');
S(26)=sym('k26');
S(27)=sym('k27');
S(28)=sym('k28');
S(29)=sym('k29');
S(30)=sym('k30');
S(31)=sym('k31');
S(32)=sym('k32');
S(33)=sym('k33');
S(34)=sym('k34');
S(35)=sym('k35');
S(36)=sym('k36');
S(37)=sym('k37');
S(38)=sym('k38');
S(39)=sym('k39');
S(40)=sym('k40');
S(41)=sym('k41');
S(42)=sym('k42');
S(43)=sym('k43');
S(44)=sym('k44');
S(45)=sym('k45');
S(46)=sym('k46');
S(47)=sym('k47');
S(48)=sym('k48');
S(49)=sym('k49');
S(50)=sym('k50');
S(51)=sym('k51');
S(52)=sym('k52');
S(53)=sym('k53');
S(54)=sym('k54');
S(55)=sym('k55');
S(56)=sym('k56');
S(57)=sym('k57');
S(58)=sym('k58');
S(59)=sym('k59');
S(60)=sym('k60');
S(61)=sym('k61');
S(62)=sym('k62');
S(63)=sym('k63');
S(64)=sym('k64');
S(65)=sym('k65');
S(66)=sym('k66');
S(67)=sym('k67');
S(68)=sym('k68');
S(69)=sym('k69');
S(70)=sym('k70');
S(71)=sym('k71');
S(72)=sym('k72');
S(73)=sym('k73');
S(74)=sym('k74');
S(75)=sym('k75');
S(76)=sym('k76');
S(77)=sym('k77');
S(78)=sym('k78');
S(79)=sym('k79');
S(80)=sym('k80');
S(81)=sym('0');
S(82)=sym('0');
S(83)=sym('0');
S(84)=sym('0');
S(85)=sym('0');
S(86)=sym('0');
S(87)=sym('0');
S(88)=sym('0');
S(89)=sym('0');
S(90)=sym('0');
S(91)=sym('0');
S(92)=sym('0');
S(93)=sym('0');
S(94)=sym('IV1');
S(95)=sym('IV2');
S(96)=sym('IV3');
S(97)=sym('IV4');
S(98)=sym('IV5');
S(99)=sym('IV6');
S(100)=sym('IV7');
S(101)=sym('IV8');
S(102)=sym('IV9');
S(103)=sym('IV10');
S(104)=sym('IV11');
S(105)=sym('IV12');
S(106)=sym('IV13');
S(107)=sym('IV14');
S(108)=sym('IV15');
S(109)=sym('IV16');
S(110)=sym('IV17');
S(111)=sym('IV18');
S(112)=sym('IV19');
S(113)=sym('IV20');
S(114)=sym('IV21');
S(115)=sym('IV22');
S(116)=sym('IV23');
S(117)=sym('IV24');
S(118)=sym('IV25');
S(119)=sym('IV26');
S(120)=sym('IV27');
S(121)=sym('IV28');
S(122)=sym('IV29');
S(123)=sym('IV30');
S(124)=sym('IV31');
S(125)=sym('IV32');
S(126)=sym('IV33');
S(127)=sym('IV34');
S(128)=sym('IV35');
S(129)=sym('IV36');
S(130)=sym('IV37');
S(131)=sym('IV38');
S(132)=sym('IV39');
S(133)=sym('IV40');
S(134)=sym('IV41');
S(135)=sym('IV42');
S(136)=sym('IV43');
S(137)=sym('IV44');
S(138)=sym('IV45');
S(139)=sym('IV46');
S(140)=sym('IV47');
S(141)=sym('IV48');
S(142)=sym('IV49');
S(143)=sym('IV50');
S(144)=sym('IV51');
S(145)=sym('IV52');
S(146)=sym('IV53');
S(147)=sym('IV54');
S(148)=sym('IV55');
S(149)=sym('IV56');
S(150)=sym('IV57');
S(151)=sym('IV58');
S(152)=sym('IV59');
S(153)=sym('IV60');
S(154)=sym('IV61');
S(155)=sym('IV62');
S(156)=sym('IV63');
S(157)=sym('IV64');
S(158)=sym('IV65');
S(159)=sym('IV66');
S(160)=sym('IV67');
S(161)=sym('IV68');
S(162)=sym('IV69');
S(163)=sym('IV70');
S(164)=sym('IV71');
S(165)=sym('IV72');
S(166)=sym('IV73');
S(167)=sym('IV74');
S(168)=sym('IV75');
S(169)=sym('IV76');
S(170)=sym('IV77');
S(171)=sym('IV78');
S(172)=sym('IV79');
S(173)=sym('IV80');
S(174)=sym('0');
S(175)=sym('0');
S(176)=sym('0');
S(177)=sym('0');
S(178)=sym('0');
S(179)=sym('0');
S(180)=sym('0');
S(181)=sym('0');
S(182)=sym('0');
S(183)=sym('0');
S(184)=sym('0');
S(185)=sym('0');
S(186)=sym('0');
S(187)=sym('0');
S(188)=sym('0');
S(189)=sym('0');
S(190)=sym('0');
S(191)=sym('0');
S(192)=sym('0');
S(193)=sym('0');
S(194)=sym('0');
S(195)=sym('0');
S(196)=sym('0');
S(197)=sym('0');
S(198)=sym('0');
S(199)=sym('0');
S(200)=sym('0');
S(201)=sym('0');
S(202)=sym('0');
S(203)=sym('0');
S(204)=sym('0');
S(205)=sym('0');
S(206)=sym('0');
S(207)=sym('0');
S(208)=sym('0');
S(209)=sym('0');
S(210)=sym('0');
S(211)=sym('0');
S(212)=sym('0');
S(213)=sym('0');
S(214)=sym('0');
S(215)=sym('0');
S(216)=sym('0');
S(217)=sym('0');
S(218)=sym('0');
S(219)=sym('0');
S(220)=sym('0');
S(221)=sym('0');
S(222)=sym('0');
S(223)=sym('0');
S(224)=sym('0');
S(225)=sym('0');
S(226)=sym('0');
S(227)=sym('0');
S(228)=sym('0');
S(229)=sym('0');
S(230)=sym('0');
S(231)=sym('0');
S(232)=sym('0');
S(233)=sym('0');
S(234)=sym('0');
S(235)=sym('0');
S(236)=sym('0');
S(237)=sym('0');
S(238)=sym('0');
S(239)=sym('0');
S(240)=sym('0');
S(241)=sym('0');
S(242)=sym('0');
S(243)=sym('0');
S(244)=sym('0');
S(245)=sym('0');
S(246)=sym('0');
S(247)=sym('0');
S(248)=sym('0');
S(249)=sym('0');
S(250)=sym('0');
S(251)=sym('0');
S(252)=sym('0');
S(253)=sym('0');
S(254)=sym('0');
S(255)=sym('0');
S(256)=sym('0');
S(257)=sym('0');
S(258)=sym('0');
S(259)=sym('0');
S(260)=sym('0');
S(261)=sym('0');
S(262)=sym('0');
S(263)=sym('0');
S(264)=sym('0');
S(265)=sym('0');
S(266)=sym('0');
S(267)=sym('0');
S(268)=sym('0');
S(269)=sym('0');
S(270)=sym('0');
S(271)=sym('0');
S(272)=sym('0');
S(273)=sym('0');
S(274)=sym('0');
S(275)=sym('0');
S(276)=sym('0');
S(277)=sym('0');
S(278)=sym('0');
S(279)=sym('0');
S(280)=sym('0');
S(281)=sym('0');
S(282)=sym('0');
S(283)=sym('0');
S(284)=sym('0');
S(285)=sym('0');
S(286)=sym('1');
S(287)=sym('1');
S(288)=sym('1');

for i=1:700
t1(i)=S(66)+S(91)*S(92)+S(93)+S(171);
t2(i)=S(162)+S(175)*S(176)+S(177)+S(264);
t3(i)=S(243)+S(286)*S(287)+S(288)+S(69);
S=[t3(i),S(1:92),S(94:288)];  
S=[S(1:93),t1(i),S(94:176),S(178:288)];  
S=[S(1:177),t2(i),S(178:287)];  
i
end

s=char(S(1));
fid=fopen('c:\char2.txt','w');
fprintf(fid,'%s',s);
toc;
clear S;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
感谢生活让我慢慢学会控制自己的情绪，学会包容与忍耐，学会不再去期待
1楼 2011-12-05 19:52:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★ ★
cw277(金币+8): 谢谢高手指点~~ 多给点 2011-12-07 23:19:58
cw277: 回帖置顶 2011-12-07 23:20:05
xzhdty(金币+2): 欢迎常来程序语言看看 2011-12-08 16:22:09
在我的电脑上运行出结果了，楼主程序的结果存放在附件char2.txt中。
我修改之后，运行的结果存放在char2_version_2.txt中。
修改后的程序是program_version_2.txt。

针对matlab的提示，我对楼主的程序进行了修改，并对运行时间进行了比较。因为运行时间比较长，只比较了一次，所以时间比较的结果可能不准确。

我主要简化了S，t1，t2，t3的初始化（楼主对S的初始化太费时费力，还容易出错），其它部分与楼主源程序一样，除了输出路径改为E（我不太习惯在C盘放文件）。
CODE:
%   修改之前的运行时间
%   Elapsed time is 427.892359 seconds.
%   修改之后的运行时间
%   Elapsed time is 333.992927 seconds.
clc
clear
tic;

S = zeros(1, 288);
S = sym(S);
for ii = 1 : 80
    S(ii) = sym(['k', num2str(ii)]);
end
%   S(81:93) = sym('0');
for ii = 94 : 173
    S(ii) = sym(['IV', num2str(ii - 93)]);
end
%   S(174 : 285) = sym('0');
S( 286 : 288 ) = sym('1');

t1 = zeros( 700, 1 );
t1 = sym(t1);
t2 = t1;
t3 = t1;

for i=1:700
t1(i)=S(66)+S(91)*S(92)+S(93)+S(171);
t2(i)=S(162)+S(175)*S(176)+S(177)+S(264);
t3(i)=S(243)+S(286)*S(287)+S(288)+S(69);
S=[t3(i),S(1:92),S(94:288)];  
S=[S(1:93),t1(i),S(94:176),S(178:288)];  
S=[S(1:177),t2(i),S(178:287)];  
i
end

s=char(S(1));
fid=fopen('E:\char2.txt','w');
fprintf(fid,'%s',s);
toc;
clear S;
赞一下(1人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : program_version_2.txt
    • 2011-12-07 19:45:46, 796 bytes
  • 附件 2 : char2.txt
    • 2011-12-07 19:46:05, 3.27 M
  • 附件 3 : char2_version_2.txt
    • 2011-12-07 19:46:27, 3.27 M
往事并不如烟
9楼2011-12-07 19:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

lshok

铁杆木虫 (著名写手)
★
jjdg(金币+1): 感谢参与 2011-12-06 02:00:03
不太明白，不过你代码最后没有关闭文件
赞一下(1人) 回复此楼
2楼2011-12-05 22:44:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★
cw277(金币+1): 这个可以考虑。。 其实主要原因是在 s=char(S(1)); 这一句。 S(1)的表达式太多了，以致于时间花很久 2011-12-06 22:57:04
jjdg(金币+1): 感谢参与 2011-12-07 01:54:13
在PC上，二进制文件和格式化文件的打开是有区别的，打开文件的语句中需要指定文件为文本格式
fid=fopen('c:\char2.txt','wt');
楼主试试，看错误是不是出在这里。
给楼主提个建议，如果程序的其它地方不需要用t1,t2,t3，那么在循环中将t1(i),t2(i),t3(i)改为t1,t2,t3，这样可以节省大量的空间和时间。
如果要用到，那么在循环之前预定义一下比较好。
赞一下(1人) 回复此楼
往事并不如烟
3楼2011-12-06 12:20:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★
cw277(金币+1): 其实主要原因是在 s=char(S(1)); 这一句。 S(1)的表达式太多了，以致于时间花很久 2011-12-06 22:57:43
jjdg(金币+1): 感谢参与 2011-12-07 01:54:06
如果是因为太长而无法写入的话，系统肯定会有提示的
赞一下(1人) 回复此楼
往事并不如烟
4楼2011-12-06 19:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +18 rare12345 2026-03-18 18/900 2026-03-19 22:36 by 学员8dgXkO
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 0703化学调剂，求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
关闭
请填处理意见
关闭