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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 如何算第四、五列
131/2返回列表12下一页
查看: 844  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

woqunina

新虫 (初入文坛)

[求助] 如何算第四、五列

94372.000     74670.000           0.000        1491.000            1491.000
          94372.000     74670.000          50.000        1493.000            1496.000
          94372.000     74670.000         100.000        1494.000            1793.000
          94372.000     74670.000         150.000        1497.000            1792.000
          94372.000     74670.000         200.000        1491.000            1795.000
          94372.000     74670.000         250.000        1497.000            1793.000
          94372.000     74670.000         300.000        1496.000            1794.000
          55688.500     65415.340           0.000        1571.000            1553.000
          55688.555     65415.340          50.000        1573.000            1553.000
          55688.555     65415.340         100.000        1574.000            1763.000
          55688.555     65415.340         150.000        1577.000            1773.000
          55688.555     65415.340         200.000        1571.000            1763.000
          55688.555     65415.340         250.000        1577.000            1776.000
          55688.555     65415.340         300.000        1576.000            1774.000
          81677.912     81412.098           0.000        1571.000            1553.000
          81677.912     81412.098          50.000        1573.000            1553.000
          81677.912     81412.098         100.000        1574.000            1763.000
          81677.912     81412.098         150.000        1577.000            1773.000
          81677.912     81412.098         200.000        1571.000            1763.000
          81677.912     81412.098         250.000        1577.000            1776.000
          81677.912     81412.098         300.000        1576.000            1774.000
              .
              .
              .
一个数据文件“A.txt” 一共有五列，几千行。所有数据每7行为一组，即（1-7行）具有相同的第一列，第二列。
（8-14行）具有相同的第一列，第二列。。。。
第三列从0以50递增直到第7行的300。然后从第8行的0 增加到第14行300 。。。         
第四列和第五列是两个实验值（无规律）。
  
  我想通过matlab运算
  输出一个文件“B.txt”，他的格式如下：
                 
  第一列二列,三列不变，     第四列的计算公式如下                                   第五列数据
                                              【第四列第n行记做" B(n,4)" 】                                   
                                                B(1,4)=A(1,4)                                                     B(1,5)=A(1,4)/A(1,5)
                                                B(2,4)=[50*A(2,4)+  0*B(1,4)]/50                         B(2,5)=A(2,4)/A(2,5)
                                                B(3,4)=[50*A(3,4)+ 50*B(2,4)]/100                       B(3,5)=A(3,4)/A(3,5)
                                               B(4,4)=[50*A(4,4)+100*B(3,4)]/150                        。
                                               B(5,4)=[50*A(5,4)+150*B(4,4)]/200                         。
                                               B(6,4)=[50*A(6,4)+200*B(5,4)]/250                          。
                                              B(7,4)=[50*A(7,4)+250*B(6,4)]/300                          。
                                              B(8,4)= A(8,4)                                                         。
                                              B(9,4)= [50*A(9,4)+   0*B(8,4)]/50               
                                             B(10,4)=[50*A(10,4)+ 50*B(9,4)]/100              
                                             B(11,4)=[50*A(11,4)+100*B(10,4)]/150            
                                             B(12,4)=[50*A(12,4)+150*B(11,4)]/200            
                                                B(13,4)=[50*A(13,4)+200*B(12,4)]/250            
                                             B(14,4)=[50*A(14,4)+250*B(13,4)]/300            
                                               B(15,4)= A(15,4)                                          
                                               。                                                            
                                               。                                                            
                                              。                                                            
                                            以此类题                                                      
                                                                                                   
