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

登录注册

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

24小时热门版块排行榜

返回列表

woqunina

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
红花: 1
帖子: 7
在线: 1.4小时
虫号: 1517630
注册: 2011-11-30
专业: 应用地球物理学

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

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

» 猜你喜欢

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

请教怎样添加error bars 已经有9人回复
有没有从计算转向实验的同学已经有12人回复
从粗产物的氢谱怎么来算产率？已经有28人回复
【分享】东方文化经典讲读Mp3系列（四书五经类）(转)【已搜索无重复】已经有12人回复
FAO列出的五大健康食品已经有4人回复
PWSCF计算PDOS时，选哪列数据？已经有15人回复
如何计算四氯化硅水解的反应热已经有2人回复
诚求高手指导下EXCEL里的数据处理问题（把五排四列的数据排成一列）已经有5人回复
【求助】MS计算CPU占用率问题，求教已经有17人回复
重奖征集国内国外做理论计算课题组的信息已经有44人回复

1楼 2011-11-30 19:47:13

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

woqunina

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
红花: 1
帖子: 7
在线: 1.4小时
虫号: 1517630
注册: 2011-11-30
专业: 应用地球物理学

有没有会啊

回复此楼

2楼2011-12-01 19:06:49

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

wangyuewu08

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

【答案】应助回帖

★ ★
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

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

【答案】应助回帖

★ ★
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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
红花: 1
帖子: 7
在线: 1.4小时
虫号: 1517630
注册: 2011-11-30
专业: 应用地球物理学

感谢你的回复，但是不对啊，运算不出来。b = data; 是否应该改成b=a；但是再往下也不对。麻烦在看看，谢谢

赞一下

回复此楼

5楼2011-12-04 17:53:22

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

wangyuewu08

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

【答案】应助回帖

★
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

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

★
xiegangmai(金币+1): 谢谢参与！ 2011-12-04 23:45:54

最后一个块注释弄错了，百分号都要放在大括号的前面，应该是：
%{
注释内容
%}

赞一下(1人)

回复此楼

往事并不如烟

7楼2011-12-04 20:05:11

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

wangyuewu08

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

【答案】应助回帖

★
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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
红花: 1
帖子: 7
在线: 1.4小时
虫号: 1517630
注册: 2011-11-30
专业: 应用地球物理学

非常感谢楼上的大力支持
现在豁然开朗！
谢谢！
现将最后完成的程序发上。
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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
红花: 1
帖子: 7
在线: 1.4小时
虫号: 1517630
注册: 2011-11-30
专业: 应用地球物理学

送鲜花一朵

лл

回复此楼

10楼2011-12-05 10:46:21

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

相关版块跳转我要订阅楼主 woqunina 的主题更新

返回列表