| 查看: 740 | 回复: 4 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
heio金虫 (小有名气)
|
[求助]
matlab求助
|
||||||
|
本人打算用matlab处理数据中的趋势项问题 但无法运行 跪求高手指点 程序如下 %滑动平均法消除趋势项 clear clc close all hidden fni=input('输入数据文件名称','s'); fid=fopen(fni,'r') sf=fscanf(fid,'%f',1); %读入采样频率值 l=fscanf(fid,'%d',1); %读入滑动阶次 m=fscanf(fid,'%d',1); %读入平滑次数 fno=fscanf(fid,'%d',1); %读入输出数据文件名 x=fscanf(fid,'%f',inf); %读入时程数据存成列向量 status=fclose(fid); n=length(x); t=(0:1/sf : (n-1)/sf); b=ones(1,l); a=[b*x(1),x,b*x(n)]; b=a; for k=1:m for j=l+1:n-1 b(j)=mean(a(j-1:j+1)); end a=b; end y=x(1:n)-a(l+1:n+l); plot(t,x,':',t,y,t,a(l+1:n+l),'-.'); xlabel('时间(s)'); ylabel('位移 mm'); legend('输入','输出','趋势'); grid on; fid=fopen(fno,'w'); for k=1:n fprintf(fid,'%f%f\n',t(k),y(k)); end status=fclose(fid); [ Last edited by heio on 2012-8-9 at 10:29 ] |
回复此楼» 猜你喜欢
谈谈两天一夜的“延安行”
已经有15人回复
-
EST投稿状态问题
已经有6人回复
-
职称评审没过，求安慰
已经有15人回复
-
垃圾破二本职称评审标准
已经有11人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有16人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
小新爱小白
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1365.5
- 红花: 1
- 帖子: 80
- 在线: 51.8小时
- 虫号: 1919956
- 注册: 2012-08-01
- 性别: GG
- 专业: 工程管理
【答案】应助回帖
★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2012-08-09 18:11:02
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2012-08-09 18:11:02
heio
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1363.9
- 散金: 41
- 帖子: 167
- 在线: 66.1小时
- 虫号: 1334666
- 注册: 2011-06-30
- 性别: GG
- 专业: 光学
送鲜花一朵小新爱小白
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1365.5
- 红花: 1
- 帖子: 80
- 在线: 51.8小时
- 虫号: 1919956
- 注册: 2012-08-01
- 性别: GG
- 专业: 工程管理