其中第五列数据，判断的第n行数据B(n,5)是否介于0.7-1.6之间， 如果是则输出实际值。
如果 B(n,5）<0.7,那么 B(n,5)=A（n,5)/4+0.525,   如果B(n,5）>1.6,那么B(n,5)=, A(n,5)/4+1.2

最后输出   “ B.txt   ”




clc;
clear;
dat=dlmread('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[M,N]=size(dat);
for i=1:M

  if dat(i,3)=0  
    dat(1,4)= dat(1,4)
     else ????
   
       fprintf(fp,'%20.2f   %20.2f  %20.2f  %20.2f   %20.2f \n',dat(i,1),dat(i,2),dat(i,3),dat(i,4),dat(i,5));
        end
    end
end
fclose(fp);

第四列和第五列我不知道该怎么办了，请大家帮忙给看看，谢谢

[ Last edited by woqunina on 2011-11-30 at 21:24 ]
回复此楼

» 本帖已获得的红花（最新10朵）

woqunina

» 猜你喜欢

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

woqunina

新虫 (初入文坛)
有没有会啊
回复此楼
2楼2011-12-01 19:06:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-12-04 18:58:41
我不知道怎么在回帖中粘贴代码，就直接复制过来了
a表示A，b表示B，建议楼主少用大写字母，大写字母一般都是留给外部变量使用的。我不是很清除格式化输出，感觉楼主比我更了解
clc
clear
a = dlmread('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[m, n]=size( a );
b = data;   %因为有三列值不用改变，可以将原始数据直接赋给b
for i = 1 : 7 : m
     if dat(i, 3) <= 10^(-5)
         b( i,   4 ) = a( i, 4 );
         b( i+1, 4 ) = ( 50 * a( i+1, 4 ) +    0  * b( i,     4 ) ) / 50;
         b( i+2, 4 ) = ( 50 * a( i+2, 4 ) +  50  * b( i+1, 4 ) ) / 50;
         b( i+3, 4 ) = ( 50 * a( i+3, 4 ) + 100 * b( i+2, 4 ) ) / 50;
         b( i+4, 4 ) = ( 50 * a( i+4, 4 ) + 150 * b( i+3, 4 ) ) / 50;
         b( i+5, 4 ) = ( 50 * a( i+5, 4 ) + 200 * b( i+4, 4 ) ) / 50;
         b( i+6, 4 ) = ( 50 * a( i+6, 4 ) + 250 * b( i+5, 4 ) ) / 50;
         b( i,     5 ) = a( 1,    4 ) / a( 1, 5 );
         b( i+1, 5 ) = a( i+1, 4 ) / a( i+1, 5 );
         b( i+2, 5 ) = a( i+2, 4 ) / a( i+2, 5 );
         b( i+3, 5 ) = a( i+3, 4 ) / a( i+3, 5 );
         b( i+4, 5 ) = a( i+4, 4 ) / a( i+4, 5 );
         b( i+5, 5 ) = a( i+5, 4 ) / a( i+5, 5 );
         b( i+6, 5 ) = a( i+6, 4 ) / a( i+6, 5 );
     end
end

%  判断B( n ,5 )
for i = 1 : m
     if b( i, 5 ) < 0.7
         b( i, 5 ) = a( i, 5 ) / 4 + 0.525;
     else b( i, 5 ) > 1.6
         b( i, 5 ) = a( i, 5 ) / 4 + 1.2;
     end
end
赞一下(1人) 回复此楼
往事并不如烟
3楼2011-12-03 20:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★ ★
xiegangmai(金币+2): 辛苦了！ 2011-12-04 18:58:51
修改一下，效率提高了一些
CODE:
clc
clear
a = dlmread('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[m, n]=size( a );
b = data;   %因为有三列值不用改变，可以将原始数据直接赋给b
for i = 1 : 7 : m
     if dat(i, 3) <= 10^(-5)
         b( i,   4 ) = a( i, 4 );
         b( i+1, 4 ) = ( 50 * a( i+1, 4 ) +    0  * b( i,     4 ) ) / 50;
         b( i+2, 4 ) = ( 50 * a( i+2, 4 ) +  50  * b( i+1, 4 ) ) / 50;
         b( i+3, 4 ) = ( 50 * a( i+3, 4 ) + 100 * b( i+2, 4 ) ) / 50;
         b( i+4, 4 ) = ( 50 * a( i+4, 4 ) + 150 * b( i+3, 4 ) ) / 50;
         b( i+5, 4 ) = ( 50 * a( i+5, 4 ) + 200 * b( i+4, 4 ) ) / 50;
         b( i+6, 4 ) = ( 50 * a( i+6, 4 ) + 250 * b( i+5, 4 ) ) / 50;
      end
end

%  用matlab的数组运算计算B的第5列
b( :, 5 ) = a( :, 4 ) ./ a( :, 5 );

%  判断B( n ,5 )
for i = 1 : m
     if b( i, 5 ) < 0.7
         b( i, 5 ) = a( i, 5 ) / 4 + 0.525;
     else b( i, 5 ) > 1.6
         b( i, 5 ) = a( i, 5 ) / 4 + 1.2;
     end
end

赞一下(1人) 回复此楼
往事并不如烟
4楼2011-12-04 12:12:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woqunina

新虫 (初入文坛)
感谢你的回复，但是不对啊，运算不出来。b = data; 是否应该改成b=a；但是再往下也不对。麻烦在看看，谢谢
赞一下 回复此楼
5楼2011-12-04 17:53:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★
xiegangmai(金币+1): 谢谢参与！ 2011-12-04 23:45:39
引用回帖:
5楼: Originally posted by woqunina at 2011-12-04 17:53:22:
感谢你的回复，但是不对啊，运算不出来。b = data; 是否应该改成b=a；但是再往下也不对。麻烦在看看，谢谢

的确，不好意思弄错了，还有第一个 if 语句中的 dat 也要改为 a
我又修改并检查了一遍，逻辑上应该没有错误了。
如果用的是C语言，就不会出现这种名称弄错的问题了
CODE:
clc
clear;
a = dlmread('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[m, n]=size( a );
b = a;   %因为第1,2,3列值不用改变，可以将原始数据直接赋给b
for i = 1 : 7 : m
     if a(i, 3) <= 10^(-5)
         b( i  , 4 ) = a( i, 4 );
         b( i+1, 4 ) = ( 50 * a( i+1, 4 ) + 0   * b( i  , 4 ) ) / 50;
         b( i+2, 4 ) = ( 50 * a( i+2, 4 ) + 50  * b( i+1, 4 ) ) / 50;
         b( i+3, 4 ) = ( 50 * a( i+3, 4 ) + 100 * b( i+2, 4 ) ) / 50;
         b( i+4, 4 ) = ( 50 * a( i+4, 4 ) + 150 * b( i+3, 4 ) ) / 50;
         b( i+5, 4 ) = ( 50 * a( i+5, 4 ) + 200 * b( i+4, 4 ) ) / 50;
         b( i+6, 4 ) = ( 50 * a( i+6, 4 ) + 250 * b( i+5, 4 ) ) / 50;
     end
end

%  用matlab的数组运算计算B的第5列
b( :, 5 ) = a( :, 4 ) ./ a( :, 5 );

%  判断并修改B的第5列
index = b( :, 5 ) < 0.7;
b( index ) = a( index ) / 4 + 0.525;
index = b( :, 5 ) > 1.6;
b( index ) = a( index ) / 4 + 1.2;

%{尽量不用循环
for i = 1 : m
     if b( i, 5 ) < 0.7
         b( i, 5 ) = a( i, 5 ) / 4 + 0.525;
     elseif b( i, 5 ) > 1.6
         b( i, 5 ) = a( i, 5 ) / 4 + 1.2;
     end
end
}%
赞一下(1人) 回复此楼
往事并不如烟
6楼2011-12-04 20:01:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)
★
xiegangmai(金币+1): 谢谢参与！ 2011-12-04 23:45:54
最后一个块注释弄错了，百分号都要放在大括号的前面，应该是：
%{
注释内容
%}
赞一下(1人) 回复此楼
往事并不如烟
7楼2011-12-04 20:05:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

【答案】应助回帖

★
xiegangmai(金币+1): 谢谢参与！ 2011-12-04 23:46:04
楼主能够把输入数据和正确的输出数据给我吗？这样可以方便我调试程序，
第4列的求解公式的最后一个数字看错了，没有注意到它们按50递增
把第4列的求解公式尽可能化简了，以提高效率。
CODE:
clc
clear
a = dlmread('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[m, n]=size( a );
b = a;   %因为第1,2,3列值不用改变，可以将原始数据直接赋给b，简化后面的计算
for i = 1 : 7 : m
     %  b( i  , 4 ) = a( i, 4 );
     %  b( i+1, 4 ) = ( 50 * a( i+1, 4 ) + 0   * b( i  , 4 ) ) /  50;
     
     %  b( i+2, 4 ) = ( 50 * a( i+2, 4 ) + 50  * b( i+1, 4 ) ) / 100;
     b( i+2, 4 ) = 0.5 * ( a( i+2, 4 ) + b( i+1, 4 ) );
     
     b( i+3, 4 ) = ( 50 * a( i+3, 4 ) + 100 * b( i+2, 4 ) ) / 150;
     
     %  b( i+4, 4 ) = ( 50 * a( i+4, 4 ) + 150 * b( i+3, 4 ) ) / 200;
     b( i+4, 4 ) = 0.25 * a( i+4, 4 ) + 0.75 * b( i+3, 4 );
     
     %  b( i+5, 4 ) = ( 50 * a( i+5, 4 ) + 200 * b( i+4, 4 ) ) / 250;
     b( i+5, 4 ) = 0.2 * a( i+5, 4 ) + 0.8 * b( i+4, 4 );
     
     b( i+6, 4 ) = ( 50 * a( i+6, 4 ) + 250 * b( i+5, 4 ) ) / 300;
end

%  用matlab的数组运算计算B的第5列
b( :, 5 ) = a( :, 4 ) ./ a( :, 5 );

%  判断并修改B的第5列
index = b( :, 5 ) < 0.7;
b( index ) = a( index ) / 4 + 0.525;
index = b( :, 5 ) > 1.6;
b( index ) = a( index ) / 4 + 1.2;
赞一下(1人) 回复此楼
往事并不如烟
8楼2011-12-04 20:29:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woqunina

新虫 (初入文坛)
非常感谢楼上的大力支持
现在豁然开朗！
谢谢！
现将最后完成的程序发上。
clc
clear
a = load('D:\A.dat');
filename='D:\B.txt';
fp = fopen(filename,'w');
[m, n]=size( a )
b = a;
for i = 1 : 7 : m
     %  b( i  , 4 ) = a( i, 4 );
     %  b( i+1, 4 ) = ( 50 * a( i+1, 4 ) + 0   * b( i  , 4 ) ) /  50;
     
     %  b( i+2, 4 ) = ( 50 * a( i+2, 4 ) + 50  * b( i+1, 4 ) ) / 100;
     b( i+2, 4 ) = 0.5 * ( a( i+2, 4 ) + b( i+1, 4 ) );
     
     b( i+3, 4 ) = ( 50 * a( i+3, 4 ) + 100 * b( i+2, 4 ) ) / 150;
     
     %  b( i+4, 4 ) = ( 50 * a( i+4, 4 ) + 150 * b( i+3, 4 ) ) / 200;
     b( i+4, 4 ) = 0.25 * a( i+4, 4 ) + 0.75 * b( i+3, 4 );
     
     %  b( i+5, 4 ) = ( 50 * a( i+5, 4 ) + 200 * b( i+4, 4 ) ) / 250;
     b( i+5, 4 ) = 0.2 * a( i+5, 4 ) + 0.8 * b( i+4, 4 );
     
     b( i+6, 4 ) = ( 50 * a( i+6, 4 ) + 250 * b( i+5, 4 ) ) / 300;
end

%  用matlab的数组运算计算B的第5列
b( :, 5 ) = a( :, 4 ) ./ a( :, 5 );

%  判断并修改B的第5列
index = b( :, 5 ) < 0.7;
b( index ) = a( index ) / 4 + 0.525;
index = b( :, 5 ) > 1.6;
b( index ) = a( index ) / 4 + 1.2;
for i = 1 : m
    fprintf(fp,'%20.2f   %20.2f  %20.2f  %20.4f   %20.4f \n',b(i,1),b(i,2),b(i,3),b(i,4),b(i,5));
end
      
fclose(fp);
赞一下 回复此楼
9楼2011-12-05 10:44:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woqunina

新虫 (初入文坛)
送鲜花一朵
лл
回复此楼
10楼2011-12-05 10:46:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 woqunina 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 8rmuugja8q 2026-02-22 7/350 2026-02-23 09:44 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:35 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 usprnugpzw 2026-02-21 11/550 2026-02-23 09:24 by w4l55oybr1
[教师之家] 为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 8/400 2026-02-23 09:19 by 欢乐颂叶蓁
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 5/250 2026-02-23 08:04 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 6/300 2026-02-23 07:59 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +9 瞬息宇宙 2026-02-20 12/600 2026-02-21 10:39 by 欢乐颂叶蓁
信息提示
关闭
请填处理意见
关闭