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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】函数使用问题
31/1返回列表
查看: 334  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jacketdanster

金虫 (小有名气)

[交流] 【求助】函数使用问题 已有2人参与

原来在6.0下面可以运行的函数，装了7.0a以后运行不了，请帮忙分析一下，谢谢

提示出错如下：

??? Error using ==> fopen
Invalid filename.

Error in ==> pretreatment at 6
fid=fopen([pname fname],'r');

我运行的函数如下：
clc,clear
format short e

% Open a File and Import Original Data.
[fname,pname]=uigetfile('*.txt','Open a file:');
fid=fopen([pname fname],'r');
fgetl(fid);% Here 'fgetl' command is used to jump to read next line.
fgetl(fid);
fgetl(fid);
a=fscanf(fid,'%f',[7 inf]);
fclose(fid);
a1=a';
Dtime=a1(:,1);Ddtg=a1(:,2);Dtg=a1(:,3);Ddhf=a1(:,4);Dhf=a1(:,5);Ddt=a1(:,6);Dt=a1(:,7);

%hl1=line(Dt,Dtg,'color','k');
%ax1=gca;
%set(ax1,'Xcolor','k','ycolor','k')
%ax2=axes('position',get(ax1,'Position'), ...
%    'XAxisLocation','top',...
%    'YAxisLocation','right',...
%    'Color','none',...
%    'XColor','b','YColor','b');
%hl2=line(Dt,Ddtg,'Color','b','Parent',ax2);
%grid on;
[haxes,hline1,hline2]=plotyy(Dt,Dtg,Dt,Ddtg);grid on;
xlabel('Temperature ( ^oC )')
axes(haxes(1))
ylabel('Thermogravimetry (%)')
axes(haxes(2))
ylabel('DTG (d%)')



% Select the Temperature Range for Calculation
prompt={'Enter the Initial Temperature /C:','Enter the Final Temperature /C:'};
def={'20','200'};
num_line=1;
answer=inputdlg(prompt,'Please Input the Temperature Range for Calculate!',num_line,def);
T(1)=str2num([answer{1}]);
T(2)=str2num([answer{2}]);


r=1;
while Dt(r)     r=r+1;
end
    rs=r;
while Dt(r)     r=r+1;
end
    re=r-1;
Dalfa=(Dtg(rs:re)-Dtg(rs))/-13.73;%%%%%%%%%
Dt=Dt(rs:re)+273.15;
Dtime=Dtime(rs:re);
Ddt=Ddt(rs:re);
Ddalfa=Ddtg(rs:re)/-13.73;

% Program for interpolation and calculate Ddtg & Ddt
rn=1;
for i=0.10:0.02:0.75 %%%%%%
    Dalfa1(rn,1)=i;
    r=1;
    ri=0;
    rf=0;
    n=size(Dalfa,1);
    while Dalfa(r)         r=r+1;
    end
    rt=r;
    if rt<12
        ri=1;
        rf=rt+11;
    elseif rt>n-11
        ri=rt-12;
        rf=n;
    else
        ri=rt-12;
        rf=rt+11;
    end
   
    % calculate for interpolation of temperature
    %p=polyfit(Dt(ri:rf),Dalfa(ri:rf),2)
    p=[Dt(ri:rf).^2 Dt(ri:rf) ones(rf-ri+1,1)]\Dalfa(ri:rf);
    %p1=polyval(p,Dt(ri:rf));
    p(3)=p(3)-Dalfa1(rn,1);
    DtI=roots(p);
    if abs(DtI(1)-Dt(rt))<=abs(DtI(2)-Dt(rt))%DtI(1)>=Dt(rt-2) & DtI(1)<=Dt(rt+1)
        Dt_new(rn,1)=DtI(1);
    else
        Dt_new(rn,1)=DtI(2);
    end
   
    % Calculate the differential temperature
    pat=[Dtime(ri:rf).^2 Dtime(ri:rf) ones(rf-ri+1,1)]\Dalfa(ri:rf);
    pat(3)=pat(3)-Dalfa1(rn,1);
    DtimeI=roots(pat);
    if abs(DtimeI(1)-Dtime(rt))<=abs(DtimeI(2)-Dtime(rt))
        Dtime_new(rn,1)=DtimeI(1);
    else
        Dtime_new(rn,1)=DtimeI(2);
    end
    Ddtg_new(rn,1)=(2*pat(1)*Dtime_new(rn,1)+pat(2))*60;
   
   
    % calculate for interpolation of time
    pt=[Dtime(ri:rf).^2 Dtime(ri:rf) ones(rf-ri+1,1)]\Ddt(ri:rf);
   
    % Calculate the differential temperature
    Ddt_new(rn,1)=polyval(pt,Dtime_new(rn,1));
    %Ddt_new(rn,1)=(2*pt(1)*Dtime_new(rn,1)+pt(2))*60;
    rn=rn+1;
    rt=0;
end

% Export the Result
c=[Dalfa1 Dt_new Ddtg_new Ddt_new];
[filename,pathname]=uiputfile('*.txt','Save the Result File Name');
fid=fopen([pathname filename],'wt');
fprintf(fid,'Temperature Range: %6.2f - %6.2f C\n',T);
fprintf(fid,' alfa       T(K)    dalfa/dt(min-1)   Rate\n');
fprintf(fid,'%5.4f     %8.4f     %8.6f    %8.4f\n',c');
fclose(fid);
clf
figure(1)
plot(Dt,Dalfa,Dt_new,Dalfa1,'o',Dt,Ddalfa,Dt_new,Ddtg_new,'^')
回复此楼

» 猜你喜欢

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

anyuezhiji

银虫 (正式写手)

星空行者
★ ★
jacketdanster(金币+1):谢谢参与
nono2009(金币+1):辛苦了。 2010-04-23 15:22
引用回帖:
clc,clear
format short e

% Open a File and Import Original Data.
[fname,pname]=uigetfile('*.txt','Open a file:');
fid=fopen([pname fname],'r');

报错报的是这一行
fid=fopen([pname fname],'r');

但是我测试前面的代码都没有问题

打开uigetfile函数会让你选择文件 是不是选了
选的文件的文件名最好是英文 且不含一些特殊字符
赞一下(2人) 回复此楼
暗月下没有留下风的痕迹，但它已经寂然飘逝。。By&amp;amp;lt;暗月之寂&amp;amp;gt;:tiger38:
2楼2010-04-23 14:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jacketdanster

金虫 (小有名气)
引用回帖:
Originally posted by anyuezhiji at 2010-04-23 14:52:30:


报错报的是这一行
fid=fopen([pname fname],'r');

但是我测试前面的代码都没有问题

打开uigetfile函数会让你选择文件 是不是选了
选的文件的文件名最好是英文 且不含一些特殊字符

谢谢
回复此楼
3楼2010-04-23 15:42:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jacketdanster 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