24小时热门版块排行榜

>论坛更新日志 (4579)
>文献求助 (470)
>虫友互识 (376)
>导师招生 (311)
>考博 (185)
>招聘信息布告栏 (164)
>硕博家园 (157)
>休闲灌水 (115)
>论文道贺祈福 (114)
>博后之家 (100)
>考研 (88)
>论文投稿 (78)
>基金申请 (67)
>教师之家 (67)
>公派出国 (62)
>找工作 (55)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wang1273

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 151.4
散金: 3
帖子: 600
在线: 137.4小时
虫号: 928886
注册: 2009-12-15
专业: 分析仪器与试剂

[求助] 求助用matlab编一个小程序

我想用matlab求光谱矩阵的极大线性无关向量组，请各位帮忙编一个程序，谢谢！

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

matlab的plot画图（如何画出与取值数目不一样的标记符）已经有4人回复
MATLAB画图显示颜色与打印结果颜色显示有一部分不一致已经有6人回复
MATLAB绘图-直线不平滑已经有6人回复
请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！已经有9人回复
一段matlab的小程序，不明白它的意思已经有5人回复
这个c++小程序编译的错误很奇怪~~ 已经有5人回复
用matlab来实现分类汇总已经有3人回复
自己写的一个小程序，出现问题不知道怎么解决已经有9人回复
有个概率公式，有积分。用啥软件可以实现呢？已经有7人回复
matlab程序问题，总是出错，麻烦大家看看问题在那里？已经有6人回复
请帮忙编写一个perl程序处理文本已经有26人回复
用matlab编一个计算稳态概率的小程序已经有4人回复
matlab编一个小程序解三角形已经有4人回复
c语言求助，一个小程序详细讲解已经有13人回复
请帮忙用matlab求矩阵特征值已经有8人回复
求助大侠帮忙，有关matlab模拟数据已经有29人回复
大家帮我看一个fortran的程序，我总是计算不出正确的结果已经有20人回复
如何用fortran编写高斯白噪声程序已经有6人回复
【求助】如何用Matlab求联合概率曲线已经有3人回复
【求助】请教一个matlab程序，谢谢帮忙！！已经有15人回复
【求助】请问一下fortran的可视化编程已经有7人回复
【求助】双倒数法求米氏常数中的稳态电流是CV测量得到的么？已经有5人回复
【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码已经有8人回复

1楼2011-12-28 10:43:36

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

zhfzh

木虫 (正式写手)

应助: 19 (小学生)
金币: 6620.5
散金: 200
红花: 4
帖子: 685
在线: 665小时
虫号: 715595
注册: 2009-03-05
专业: 原子和分子物理

【答案】应助回帖

★ ★ ★
wang1273(金币+5): ★★★★★最佳答案 2011-12-30 14:58:38
余泽成(金币+3): 谢谢参与应助！ 2011-12-30 20:47:50

%求矩阵的最大无关组
%数据文件格式说明
%第一行行数（n）  列数（m）
%以后矩阵元素
%************代码段****************************************************
clc    %清屏
clear all %清除变量
%***************读入数据***************************************************
%打开文件，读取数据
[Filename,Pathname]=uigetfile('*.*','请打开数据文件');
fid=fopen([Pathname Filename],'r');
[data,count]=fscanf(fid,'%d',inf);
fclose(fid);
%**********************判断数据是否正确*********************************
[rdata,ldata]=size(data);
row=data(1,1);
col=data(2,1);
if ~(rdata)==0
if  row*col+2==rdata
      fprintf('\n\n^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^\n');
      fprintf('          @数据导入成功！@\n');
      fprintf('       矩阵的维数为%dx%d\n\n',row,col);
      fprintf('^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^ ^_^\n\n');
else
      fprintf('\n\nvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n');
      fprintf('          @数据文件格式错误！@\n');
      fprintf('vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n\n');
      break;
end
else
fprintf('\nvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n');
fprintf('             @文件打开失败!@\n');
fprintf('vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n\n');
break;
end
%*************处理矩阵*************************************************
%把数据放入mdata
mdata=zeros(row,col);
start=3;
for i=1:row
mdata(i,1:col)=data(start:start+col-1,1)';
start=start+col;
end
%求极大无关组
d_mdata=rref(mdata);
c=zeros(col,1);
for i=1:row
for j=1:col
      if d_mdata(i,j)~=0
         c(j,1)=1;
         break
      end
end
end
d_data=[];
k=1;
for i=1:col
if c(i,1)==1
      d_data(:,k)=mdata(:,i);
      k=k+1;
end
end
fprintf('行最简形为：\n');
d_mdata
fprintf('一个极大无关组为：\n');
d_data

睡前写了下，不会上传文件，自己复制吧