| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 546 人关注过本帖
标题:求高手检查
只看楼主 加入收藏
nba123344
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2012-9-15
收藏
 问题点数:0 回复次数:0 
求高手检查
求检查    我这个程序运行不出来:
clc;
[Filename,Pathname]=uigetfile('','选择谱数据');             %输入能谱数据
fid=fopen([Pathname Filename],'r');
array_z=fscanf(fid,'%d',inf);
fclose(fid);
count=length(array_z);
x=linspace(0,count,count);
plot(x,array_z)
array_feng=[3603,1029,545];
array_dao=[703,1505,5198];
FWHM=zeros(1,3);
E=zeros(1,3);
for i=1:3
    E(i)=array_dao(i)*0.0003396;
    for j=array_dao(i)-10:array_dao(i)
        if array_z(j)>array_feng(i)/2
            Ib=array_z(j);
            Eb=0.0003396*j;
            Ia=array_z(j-1);
            Ea=0.0003396*(j-1);
        end
    end
    for k=array_dao(i):array_dao(i)+10
        if array_z(k)<array_feng(i)/2
            Id=array_z(k);
            Ed=0.0003396*k;
            Ic=array_z(k-1);
            Ec=0.0003396*(k-1);
        end
    end
    I=array_dao(i)*0.0003396;
    FWHM(i)=Ed+(I-Id)/(Ic-Id)*(Ec-Ed)-(Ea+(I-Ia)/(Ib-Ia)*(Eb-Ea));
end

L1='a+b*sqrt(E(1)+c*E(1)*E(1))=FWHM(1)';
L2='a+b*sqrt(E(2)+c*E(2)*E(2))=FWHM(2)';
L3='a+b*sqrt(E(3)+c*E(3)*E(3))=FWHM(3)';
g=solve(L1,L2,L3)



但是我检查出只是这部分程序出错误了:
L1='a+b*sqrt(E(1)+c*E(1)*E(1))=FWHM(1)';
L2='a+b*sqrt(E(2)+c*E(2)*E(2))=FWHM(2)';
L3='a+b*sqrt(E(3)+c*E(3)*E(3))=FWHM(3)';
g=solve(L1,L2,L3)


而我把它换成:
L1='a+b*sqrt(0.2387+c*0.2387*0.2387)=0.0017 ';
L2='a+b*sqrt(0.5111+c*0.5111*0.5111)=0.0101';
L3='a+b*sqrt(1.7652+c*1.7652*1.7652)=0.0076';
g=solve(L1,L2,L3)
却算的出来          两者是等价的算出来检查过      



[ 本帖最后由 nba123344 于 2012-9-15 09:49 编辑 ]
搜索更多相关主题的帖子: 检查 count 
2012-09-15 09:46
快速回复:求高手检查
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.023163 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved