版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

yhyffrshuh

木蟲 (正式寫手)

流浪者一枚

應助: 2 (幼兒園)
金幣: 1926.8
散金: 745
紅花: 9
帖子: 807
在線: 492.1小時
蟲號: 427865
注冊: 2007-07-29
性別: MM
專業(yè): 基礎物理學

[求助] 兩個for循環(huán)出錯

問題：設置N為3，F(xiàn)為5，A97中的每一行（一個樣本）的到的平滑數(shù)據(jù)、一階導數(shù)和二階導數(shù)均是是1*104的矩陣；有97個樣本，得到的數(shù)據(jù)應該是97*104的矩陣，但按下面的程序運行得到的平滑數(shù)據(jù)和一階導數(shù)，二階導數(shù)為 1*9897，這顯然不對啊

程序代碼如下：
clear,clc
load A97 %導入光譜,第一列為編號，最后一列為濃度值
dx=1;
x=215:321; %波長范圍是215nm-321nm
a=length(x); %x的數(shù)據(jù)個數(shù)
k=0:96; % k為樣本編號
SG0=[];
SG1=[];
SG2=[];
N=input('請輸入擬合次數(shù)');
F=input('設置窗口參數(shù)(奇數(shù))');
for k=k+1
y=A97(k,2:108); %第k個樣本 215nm-321nm對應的吸光度
[b,g]=sgolay(N,F);
Halfwin=((F+1)/2)-1;
for n=(F+1)/2:a-(F+1)/2,
%SG平滑
SG0(n)=dot(g(:,1),y(n-Halfwin:n+Halfwin));
SG0=[SG0,SG0(n)];
%一階差分求導
SG1(n)=dot(g(:,2),y(n-Halfwin:n+Halfwin));
SG1=[SG1,SG1(n)];
%二階差分求導
SG2(n)=2*dot(g(:,3)',y(n-Halfwin:n+Halfwin))';
SG2=[SG2,SG2(n)];
end

SG1=SG1/dx;
SG2=SG2/(dx*dx);

subplot(3,1,1);
plot(SG0)
legend('S-G Smooth');

subplot(3,1,2);
plot(SG1)
legend('S-G Smooth 1st derivative');

subplot(3,1,3);
plot(SG2)
legend('S-G Smooth 2st derivativa');
end

回復此樓

» 猜你喜歡

288求調(diào)劑已經(jīng)有5人回復
281求調(diào)劑（0805）已經(jīng)有19人回復
材料考研調(diào)劑已經(jīng)有3人回復
346求調(diào)劑[0856] 已經(jīng)有6人回復
294求調(diào)劑材料與化工專碩已經(jīng)有9人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有12人回復
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑已經(jīng)有5人回復
085601材料工程專碩求調(diào)劑已經(jīng)有7人回復
321求調(diào)劑已經(jīng)有6人回復
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復

Allhumanwisdomissummedupintwowords:waitandhope.

1樓 2012-04-07 12:41:50

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

做消費行為研究，兩個問題得分相乘得出來的分數(shù)是不是連續(xù)變量？已經(jīng)有4人回復
intel fortran局部變量自動更新已經(jīng)有19人回復
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧已經(jīng)有105人回復
【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離已經(jīng)有34人回復
【求助】在進行fluent模擬求解過程中，突然提醒出錯，然后進入死循環(huán)，怎么回事？已經(jīng)有5人回復

點擊這里搜索更多相關資源

科研從小木蟲開始，人人為我，我為人人

相關版塊跳轉(zhuǎn) 我要訂閱樓主 yhyffrshuh 的主題更新