1樓:匿名使用者
a = 0;
for i = 1 : 10
a = a + i ^ 2;
enda = 0;
i = 1;
while i <= 10
a = a + i ^ 2;
i = i + 1;end
2樓:
function ans = for(n)sum = 0;
for i = 1:n
sum = sum+i^i;
i = i+1;
endans = sum;
function ans = while(n)sum = 0;
i = 1;
while i<=n
sum = sum+i^i;
i = i+1;
endans = sum;
根據matlab的已學知識,寫一個matlab函式m檔案,具體要求如下:
3樓:白雲裳愛戰天風
%matlab聯絡之——利用最小二乘法擬合函式
clear;clc
%原始離散點輸入
x=[-1.0 -0.5 0.0 0.5 1.0 1.5 2.0];
y=[-4.447 -0.452 0.551 0.048 -0.447 0.549 4.552];
%原始離散點作圖
figure(1)
plot(x,y,'o','linewidth',2)
%擬合階數定義
n=3+1;
p=ones(n,length(x));
%最小二乘法擬合
p(1,1:length(x))=ones(1,length(x));
pp(1,n)=[1];
for h=2:n
p(h,:)=zeros(1,length(x));
a(h-1)=sum((x.*p(h-1,:)).*p(h-1,:))/sum(p(h-1,:).^2);
if h-2>0
a(h-2)=sum(p(h-1,:).*p(h-1,:))/sum(p(h-2,:).^2);
else
p(h,:)=p(h,:)+(x-a(h-1)).*p(h-1,:);
pp(h,n-1:n)=[1 -a(h-1)];
continue
endp(h,:)=p(h,:)+(x-a(h-1)).*p(h-1,:)-a(h-2)*p(h-2,:);
pp(h,n-h+1:n)=conv([1 -a(h-1)],pp(h-1,n-h+2:n))-a(h-2)*[zeros(1,2) pp(h-2,n-h+3:n)];
end%相關評價引數的求取
a=zeros(n,n);
for h=1:n
a(h,h)=sum(p(h,:).^2);
endb=ones(n,1);
for h=1:n
b(h)=p(h,:)*y';
endfor h=1:n
a(h)=b(h)/a(h,h);
endp_val=0;
for h=1:n
p_val=p_val+a(h)*p(h,:);
endp_a=zeros(1,n);
for h=1:n
p_a=p_a+a(h)*pp(h,:);
ende=sum(y.^2)-y*p_val'
%擬合函式畫圖,同離散點一起顯示
x=linspace(x(1),x(length(x)),10*length(x));
p_val=polyval(p_a,x);
figure(1)
hold on
plot(x,p_val,'r','linewidth',3)
legend('原始離散點','最小二乘法擬合曲線')
xlabel('x'),ylabel('y')
title('最小二乘法擬合示範結果')
matlab 怎樣執行m檔案,matlab怎麼呼叫m檔案?
執行的方式有2種,一是在命令中心輸入.m檔案的名字 如果有引數,就需要給出引數 二是在.m檔案編輯環境中直接執行,一般是在debug選單項中選run或者直接按f5 執行的結果也是在命令中心顯示 要是是畫圖的話,就是figure視窗了 直接用dos函式執行就可以了。舉例說,如果批處理檔名 test b...
怎樣將將寫好的m檔案加到matlab庫函式中
在file選單下面有 set path選項,將你的函式所在資料夾加入到對應matlab搜尋路徑下就行了 在一個m檔案中寫多個函 copy數,每個函式都以function開始,結束可以帶end,也可以不帶 這種情況下,出現另一個function或者到達檔案結尾表明前一個function結束 第一個函式...
怎麼用matlab程式設計讓T 1 2
1 開啟matlab軟體。2 a zeros n 該函式產生n n的全零矩陣。例如生成4 4的全零矩陣如下。3 a zeros m,n 該函式生成m n的全零矩陣。例如生成2 3的矩陣。4 a zeros size b 該函式產生和矩陣b維數相同的矩陣。例如生成與b維數相同的矩陣c。5 a 1.ze...